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ABSTRACT 



A Search for the Singlet-P State h c {l l P x ) of 
Charmonium Formed in Proton- Antiproton 
Annihilations at Fermilab Experiment E835p 



David Noah Joffc 



We present the results of a search for the spin-singlet P-wave state /i c (l 1 Pi) of 
charmonium formed through proton-antiproton annihilation at Fermilab experiment 
E835. The decay channels which were studied were pp — > J ftp + X — > e + e~ + X, 
pp — > J/i/j + 7T° — > e + e~ + 77, pp — > J/ip + 7r°7r° — > e + e~ + 47, and the neutral 
channel — > 7^7 — > (77)7. The decay — > J/^7 — ► e + e~7, into which : Px decay 
is forbidden by C-parity conservation, was also examined for comparison. 

The 90% confidence upper limits for the decay channels studied in the mass 
range 3525.1-3527.3 MeV for a 1 P\ resonance with a presumed width of 1.0 MeV 
were determined to be B(pp P x ) x B( 1 P 1 -> J/if> + X) < 1.8 x 10~ 7 , B(pp -> 
1 P 1 ) x B( 1 P 1 -> J/^ + 7T°) < 1.2 x 10- 7 , and B(pp Pi) x P( 1 P 1 -> J/^7) < 
1.0 x 10~ 7 . No evidence for a 1 P\ enhancement was observed in either of the two 
additional reactions studied; pp — > J/ip + 7r°7r° — > e + e~+47 and pp — > 7^7 — > (77)7. 
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Chapter 1 
Introduction 



In the Standard Model of elementary particle physics, all matter is described as 
being composed of quarks and leptons, and the forces between them are described 
as being mediated by four gauge bosons; photons mediate the electromagnetic force, 
gluons mediate the strong hadronic or nuclear force, and W ± and Z° bosons mediate 
the weak force. The quarks participate in all three forces, the charged leptons in the 
electromagnetic and weak forces, while the neutral leptons (neutrinos) participate 
only in the weak force. In this description of the natural world, these constituents are 
themselves structureless or pointlike; they are considered to be truly fundamental. 
Although the idea of explaining the world in terms of a series of ultimate constituents 
goes back at least as far as the ancient Greeks with the atomic theories of Democritus 
and Epicurus, the Standard Model is by far the most successful such description, 
effectively including all basic physical phenomena with the exception of gravity. The 
quarks and leptons which make up the basis of the Standard Model are shown in 
Table 1.1. 

In the Standard Model, all strongly interacting, or hadronic, matter (and thus 

1 



2 



Table 1.1: Components of the Standard Model: the three generations of leptons and 
quarks, and the gauge bosons. The masses for quarks are the so-called 'constituent 
quark' masses. [T] 



Generation 


Quarks 




Lept 


ons 






Charj 


;e Mass 




Charge 


Mass 


I 


d 


-1/3 


~ 300 MeV 


e 


-1 


0.511 MeV 




u 


+2/3 


~ 300 MeV 







< 15 eV 


II 


s 


-1/3 


~ 500 MeV 


^~ 


-1 


105.66 MeV 




c 


+2/3 


~ 1500 MeV 







< 170 keV 


III 


b 


-1/3 


~ 4700 MeV 




-1 


1777.05 MeV 




t 


+2/3 


173.8 GeV 


V T 





< 18.2 MeV 


Gauge 










w ± 


±1 


80.42 GeV 


Bosons 


9 








z° 





91.19 GeV 



the vast majority of the matter of which our world is made) is composed of quarks, 
which interact via the exchange of gluons. The theory governing this quark-gluon 
interaction is known as Quantum Chromodynamics, or QCD. Many of the basic 
concepts of QCD date back to the 1960's, and the term 'quark' itself was coined by 
Gell-Mann in 1964 [2] in his discription of the 'eight-fold way' of understanding the 
symmetries among hadrons. Three quarks, 'up', 'down' and 'strange' were initially 
proposed to explain all hadronic states then observed. In this three quark model, 
the lowest-mass hadrons are explained as bound states of three quarks (baryons), 
three antiquarks (antibaryons), or a quark- ant iquark pair (mesons). The lowest- 
mass baryons (and antibaryons) are arranged into an octet of J p = | + states and 



3 



A" A A + A* 




a 

Figure 1.1: SU(3) classification of light quarks, and of hadrons into a J p = 0~ 
meson nonet, a spin- 1/2 baryon octet, and a spin-3/2 baryon decuplet. 

a decuplet of J p = | + states, where J and P are the spin and parity quantum 
numbers (see Figure 1.1). The lowest mass mesons are arranged into two nonets 
with J p = 0~ and J p = 1~. 

In 1964, Greenberg proposed the term 'color' to represent the strong interaction 
charge' [3]. Three such 'colors' were necessary to explain the existance of the corner 
members of the baryon decuplets, each of which contains three identical quarks in 
relative s-states. The Pauli exclusion principle forbids this unless the three quarks 
are made distinguishable by assigning them a new quantum number, color, making 
them antisymmetric under color exchange. Although the quarks themselves were 
assigned three colors, termed red, green, and blue, only color neutral (or white) 
objects were allowed to exist freely in nature. The interaction between quarks was 
mediated by massless vector bosons called gluons, which carry both color and anti- 
color. This color interaction between quarks and gluons forms the basis of the theory 
of Quantum Chromodynamics (QCD) which will be discussed in detail in Chapter 
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2. 



Although three quarks were sufficient to explain all hadrons which had been 
observed until the 1960's, a fourth quark was soon proposed. Bjorken and Drell first 
proposed its existance in order to allow the quarks to be grouped into two doublets 
in analogy to the two lepton doublets, (e,u e ) and (//, z^), then known to exist [I]. 
They gave it the name 'charm' quark. 

Further evidence for the existance of the fourth quark came when Cabibbo pro- 
posed [3] that the quark states actually participating in the weak interaction were 
the u and an admixture, d c = d cos 9 C + s sin 9 C , of the physical quarks d and 
s. This model successfully explained the experimentally-observed suppression of 
the strangeness-changing ( AS* = 1 ) semileptonic weak decays with respect to the 
strangeness-conserving ( AS = ) decays. The experimental value of the Cabibbo 
angle 9 C was found to be 0.25 rad. One consequence of the Cabibbo theory, however, 
was the existence of neutral currents with AS = 1 , which were not observed in 
nature. A famous paper by Glashow, Iliopoulos and MaianijB] in 1970 solved this 
problem and provided strong theoretical support for the existence of a fourth quark, 
named charm. The authors proposed that the charm quark c would participate in 
the weak interaction in a doublet with a state, s c = s cos 9 C — d sin 9 C , orthogonal 
to the d c state. Hence the two quark doublet for weak interaction would be: 



' ' \0.97d + 0.25sJ \0.97s -0.25dJ ' 
Thus s c s c + d c d c = ss + dd independently of the mixing angle 9 leading to a vanishing 






c 
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term for the neutral currents with AS = 1 , in agreement with the experimental 
observations. 

Despite these theoretical successes, the quark model, and thus the core of the 
Standard Model itself, really only became universally accepted a decade later with 
the discovery of the J/ip in 1974. 

This discovery was made essentially simultaneously by two separate groups work- 
ing on opposite sides of the U.S. The east coast group, led by Ting at the Brookhaven 
National Laboratory (BNL), reported the observation of a particle which they called 
J, in the e + e~ invariant mass in the reaction p + Be — > e + e~ + X at 28 GeV |10j . 
They reported a mass of 3.1 GeV, and a width "consistent with zero". The BNL 
results can be seen in Figure 1.2. The west coast group, led by Richter at the 
Stanford Linear Accelerator (SLAC), reported a resonance, which they called if), in 
the reactions e + e~ — > e + e~ , , hadrons jTT]. The SLAC results can be seen in 
Figure 1.3. They reported a mass of 3.105(3) GeV, and a width of < 1.3MeV. Both 
discoveries were published in the December 2, 1974 issue of the Physical Review 
Letters, and the resonance eventually came to be known as the J/ip. The discov- 
ery was quickly confirmed by e + e~ annihilation experiments at Frascati ^2] and 
DESY HB]. 

The discovery of the J/tp prompted a spate of theoretical papers within weeks 
of the announcements, the most important of which were that of Appelquist and 
Politzer ^1] and De Rujula and Glashow ^Hj which proposed the interpretation of 
the J /ip as the bound state of a charm quark and an anticharm quark. 
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With the discovery of the J/ip, and its interpretation as a charm-anticharm cc 
bound state, the second family of the quark sector of the standard model became 
experimentally established. Other cc states were soon found, starting with the ip', 
the first radial excitation of the J/ip, which was discovered at SLAC only days after 
the J/ijj was observed [T?j . 




2.6 2,75 3.0 3,25 3,5 

M,,,-, GeV 

Figure 1.2: Results of Aubert et al. (1974) indicating the narrow resonance J/i/j 
in the invariant mass distribution of e + e~ pairs produced in inclusive reactions of 
protons with a beryllium target. The experiment was carried out with the 28-GeV 
AGS at the Brookhaven National Laboratory. [IS] 
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Figure 1.3: Results of Augustin et al. (1974) showing the observation of the J/ip 
resonance of mass 3.1 GeV, produced in e + e~ annihilation at the SPEAR storage 
ring, SLAC. [II] [E| 



s 



Further evidence that the newly found J/ip and if)' states were in fact bound 
states of cc, was that their widths were soon determined to be very small, < 100 keV 
and < 300 keV respectively. Most strong interaction resonances with smaller masses 
were known to be much wider, as large as a few hundred MeV, i.e., three orders of 
magnitude larger than those of the newly observed states. This made it difficult to 
explain J/ip and ip' m terms of the u, d, and s quarks. By interpreting the J /ip and 
ip' as cc states, and appealing to the Okubo-Zweig-Iizuka (OZI) rule j7j [S] 0, the 
narrowness of the states could be easily explained. 

The OZI rule states that processes which can only be described by diagrams that 
contain disconnected lines (i.e. no quark flow) between the initial and final states 
should be strongly suppresed as compared to diagrams which contained connected 
lines (see Figure 1.4(a,b)) 

The OZI rule can be explained intuitively by noting that diagrams containing 
disconnected lines require the emission of "hard" or highly energetic gluons, which 
must carry the full four-momentum of the annihilating quark-antiquark pair (Fig. 
1.4(a)). These gluons are much less likely to be produced than the "soft" gluons 
emmitted by a quark that continues to exist in the final state, as in the connected- 
line diagram (Fig. 1.4(b)). Thus, cc states which do not have enough mass to 
decay via the lowest-energy connected diagram (D + D~), must necessarily be narrow, 
thus expaining the small widths of the J/ip and ip' resonances. After their spins 
had been determined by studying interference effects and angular distributions of 
decay products, the J / ip and ip' were assigned the quantum numbers of the photon: 
J PC = 1 , where J and P are the spin and parity and C is the charge conjugation 
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Figure 1.4: Examples of disconnected and connected diagrams. Top (a): the discon- 
nected diagram of the decay J/ip —>■ 7r _ 7r°7r + . Bottom (b): the connected diagram 
of the decay ^"(3770) -> D+D~. 

parity. 

After the discovery of the J/ip and the ip', the SLAC group continued to make 
important new observations over the next several years with succesively improved 
detectors; Mark I, Mark II, Mark III, and finally Crystal Ball. After the J ftp and ip', 
the SLAC group observed the Xcj resonances (the ++ , 1 ++ , and 2 ++ bound states 
of cc), followed by the h cc ground state, the r] c , named in analogy with the light 
quark rj. Figure 1.5 shows the Crystal Ball observation of the Xcj resonances, i] c , 
and rj' c (later shown not to be true) [TH]. The DESY, Orsay, and Frascati groups also 
made imporant contributions to the spectroscopy of cc states, which became known 



10 



collectively as charmonium. The spectrum of charmonium states is shown in Figure 
1.6. Charmonium states are labeled using the spectroscopic notation n 2s+1 Lj, where 
n is the number of nodes in the radial excitation plus one, s is the combined spin 
of the two quarks, L is the orbital angular momentum, and J is the total angular 
momentum. Parity and charge conjugation, as in any quark- ant iquark state, are 
given by P = (— 1) L+1 and C = (— 1) L+S respectively. A more abbreviated notation 
is to characterize the states just by their J PC . 




Figure 1.5: Inclusive photon spectrum from ■0' decay, as measured by Crystal Ball 
(SLAC). [HI 
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Figure 1.6: Spectrum of charmonium resonances. 
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The last charmonium state to be discovered was the T]' c , the radial excitation 
of the r] c ground state, whose existance was firmly established only in 2003, by 
Belle by our own group at CLEO|20], and BaBar |2*T] . The CLEO results are 
shown in Figure 1.7. 
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Figure 1.7: K^K ±r K T invariant mass in the reaction 77 — > K^K^ir^ from (a) CLEO 
II/II.V and (b) CLEO III data [H|. The rj' c is the enhancement at 3.643 GeV. 
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Figure 1.8: Formation of a vector charmonium resonance from e + e annihilation. 



All studies of charmonium states described so far utilized e + e~ annihilation, in 
which only the cc vector states ( J PC = 1 ) are directly formed via the intermediate 
photon with J PC = 1 (see Figure 1.8). All other states are only reached by decays, 
mostly radiative, from these vector states. 

A major departure from this technique came with the R704 experiment at the 
ISR at CERN [22], which demonstrated that high resolution charmonium spec- 
troscopy could be studied by charmonium formation in proton-antiproton annihi- 
lation (see Figure 1.9). This technique has the advantage of being able to directly 
populate cc mesons of all J PC via two and three gluon annihilations leading to C 
= +1 and C = — 1 states respectively, and was used by Fermilab experiments E760 
and E835 in their 1990, 1997, and 2000 runs. 
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Figure 1.9: Formation of a charmonium resonance by two and three gluon pp anni- 
hilation. 

Because it was a pp annihilation experiment which could directly populate the 
non- vector charmonium states, Fermilab E760 was in a unique position to study the 
Xcj states in greater detail, as well as to conduct a search for the h c or l 1 Pi(l +_ ) 
state of charmonium. 

The 1 P\ resonance, which is the singlet-P partner of the Xcj resonances, is the 
final bound state of cc whose observation remains unconfirmed. The Crystal Ball 
experiment at SLAC failed to find it in the reaction e + e~ — > ip' — > x Pitt , 1 P\ — > 
'yrjc |23j . During its 1990 run, E760 claimed observation of the 1 Pi in the reaction 
pp — > x Pi — > J/tpTi — > e + e~ based on 15.9 pb -1 of data taken. It claimed a 1 P\ 
mass of 3526.2 ±0.15 ±0.2 MeV and a width of V < 1.1 MeV (see Figure 1.10) (23- 

In 1997, an attempt was made to confirm this observation in the Fermilab experi- 
ment E835, the successor experiment to E760, using 38.9 pb _1 of data. This attempt 
was inconclusive, however, due to instabilities in the antiproton beam energy. Dur- 
ing the year 2000 run of E835 (E835p), 50.5 pb" 1 of data was taken dedicated to 
confirming the E760 observation of l P\. The analysis of these data in a variety of 
possible 1 P 1 decay channels forms the topic of this dissertation. 
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Figure 1.10: The l P\ resonance as measured by Fermilab experiment E760. [21] 



Chapter 2 of this thesis consists of a theoretical discussion of charmonium spec- 
troscopy, in terms of the quark model and QCD, and an examination of theoretical 
predictions regarding the 1 P%. In Chapter 3, the experimental set-up of E835 is de- 
scribed in detail. In Chapter 4, the data selection and results of the X P\ search 
is discussed for reactions leading to final states containing J/ip: pp — ^ P\ — > 
J/i/j + X — > e + e~ + X, pp — Pi — > J /if) + 7T° — > e + e~ + 77, and pp — s- 1 Pi — > 
J/ip + 7r°7r° — > e + e~ + 47. Chapter 4 also describes the study of the reaction 
pp — > J /if) + 7 — > e + e~7, which is forbidden to occur via 1 Pi by C-parity conserva- 
tion, and which was studied as a control. Chapter 5 describes our search for /i c ( 1 P 1 ) 
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in the reaction pp — s- 1 Pi — > rj c + 7 — > (77) + 7 in the E835p data. Data tables 
for the final event selection in the — ^ Pi — > J ftp + 7r° — > e + e~ + 77 reaction 
are given in Appendix A, and software used in determining luminosity during data 
taking is given in Appendix B. 



Chapter 2 



Theory 



2.1 Quantum Chromodynamics 

In the Standard Model of particle physics, Quantum Chromodynamics, or QCD, 
is the established theory of strong interactions. QCD is the non-Abelian gauge 
theory of quarks and gluons, which carry a strong interaction charge called 'color'. 
The theory is non-Abelian because, unlike Quantum Electrodynamics (QED), the 
established theory of electromagnetic interactions, the gauge boson mediating the 
field itself carries the fundamental charge. Thus, while photons have no electric 
charge, and do not directly couple to each other, gluons, the QCD gauge bosons 
carry both a color and and anti-color, and do directly interact. The basic QCD 
interaction is invariant under the interchange of the different colors. The gluons are 
postulated to belong to an octet representation of the symmetry group SU(3), and 
the color states of the gluons in this octet can be expressed as: 



- _ _ _ _ - rr-bb rr + bb- 2gg 
rb,rg,bg,br,gr,gb,—^=—, -j= (2.1) 
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This non-Abelian field, combined with the much larger coupling constants for 
QCD, mean that although the formalism is pattered after the Abelian QED theory, 
in practice, QED and QCD are quite different. As can be seen in Figure 2.1, the 
color-charge of the gluons means that Feynmann diagrams with three and four gluon 
vertices are allowed, as well as the two gluon-quark vertices which would be expected 
in analogy to QED. Furthermore, the strong coupling constant a s , which governs the 
quark-gluon interaction, is much larger than its electromagnetic counterpart a em , 
and the value of a s tends to infinity for small values of q 2 (large distances). Because 
of these complications, QCD, unlike QED, is a theory in which it is impossible, 
to solve problems analytically, and thus requires extremely cumbersome numerical 
calculations. QCD based computer models of varying degrees of sophistication are 
thus used to generate theoretical predictions with which the experimental data may 
be compared. A brief description of QCD and some of the more popular calculation 
techniques are given in this chapter, for further details a very good synopsis of QCD 
is given by Wilczek j2S], and a summary of QCD calculational techniques as applied 
to charmonium is given by Seth [26J. In the following we liberally borrow with 
permission from the review by Seth [2*o] . 

Charmonium is a good system to test the assumptions of QCD. Unlike the case 
of light-quark hardons, for charmonium the value of a s is sufficiently small ~ 0.3 to 
make perturburtive calculations possible. Furthermore, the relatively small binding 
energy compared to the rest mass of its constituents allow cc states to be described 
non-relativistically (with v 2 /c 2 ~ 0.25). The fact that the charmonium resonances 
are eigenstates of J PC produces symmetry conserving simplifications. Finally, the 
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Figure 2.1: The three basic QCD vertices. 

bound cc states are well separated in energy and narrow in width, as opposed to 
the light-quark resonances which have large, often overlapping, widths, creating a 
complicated and difficult spectroscopy. 

The study of charmonium also has experimental advantages over the heavier 
quarks; the top quark decays too quickly to form any QCD bound states, and the 
production cross sections for bb, or bottomonium, (particularly in pp) are much 
smaller than that for charmonium, making direct production of anything other than 
the J PC = 1 states all but impossible. 

By making precise measurements of the masses, widths, and branching ratios 
charmonium states, important information about the dynamics of the strong in- 
teraction may be extracted. For instance, by comparing the hadronic and electro- 
magnetic branching ratios of appropriately chosen charmonium states, an estimate 
of the strong coupling constant a s can be derived [21]. Unknown quantities, such 
as the squared absolute value of the wave function at the origin |^(0) | 2 , or poorly 
measured quantities, such as branching ratios between the resonance and the initial 
state, may often cancel in the ratio, leaving a s as the only unknown. For example, 
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in the case of r) c , r)' c , Xco, and Xc2, this cancelation occurs when one compares the 
branching ratio into two photons and the branching ratio into two gluons. Different 
theoretical models may also provide predictions for the radiative partial widths of 
charmonium states which may be compared to the experimental results. Examples 
of this include the electric dipole transitions of the three Xcj states to J ftp, and the 
magnetic dipole transition of J/ip to i] c , and ip' to r)' c . 



The most fundamental expression of the theory of Quantum Chromodynamics be- 
gins with the QCD Lagrangian itself. This Lagrangian, which describes the inter- 
action of quarks and gluons is: 



Here, G$ is the gluon field strength tensor, with A$ as the eight gluon fields 
(a = 1, 2, 8), q are the quark fields of six different flavors (u, d, s, c, b, t), each of 
three colors, m q are the quark masses, and V is the covariant derivative. In turn, 
the gluon field strength tensor and the covariant derivative are described by: 



2.1.1 The QCD Lagrangian 




(2.2) 



G« = d,A^ - d v Af + g s f abc A^Ai c 



(2.3) 



and 



V = d ll -igAW\W/2 



(2.4) 
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where A, the Gell-Mann matrices, and f abc are the generators and structure con- 
stants, respectively, of SU(3).The gauge constant g s , which determines the coupling 
between the quark and gluon fields, is related to the scale dependent (or "running") 
effective coupling constant of strong interaction by the relation: 

a, = I (2.5) 

This coupling constant, a s depends on the choice of the renormalization scale, 
At the energy scale /i, in the lowest order, 

" s(/i) = [11 - (2/3)71/)] x Zn(/x/A) (2 ' 6) 

where Uf is the number of quark flavors with mass less than /i, and A is the 
QCD scale parameter, which depends on the number of quark flavors. At the mass 
of the Z° boson, a s {Z°) is measured to be 0.1172 ± 0.0020 |T . 

As can be seen from the above expression for a s , QCD incorporates a unique 
property known as 'asymptotic freedom'; a s decreases as one goes to higher energies 
or smaller distances, becoming zero at fi = oo. It is this feature of QCD which 
allows one to use perturbative methods at high energies (and small distances). At 
low energies a s becomes large, and the use of perturbative QCD becomes highly 
questionable for the light (u, d, s) quarks. 

The non-linear terms in the field strength tensor makes QCD impossible to 
solve analytically, and difficult to work with even in lattice calculations. However, 
certain symmetries of the theory can be seen from the Lagrangian itself [2*5] . 

1. The discrete symmetries P, C, and T are preserved by the strong interaction. 
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2. Quark flavor is conserved by the strong interaction, as are the quantities 
which can be derived from quark flavor such as isospin, electric charge, and baryon 
number. 

3. Approximate chiral and flavor symmetries are preserved for the light u,d,s 
quarks, or in the very high energy regime, when all quark masses are comparatively 
negligible. 

It is also worth noting that neither QCD nor anything else in the Standard Model 
predicts either the number of quark families, the scale parameter A, or the quark 
masses. These must be added to the theory as input parameters. 

Using the QCD Lagrangian, rules for calculating Feynman diagrams for strong 
interaction processes may be constructed. Since this method of calculation treats 
interactions essentially as perturbations about the field-free ground state, they can 
be used to obtain predictions for reactions only at very high momenta, where the 
strong coupling constant becomes small. At small momenta the strong coupling 
constant becomes too large, for these perturbative methods (pQCD) to remain valid, 
and one is then forced to use numerical solutions such as Lattice Gauge Calculations. 

2.1.2 Potential Models 

One technique used in calculations of hadronic resonances has been to replace the 
non-Abelian guage field theory of QCD by a non-relativistic potential model. De- 
spite the fact that many have questioned the very existence of a static potential in the 
presence of an active gluon condensate [IHj, potential model calculations have been 
very successful in predicting the masses of bound states of both charmonium and 
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bottomonium. Non-relativistic models for charmonium are possible only because of 
the relatively large mass of the charm quark, which have v 2 / c 2 ~ 0.25 in their bound 
states, as opposed to v 2 /c 2 ~ 0.8 for the light quark mesons. Using corrections for 
relativistic, channel coupling, and radiative effects, the success of potential models 
for theoretical predictions even extends to some decays of charmonium states as well 
as their masses and widths. 

Spin-Independent Potentials 

As early as 1975, Appelquist and Politzer recognized that the single gluon exchange 
between a charm quark and antiquark should give rise to a Coulombic potential pro- 
portional to 1/r at small distances ^3]. They coined the name "orthocharmonium" 
for the J/ip in analogy with the 3 Si orthopositronium, and went on to extend the 
analogy to predict the existence of 1 Sq " paracharmonium" . With the announcement 
of the discovery of if)' just a week after J /ip it became possible to flesh out the poten- 
tial picture, and Appelquist and Politzer were able to predict the complete spectrum 
of bound charmonium states based on a charmonium potential which was expected 
to be intermediate between a Coulombic and a harmonic oscillator potential. The 
spectrum of charmonium states and a comparison with the positronium spectrum 
is shown in Figure 2.2. 
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Figure 2.2: A schematic showing the similarities between positronium (top) and 
charmonium (bottom) spectra [2*%] . 
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After the initial predictions of Appelquist and Politzer, another early potential 
model which was examined was a purely linear potential, used by Harrington, Park, 
and Yildez in 1975 [29 . The most popular of these early potential models, however, 
is a combination of both the linear potential and Coulombic potentials. This model 
has become known as the Cornell potential. It was first proposed by Eichten et al. 
in 1975 [30J and combines a one-gluon exchange Coulombic component, proportional 
to 1/r, which is dominant at short distances, and an additional confinement term 
proportional to r, analogous to the "string tension" of a multigluon "flux tube", 
which is dominant at large distances. This linear confinement term was added in 
recognition of the fact that colored free quarks are never observed experimentally, 
but rather are confined to color singlet hadrons, either qq mesons, qqq baryons, or 
qqq antibaryons. Other objects are possible as well, such as glueballs (gg, ggg), 
hybrids (qqg), and multi-quark states (qqqq)- These objects are allowed in QCD as 
long as the total object is color neutral. 

The Cornell Potential, with its linear confinement term, may be written as : 

V(r) = -k/r + r/a 2 (2.7) 

where the constant k is often identified with 4/3a s , the 4/3 factor being a nu- 
merical artifact of color symmetry, and the constant I /a 2 is of order ~ 1 GeV/fm. 

Subsequently, many alternative formulations of the spin-independent potential 
have been suggested. One important modification was motivated by the observation 
that the energy of the radial excitiation for S states, A[M(2 3 S\) — M(l 3 5i)], is 
nearly the same for charmonium (589 MeV) and bottomonium (563 MeV). For this 
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property to be truly independent of quark mass one must use a logarithmic potential. 
Furthermore, the 25 — IP mass difference, which is zero in a purely Coulombic 
potential, is quite significant for charmonium, with M(2 3 5i)— < M(l 3 Pj) > = 
161 MeV. These considerations lead Quigg and Rosner [31] in 1977 to propose the 
following potential: 



Incorporating this logaritmic feature, several different forms of potentials, inter- 
polating between the Coulombic and confinement parts, have been proposed. One 
idea, first used by Richardson in 1979 [32|, is to construct the potential in momentum 
space taking into account the momentum dependence of a s as given by: 



This formulation has the advantage that it contains only one parameter, Aqqd, 
and in configuration space it varies as 1/r at small distances, and as r at large 
distances. Improvements on this approach were later made by Buchmuller [33j[34j. 
In Table 2.1 we list the configuration space representations of some of the many 
empirical potentials used to fit the charmonium spectra. 

All the different potential models obtain their parameters by fitting the spectra 
of the bound states of charmonium and bottomonium. The physical radius of these 
bound states range between r = 0.2 and 1.0 fermi, and in this region all potentials, 
including the purely empirical ones due to Quigg and Rosner [3T] , and Martin J3H] are 
essentially identical. These potential models share another important feature, which 





V{q 2 ) = -(4/3)[12tt/(33 - 2n f )][q 2 x ln(l + (q 2 /A 2 T 



l 



(2.9) 
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Author 


Potential 


Applequist and Politzer [Hj 


only Coulombic 


De Rujula and Glashow [T3] 


only Coulombic 


Kang and Schnitzer |H5] 


only Linear 


Harrington [29J 


only Linear 


Eichten jSHl 


-k/r + r/a 2 , k = (4/3)a s 


Quigg and Rosner [3T] 


A ■ ln(r/ro) 


Martin [H 


A + Br 0A 


Celemaster and Georgi [37J 


-(4/3)(l/87r6)(r/n(r/r ))" 1 - cre~ r/d + ar 


Celemaster and Henyey [37] 


-(4/3)(1/(2tt) 3 ) jd 3 k ■ exp{ik ■ r){Ana s {k 2 )/k 2 ) 


Richardson [32] 


-(4/3)(1/8tt 3 6) / d 3 k ■ exp{ik ■ f)/k 2 ln{l + k 2 /A 2 ) 



Table 2.1: Some QCD based spin-independent qq potentials. (After Lucha et 
al. [38 ). The constant 1/6 = 487r 2 /(33 — 2nf). The potentials used by Buchmiiller 
et al. j33| are essentially identical to that of Celemaster and Henyey, and are nearly 
equivalent to that of Richardson [32], since a s ~ l/4irbln(k 2 /A 2 ). 



2S 




R(fm) 

Figure 2.3: Schematic showing various model potentials as a function of quark 
separation distance R. Various cc and bb levels are indicated by the vertical lines [28J. 

is that even at these very small distances, the quark-antiquark pair does not reside in 
a purely Coulombic potential, and the confinement contribution is significant even 
in the lowest lying states. A schematic showing the differences between the various 
potential models is shown in Figure 2.3. 

Combining Lattice Guage Calculation techniques with these potential models, 
static spin-independent interquark potentials have been calculated on the lattice 
in both quenched and unquenched approximations, and excellent fit to the lattice 
potential for r > 0.3 fm is obtained with the parametrization of the Cornell potential, 
with k = 0.322, and a = 2.56 GeV -1 03]. There is little difference between the 
quenched and unquenched lattice results for distance up to r = 1.5 fm. 
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Spin-dependent Potentials 

Given the Coulombic nature of the short-range qq potential, and the spin depen- 
dence of the vector Coulomb potential, consisting of spin-orbit, spin-spin, and tensor 
components, the addition of spin-dependence to the emerging potential models was 
a natural step. The development of these spin-dependent potential models quickly 
followed the spin-independent versions, starting with Eichten et al. in 1975 J3EL, and 
Henrique, Kellett and Moorhouse in 1976 [35] . A full and systematic investigation 
of the spin dependent forces was done by Eichten and Feinberg j^U] in 1981 us- 
ing a gauge-invariant formalism. Their representation of the spin dependent forces 
was done with a nonrelativistic potential, with a short-range vector exchange and 
a long-range scalar exchange. In this model, the spin dependent potential may be 
expressed as: 

V SD = V'{r)S ■ L + 4a s /(3m 2 r 3 )S ■ L 

+ {327c/9mq 2 ){Sl ■ S2)6(r) 

+ (4a s /3mg 2 r 3 ) x (3{S1 ■ R){S2 ■ R) - SI ■ S2) 

where V'(r) is the derivative of the spin independent central potential. Using the 
Cornell potential for V(r), the first two terms combine to give the overall spin-orbit 
potential: 



V SL = (l/2mgV)(4a s /r + r/a 2 )S ■ L. (2.10) 

The other two terms can be identified as the spin-spin potential, and the tensor 
potential. 
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Precise measurements of the different resonance masses, or more particularly the 
differences between them, are a very effective way to test the spin-dependence of 
the different potential models. For instance, the tensor and spin-orbit interaction 
split the masses of the Xcj(l 3 Pj), J = 0, 1,2 states (fine splitting). The spin-spin 
force splits the vector and psuedoscalar states, and this is repsonsible for the mass 
difference between J/ip and t] c , and between ip f and v( c (hyperfine splitting). A 
measurement of the deviation of the X P\ mass from the center of gravity of the Xcj 
states would indicate a departure from first order perturbation theory, since the 
spin-spin potential is a contact potential, which survives only with the finite wave 
function at the origin. Thus, this potential gives rise to hyperfine splitting between 
the triplet (S = SI + S2 = 1) and singlet (S = SI + S2 = 0) states only for L 
= S-wave states, and not for P-wave or any other higher L-states. This is the 
direct consequence of the long-range confinement potential having been assumed to 
be pure Lorentz scalar. There is no fundamental justification for this assumption, 
although support for it is observed in the results of quenched lattice calculations, 
for example, those of Bali, Schilling and Wachter [4T] . 

2.1.3 Calculational Techniques 

In 1974, Wilson (1974) |12j showed how to quantize a gauge field theory on a discrete 
lattice in Euclidean space-time preserving exact gauge invariance, and applied this 
calculational tecnique to the strong coupling regime of QCD. In these Lattice Gauge 
Calculations, space-time is replaced by a four dimensional hypercubic lattice of size 
L 3 T. The sites are separated by the lattice spacing a. In more recent calculations, 
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asymmetric lattices have been used in which the lattice spacing for time is chosen 
to be much smaller than for space, with a s as small as 0.07 fm, and a s /a t as large 
as 3. The overall size of the lattice is generally 1-4 fermis. 

The quark and gluon fields are defined at discreet points on the lattice, and 
physical problems are solved numerically by Monte Carlo simulations using powerful 
computers, requiring only the quark masses as calculational input. Early lattice cal- 
culations were done in the " quenched" approximation, in which no quark- ant iquark 
pairs are allowed to be excited from the QCD vacuum (or sea). Recently, the vast 
improvement in computers and computational techniques have made it possible to 
do unquenched calculations in which uu, dd, and even ss quark-antiquark pairs may 
be excited from the vacuum. 

A useful variant of lattice calculations for heavy quarks is the so called Non 
Relativistic QCD (NRQCD), pioneered by Lepage and colleagues JI3]. As heavy 
quarks (i.e. charm and bottom) are generally non-relativistic, renormalization group 
techniques may be used in these calculations to replace the relativistic Dirac action 
for the heavy quarks on the lattice by a non-relativistic Schrodinger action, which 
simplifies the calculations considerably. However, even with such simplifications, 
Lattice Gauge Calculations generally require immense computing resources. Yet 
despite the large amount of computation required lattice calculations are able to 
provide useful insight into static potentials for QCD, masses of charmonium states, 
and many decay characterstics. For an excellent review of lattice methods and their 
application to charmonium physics, see the Physical Report article by Bali 03]. 

Another important technique used for QCD calculations is the QCD sum rule 
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technique, which was introduced in 1979 by Shifman, Vainshtein and Zakharov [45J. 
The basic premise of the QCD sum rule technique is that the QCD vacuum is 
populated by large fluctuating fields whose strength is characterized by gluon and 
quark condensates 

< Q\(a s /n)G$G^ v \Q >, and < 0|gg|0 > 

When a pair of quarks is injected into this active vacuum, its dynamics is deter- 
mined by the characterstics of the vacuum, and the subsequent formation of hadrons 
can be reliably calculated by dispersion relations. The heavy quark charmonium and 
bottomonium states are uneffected by light quark condensates, and are sensitive only 
to the gluon condensate. The sum rule technique is nearly saturated even by the 
lowest excitations of each J PC , and is therefore next to impossible to apply to radial 
excitations or to resonances with higher orbital excitations than J = 2. 

One early spectacular success of the sum rule calculations was the correct pre- 
diction of the mass of the charmonium ground state, ^(l 1 ^) [HI]. Highly successful 
QCD sum rule calculations of IP charmonium states were made subsequently made 
by Reinders, Rubinstein and Yazaki |47j . 

Yet another QCD calculation technique is used to attempt to improve the rel- 
ativistic problems caused by the singular nature of the Coulombic potential. The 
"smearing" technique is used to spread the singularity into a small region around 
the origin. This technique is used by Godfrey and Isgur [38] . 

In addition to lattice and sum-rule predictions, prediction may also be obtained 
using perturbative QCD (pQCD), in analogy to the perturbative QED (pQED) used 
for positronium. The technique is valid for large momenta and small values of a s , 
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and has been used for charmonium annihilation calculations under the assumption 
that the qq wave function is purely color singlet, and that the annihilation is a short 
distance process. Strong radiative corrections for the electromagnetic, radiative, 
and hadronic decays of S and P wave quarkonia have also been made under this 
assumption by several authors jlH] jSO] [EE] 

It has been argued, notably by Bodwin et al. that the hadronic decays 
of charmonium, and particularly the P wave states, require taking acount of the 
possible qqg components, with qq in a color octet. The octet components may be 
determined only empirically from the existing data, but it is argued that the 3 Pi 
and 1 P\ hadronic decays take place dominantly through their octet components, 
and therefore these provide a good estimate. Using these suggestions, radiative 
corrections for octet decays have been calculated by several authors, and these have 
been summarized by Vairo (SHI- A couple of problems should be noted about pQCD 
predictions. One problem is that unlike the case of positronium, in charmonium 
m(qq) 7^ 2m(q), and thus there is an ambiguity about whether to evaluate pQCD 
expressions using m(qq)/2 or m(q), which may affect calculations of hadronic decay 
widths. The second problem relates to strong radiative corrections since, unlike 
the pQED case, a s is large enough (~ 0.3) for charm quarks that the lowest order 
gluon radiative corrections are often very large, up to 100%. Despite this, radiative 
corrections calculated with pQCD have in some cases led to favorable agreement 
with experimental results. 



2.2 Theoretical Predictions 
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A large number of theoretical predictions exist for the masses of charmonium states, 
including the /i c ( 1 Pi). The majority of these predictions are based on potential 
model calculations, and differ only in the choice of the common parameters; the 
strong coupling constant, a s , and the charm quark mass, m c . Some of these potential 
model and parameter choices are listed in Table 2.2. The Cornell potential, and the 
QCD based potentials modeled after that by Richardson are the most commonly 
used. Minor variations on these are used by several authors. 

The linear confinement potential is generally assumed to be scalar (S). Some 
authors have considered vector confinement (V) and mixtures of the two (V+S). 
Most calculations are non-relativistic. Some include relativistic corrections at the 
level of v 2 1 c 2 . The only predictions based on lattice calculations are those by Bali |H] 
and Okamoto [5*4*] . 

A summary of predictions for the properties of the 1 P\ state are given in Tables 
2.2 - 2.4. Mass and width predictions are shown in Tables 2.2 and 2.3. The x Pi 
mass is generally predicted to lie within a few MeV of the centroid of the Xcj states; 
< M(xj) >expt = 3525.3 MeV. The predictions for the total width of 1 P\ lie in 
the range of 500-1000 keV. Predictions for the partial widths of various l P\ decays 
are shown in Table 2.4. The most prominent decay of the 1 P\ is expected to be 
the radiative transition to r] c , with predicted partial widths in the range of several 
hundred keV, based on the measured width of the El transition 3 Pi — > jJ/i/j. 
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Author 


Year 


AMht 
(MeV) 


Potential 




m c 
(GeV) 


Conf. 


Eichten gU] 


1979 





Cornell 


0.341 


1.84 


S 


Ono EE! 


19* 


i2 


+1.0 










Gupta [SH 


19* 


i2 


-1.4 


Cornell 


0.392 


1.2 


s* 


McClary [SB] 


19* 


13 


+5.3 


Cornell 


0.341 


1.84 


s* 


Moxhay 


19* 


13 


~ 


Richardson 




1.5 


Tensor* 


Godfrey gH] 


19* 




+5 


Cornell* 


0.34 


1.5 


S 


Pantaleone jHEl 


19* 


i6 


-3.6 


QCD 


0.24 


1.48 


V + s 


Olsson [SH 


19* 




+5.4 + 0.8 


Cornell 






V + s 


igi |S2! P! 


19* 




+24.1 + 2.5 


QCDt 




1.506 




Pantaleone jHSj 


19* 


^8 


-1.4 


QCD 


0.33 


1.48 


V + s 


Gupta [HZ] 


19* 


^9 


-2.0 


Cornell 


0.392 


1.2 


V + s* 


Badalian |fi£] 


1990 


-3.1 + 4.9 


QCD 


< 0.39 


1.35 




Dixit jMJ 


1990 


-8.0 + 8.0 


Power Law 






V 


Galkin [70] 


1990 


+8 


Cornell 


0.52 


1.55 


v* 


Chakrabarty [7T] 


1991 


+ 70 


Power Law 


0.25 


1.5 


V 



Table 2.2: A summary of potential model predictions for the mass and width of 
the l Pi (part 1). The hyperfine mass splitting AM hf =< M(xj) > -M( l P x ) 
where < M(xj) > e xpt = 3525.3 MeV. Asterisks denote inclusion of some relativistic 
corrections. Daggers denote modified potentials. 
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Author 


Year 


AM hf 


Potential 




m c 


Conf. 






(MeV) 






(GeV) 




Fulcher [?2] 


1991 


-3.0 


QCD 


1.30 


0.54 


V ± S 


StubbejZSj 


1991 


-23 ± 13 


QCD 




1.5-1.8 




Lichtenberg |73| 


1992 


+4 


QCD 




1.82 


V 


Lichtenberg |74j 


1992 


-1.3 ± 4.2 


QCD 


< 0.44 


1.65 


s 


Beyer |ZZj 


1992 


+15 


Cornell 


~ 0.5 


1.9-2.3 


s 


Halzen [ZH] 


1992 


-0.7 ±0.2 


Cornell 


0.28 


1.2 




Chen jZH] 


1992 


-1.4 ±0.5 


QCD 


0.22 


1.48 




Eichten [HI] 


1994 


±1 


QCD 


0.31 


1.48 




Gupta jH2 


1994 


-0.9 


Cornell 


0.392 


1.2 


V ± s 


Zeng |H33 


1995 


±6 


Cornell^ 




1.53 


s 


Chen jHU 


1996 


-5± 1 


QCD 


runs 


1.478 




Bali glj 


1997 


-1.5 ± 2.5 


Lattice 


0.183 


1.33 


s* 


Okamoto [HI 


2002 


-1.5 ± 2.6 


Lattice 








Lahde [HS] 


2002 


±12 


QCD 


0.38 


1.50 


V 


Ebert jHHj 


2003 


-0.7 


Cornell 


0.314 


1.55 


V ± s* 


Recksiegel jHTj 


2003 


-0.8 ±0.8 


QCD 


runs 


1.243 





Table 2.3: A summary of potential model predictions for the mass and width of 
the 1 Pi (part 2). The hyperfine mass splitting AM hf =< M(xj) > -M^Px) 
where < M(xj) >expt = 3525.3 MeV. Asterisks denote inclusion of some relativistic 
corrections. Daggers denote modified potentials. 
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Authors 


r(^ c 7) 


r(j/^7r°) 


r(J/^7T 7T ) 


T(hadrons) 


r (total) 




(keV) 


(keV) 


(keV) 


(keV) 


(keV) 


Renard [22] 


240 






370 


500-1000 


Novikov [EE] 


975 






60-350 




McClary [21] 


485 










Kuang [Mj 




2 


4-8 


54 


395-400 


Galkin [7U] 


560 










Chemtob [H2] 




0.006 


53 






Bodwin [221 


450 






530 


980 


Chen [72] 




0.3-1.2 


4-14 


19-51 


360-390 


Chao [Hi 


385 










Casalbuoni [§U] 


450 










Ko [m 


400 


>1.6 








Gupta [H2] 


341.8 











Table 2.4: A summary of predictions for the partial widths of various 1 P\ decays. 



Chapter 3 

Experimental Apparatus 



The experimental set-up for the Fermilab E835 experiment, which follows the origi- 
nal E760 experiment, consists of a hydrogen gas-jet acting as a fixed proton target in 
the path of a circulating beam of antiprotons in the Fermilab Antiproton Accumu- 
lator. Decay products from the resulting pp annihilations are detected in a detector 
system with cylidrical geometry surrounding the interaction region. It has no mag- 
netic field, and is designed to optimally detect and identify electrons, positrons, and 
gammas. As such, it is particularly suitable for the study of those states whose 
decay involves the vector (1 ) states of charmonium, J ftp and ip' which have sig- 
nificant branching fractions into e + e~ final state, and/or hadrons which can decay 
into multiple gammas, in particular 7T° and rj, which have large branching fractions 
for decay into two photons. Thus, E835 is particularly suited to the search for h c , 
the 1 Pi resonance of charmonium, one of whose principal decay modes is expected 
to be pp — s- 1 Pi — > J/ip + n° — > e + e~ + 77. The study of this decay is the primary 
objective of this dissertation. Several other possible decay channels of h c were also 
studied, and are also described. In this chapter we describe the experimental set-up 
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in some detail. A fuller discussion of the E835 experimental setup can be found 
in jnS]- E835 had two runs, one in 1997 with 141.4 pb _1 of luminosity, and one in 
2000 (also called E835p) with 113.2 pb _1 of luminosity. This dissertation is devoted 
to the search for h c in the year 2000 data, although we also refer to the E760 and 
E835 (1997) data. 

The excitation of a charmonium resonance may be described in terms of the 
Breit Wigner formula for resonance cross sections: 



where Mr is the mass of the resonance, Tr is the total width of the resonance, 
B in = B(pp —>■ R) and B out = B(R — > finalstate) are the branching ratios to the 
initial and final states, J is the spin of the resonance, m p is the proton/antiproton 
mass, and y/s is the energy in the center of mass frame. 

Charmonium resonances at Fermilab E835 are studied by sweeping the pp center 
of mass energy through the resonance region, and measuring the yield of various 
decay channels as a function of center of mass energy. An example of such a scan is 
shown in Figure 3.1. 

The antiproton beam is stochastically cooled so as to be as monochromatic as 
possible, but it still maintains a small but finite momentum spread, which must 
be convoluted with the Breit Wigner resonance cross section and multiplied by the 
efficiency and acceptance of the detector to obtain the measured cross section: 



vr(2J+ 1) 
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^5 1 



Solid: Fit to data 
Dashed: E „ distribution 

cm 




3555 3560 3565 

E rm (MeV) 



Figure 3.1: A scan of the Xc2 resonance of charmonium using the E835 detector. The 
half width of the Xc2 resonance is ~ 1.8 MeV. The center of mass energy distibution 
of the beam, with typical half width of ~ 500 keV, is shown by the dotted curve. 



where G(ys' — yfs) is the function describing the center of mass energy distri- 
bution of the antiproton beam as a function of the energy vF, when ^Js is the 
center of mass energy at which the cross section is being measured, and e and a 
are the detector efficiency and acceptance respectively. In general, there is also a 
background cross section <7bkg with a relatively slow variation with center of mass 
energy, and the Breit Wigner resonance sits on top of this background. 

The number of events iV, observed at a given energy point J~sl is given by: 



Ni = / d(t) x [a(y/sl) + a bkgd (yfsl)}dt 



(3.3) 



where Ci(t) is the instantaneous luminosity at energy point i, <Jbk g d is the back- 
ground cross section at that energy, and (j(y/s~i) is given by equation 3.1. The 
background cross section is generally measured with data taken several half-widths 
away from the resonant energy. Resonance parameters such as mass, width, and 
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branching ratios may be calculated from measurements of events JVj for a particular 
set of decay products at various energy points ^/sl across the resonance region. 

The plot of cross section versus energy is known as the excitation curve, and the 
area under this curve is given by: 

f°° IT 

A = j a(y/s)dy/s = -a peak T R (3.4) 

where, 

4?r(2 J + l)B in B out ae 
apeak ~ Ml- Ami 

For relatively broad resonances which have a width much larger than the spread 

of the beam, Y R may be measured directly from the shape of the excitation curve. 

Even when the width of the beam affects the shape of the excitation curve, causing 

it to differ from a pure Breit-Wigner for narrow resonances, the area under the curve 

is conserved. 

3.1 The Fermilab Antiproton Accumulator 

Antiprotons for experiment E835 are produced in pp collisions using protons from 
the Fermilab Main Injector. Protons are initially injected into the accelerator system 
as H ions and boosted to 750 keV by a Cockroft- Walton accelerator. They are 
then injected into the Fermilab Linac, where they are accelerated to 400 MeV. 
At this point they are made to pass through a carbon foil and stripped of their 
electrons, transforming them from H~ ions to bare protons. These protons are then 
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accelerated in two rings. The first, known as the Booster ring, brings them to an 
energy of 8 GeV, and the second, the Fermilab Main Injector, which brings them to 
an energy of 120 GeV. 

These 120 GeV protons from the Main Injector are then focused onto a nickel 
target, where they produce a wide variety of secondary particles. For every mil- 
lion protons hitting the target, approximately 20 antiprotons are produced, with 
energies which have a distribution with a peak around 8.9 GeV. These, along with 
other negatively charged particles, are focused by a lithium lens towards a mag- 
net, which then directs them to the Debuncher Ring. Negatively charged pions and 
muons which accompany the antiprotons decay on the way to the debuncher, while 
electrons miss the beam aperture as they lose energy due to synchrotron radiation. 
This results into an essentially pure antiproton beam reaching the Debuncher, the 
primary purpose of which is to reduce the initial large divergence of the captured 
antiprotons, and use the stochastic beam-cooling systems to create a brighter beam. 
This beam is injected into the Antiproton Accumulator Ring which is located be- 
side the Debuncher Ring. The various elements of the Fermilab accelerator system 
relevant to the production of antiprotons are shown in Fig. 3.2. 

The Antiproton Accumulator is designed to accumulate antiprotons, a process 
which is called stacking, for use in the Tevatron . Antiprotons are stacked in the 
accumulator at a rate of ~ 2 — 4 x 10 10 p/hour, until a beam of approximately 6 x 
10 11 p has been collected. The Antiproton accumulator also reduces the momentum 
spread of the beam by a process known as stochastic cooling (HHl- Once enough 
beam has been accumulated and cooled, it can either by extracted for the Tevatron, 
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Fermilab Antiproton Source 




Figure 3.2: Schematic of the Fermilab Antiproton Source and its position relative 
to the FNAL Main Ring. The location of the E835 experiment in the Antiproton 
Accumulator is also shown. 

or, when the E835 experiment is running, decelerated from the momentum of 8.9 
GeV to the momenta needed by E835, which ranged from ~ 5 — 7 GeV/c for the 
E835 year 2000 run. This deceleration is done with an RF cavity operating at the 
second harmonic of the beam revolution frequency. This cavity has a maximum 
RF voltage of 3 kV, allowing for a deceleration rate of 20 MeV/s. The deceleration 
is performed in several steps constituting a 'ramp'. The momenta to which the 
p beam must be brought for the various charmoinum states are shown in Table 
3.1. The beam parameters are measured and corrected at the end of each ramp, 
and the magnet settings of the dipoles and focusing quadrupoles of the Antiproton 
Accumulator are appropriately adjusted. 

There are 48 horizontal and 42 vertical Beam Position Monitors (BPMs) posi- 
tioned around the Accumulator Ring. 

There are 38 dipole magnets that bend the beam horizontally around the ring. 
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cc State 


E cm (GeV/c 2 ) 


Pbeam (GeV/c) 


Vc 


2.9797 


3.6919 


J/4 


3.0969 


4.0657 


Xo 


3.4151 


5.1931 


Xi 


3.5105 


5.5502 


'Pi 


3.5262 


5.6009 


X2 


3.5562 


5.7246 


1p' 


3.6860 


6.2321 



Table 3.1: Masses of charmonium states and the necessary p beam momenta. Masses 
are from PDG2002 0. 

There are also 48 horizontal and 42 vertical Beam Position Monitors (BPMs) po- 
sitioned around the Accumulator Ring. They are positioned so that the orbit dis- 
placement created by each individual magnet is measured by at least one BPM. 

3.1.1 Stochastic Cooling 

Once the desired beam momentum is reached the beam is again stochastically cooled. 
There are two types of cooling, transverse and longitudinal. The transverse cooling 
reduces the growth of the beam which occurs due to multiple scattering of the an- 
tiprotons with the target and with residual gas in the Accumulator. The transverse 
cooling reduces the size of the beam so that at the gas jet target 95% of the beam is 
contained in an approximately circular region of radius 2.45 mm. The longitudinal 
cooling reduces the momentum spread in the beam and achieve the small beam en- 
ergy spread shown in Figure 3.1. A distribution of the center of mass energy spread 
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of the beam for all runs in the year 2000 E835 1 P\ data stacks is shown in Figure 
3.3. It ranges from a = 0.2 — 0.55 MeV, with an average of 0.32 MeV. Both the 
transverse and longitudinal cooling systems are explained in the following. 

20 
•8 

6 
'4 

2 

•o 

8 
6 
4 



"0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 
a*u* (MeV) 

Figure 3.3: Distribution of the center of mass energy spread (r.m.s.) for all runs of 
the year 2000 E835 1 Pi data stacks. 



Transverse Stochastic Cooling 

As a p passes one of the several pickup electrodes positioned around the ring, its 
deviation x from the central orbit position is detected. A correction can then be 
applied by transmitting a signal to a kicker electrode, which is located an odd 
number of quarter-wavelengths of betatron oscillation downstream from the pickup. 
The kick is timed so that it is delivered when the particle detected by the pickup 
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arrives at the kicker. It causes the antiproton to have a transverse position x — gx 
at the pickup, where g is the system gain. For a "beam" which consists of a single 
particle, a single kick would be enough to correct its orbit. However, since we are 
dealing with a beam of N antiprotons, each of which affects the motion of the others, 
the effect of each kick which is delivered by the kicker is smaller. Furthermore, the 
presence of the other particles means that the pickup detects the mean deviation of 
a portion of the beam, and delivers an appropriate kick. The effect of this, along 
with the fact that the system gain g cannot be exactly 1, is that cooling the beam 
requires not one, but many, kicks. The cooling principle, though, is applicable for a 
beam of any size. 

Longitudinal Stochastic Cooling 

Transverse cooling, as described above, decreases the physical size of the beam by 
decreasing the amplitude of the betatron oscillations, but this has only a marginal 
effect on the momentum distribution of the beam. For that purpose, longitudinal, 
or momentum, cooling must be done. 

In momentum cooling, it is necessary to detect variations, Ap, from the mean, 
or central, beam momentum < p >. The mechanism for momentum cooling is 
similar in nature to that used for transverse cooling. In this case, however, a band- 
pass or "notch" filter is used in the pickup-kicker network, so that particles nearest 
the central momentum, which corresponds to the central frequency, are the least 
affected. That is, the presence of the filter allows a positive correction to the slightly 
low frequency particles, and a negative kick to the slightly high frequency ones, while 
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leaving the particles near the central orbit frequency alone. 

Transverse cooling in the Debuncher reduces the emittance of the p beam from 
~ 207r to ~ 7tt mm-rad. Longitudinal cooling reduces Ap/p from the ~ 0.2% 
achieved by the RF to ~ 0.09% [2E]- At this point, the p beam is transferred into 
the Accumulator Ring, and the Debuncher is ready to accept a new p batch. A 
similar cooling processes is then performed in the Accumulater Ring. 

The vacuum in the Accumulator Ring is typically of the order 3 x 10~ 10 torr, 
giving a fully stacked beam a lifetime (£1/2) of approximately 1000 hours with the 
jet target off. With the jet target running, the beam lifetime is reduced substatially, 
to ~ 50 hours. The E835 target and detector is located in a low dispersion region 
of the Accumulator, where a momentum spread of Ap/p of 10~ 4 corresponds to a 
longitudinal displacement of only 50 ^m. 

3.1.2 Measurement of Beam Center of Mass Energy 

The Lorentz-invariant center of mass energy of the beam/target interaction is de- 
termined from the antiproton energy and momenta in the lab frame by the relation: 



E cm = (Ep(lab) + m p C 2 f 



2 2 
- C Pl ab 



(3.6) 



Since, 




(3.7) 
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we obtain 



E 2 cm = 2m p c 2 (m p c 2 + E p{lab) ). 



(3.8) 



EpQ a b) can also be written as: 



E p{ i ab ) = 'jrripC 2 = m p c 2 / Jl- (v^/c 2 ) 



(3.9) 



Thus, the measurement of the center of mass energy depends only on the mass 
of the proton/antiproton and the velocity of the antiproton in the lab. Since the 
proton/antiproton mass is known to be 938.271998 MeV±38 eV Q, the precision of 
the energy measurement rests depends only on the precision of the determination 
of the antiproton velocity. This velocity is obtained by multiplying the frequency of 
circulation of the antiprotons in the accumulator with their orbit length: 



The uncertainty in the E p i ab measurement can then be calculated by dfferenti- 
ating this expression with respect to f p and L orbit to give: 




(3.10) 



Thus 




(3.11) 
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S J2Z*L = Plj}( S Jl + S J?**) (3.12) 

Ep(lab) Jp J-'orbit 

where ftp and 7p are the relativistic factors for the antiproton in the lab frame, 
and 5 Epilab) is the uncertainty on Ep(i a b). 

The corresponding expression for SE cm can be similarly be written as: 

5E cm _ fiplp ^ 5fp SLprbu , ^ ^ 

E C m 2(1 +7p) fp L or i,it 

As an example of the effect of the uncertainties on the L^ut and fp measurements 
on the E cm uncertainty, at the ip' center of mass energy, dE cm /df = 113.2 keV/Hz 
and dE cm /dL = 149.3 keV/mm [HZj. The ip' center of mass energy is of crucial 
importance for this calculation, since this is the energy at which the orbit length is 
measured. Typical values for fp and L or ut are 0.6 MHz and 474 m, while typical 
uncertainties are in the range of 5f ~ 0.1 Hz and 5L « 1.2 mm, which leads to: 

^ = 0.6 x 10- 7 , ^ = 2.5 x 10~ 6 . (3.14) 
fp L 

Thus, the uncertainty in the measurement of the center of mass energy is dom- 
inated by the uncertainty determination of the orbit length. Measurement of /, L, 
5f, and 5L is described in the following subsections. Using these measurements we 
may constrain the uncertainty in the center of mass energy to within 180 keV for 
■s/s in the ip' region, and 60 keV for \fs in the J/ip region j^Tj. These uncertainties 
are further reduced by scanning the ip' resonance and calibrating the orbit length 
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using the accepted value for the if)' center of mass energy, as will be described in 
Section 3.1.4. 

3.1.3 Measurement of fp 

The orbit frequency of the beam is determined using the Schottky noise spectrum, 
which is created by the sum of the pulses generated by the passage of each particle 
in the vicinity of the pickup. The Schottky noise power spectrum is measured by 
using a Schottky pickup and a network analyzer. This power spectrum P(f), is 
proportional to the frequency spectrum of the orbiting antiprotons (dN/df), or one 
of the harmonics of that frequency spectrum, by the relation: 

P(f)=2n(ef) 2 dN/df (3.15) 
where e is the electric charge. 

An example of the Schottky noise power spectrum is shown in Figure 3.4 for the 
127th harmonic. This power spectrum is fit in order to determine the peak value to 
within 10 Hz, corresponding to an uncertainty in the fundamental frequency of less 
than 0.1 Hz. Dividing the central frequency and the width of the peak by 127 for this 
plot gives an orbit frequency of (625.366 ± 0.004) kHz, and thus Sf/f ~ 10~ 7 [93J. 
Note that the y-axis in Fig. 2 is plotted in a logarithmic scale. 

The exact relation between the spectrum of the beam momenta and the fre- 
quency spectrum (measured via the Schottky noise power spectrum) is given by the 
relation P?j : 
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Figure 3.4: Measured Schottky noise power spectrum of the 127* harmonic of the 
antiproton beam with momentum 6232 MeV/c (the resonance). The central fre- 
quency is 79.421474 MHz 93 . An uncertainty Sf of ~ 10 Hz in the 127 th harmonic 
corresponds to a Sf of ~ 0.1 Hz in the fundamental. 



^ = - 1 x d l . (3.16) 
Pp V f 

where 7] is a parameter known as the slip factor, and is given by: 



,= 1-1, (3.17) 

where 7t, the gamma factor at the transition energy, is a parameter determined 
by the machine lattice [98J [99J. Since df/f can be well measured directly from the 
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Schottky spectrum, to obtain the beam momenta we must determine 77 as a function 
of beam energy. This can be done with two different techniques [93j . 

The first technique is to determine 77 from the synchrotron frequency as a function 
of peak RF voltage. With RF power on, the energy of the orbiting antiprotons will 
oscillate around the central energy with a characteristic synchrotron frequency f s 
given by: 



where V r f and Jrf are the peak RF voltage and RF frequency, E is the beam 
energy, h is the harmonic number, and <p s is the synchronous phase, which is 
above transition (7 > 7 t ) and 7r below transition. The synchrotron frequency may 
be determined to a precision of 1%, and the uncertainty in the measurement of 77 is 
dominated by the uncertainty in the RF voltage, which is of the order of 5%. 

The second technique is to determine r\ from j t (eq. 3.17). j t m ay be deter- 
mined by varying the dipole magnetic fields in the absence of RF and measuring 
the resulting change in /. The gamma factor at the transition energy may then be 
determined by using the relation: 




(3.18) 



E 2nh p 2 



dB/B = ^df/f 



(3.19) 



with the uncertainty in 77 related to the uncertainty in 7 t by: 



5 V = {2hl)8 lt 



(3.20) 
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As this method requires knowledge of the magnetic fields in the dipoles, and 
these have large systematic errors, it is used primarily to cross check the result from 
the first technique. 

3.1.4 Measurement of L or ut 

The Measurement of the orbit length; L OT ut is done by measuring the difference in 
length between the current orbit and a reference orbit using a system of 48 Beam 
Position Monitors. The reference orbit is chosen to be that which has a center 
of mass energy at the peak of an easily detectable resonance whose mass is well 
known. E835 uses the ij}' resonance, which has a mass of (3686.111 ± 0.025 ± 0.009) 
MeV |100j . and which is observed through its inclusive decays into J/ip followed by 
the subsequent decay J/xjj — > e + e~, as well as its direct decay ip' — > e + e _ , which 
give clear signal e + e~ signals in the detector. 

The reference length L may then be determined from the mass of the ifi' resonance 
and the measured beam frequency at the ip' energy using the relation: 



which may be derived from Equations 3.12 and 3.13 for the case where the 
center of mass energy is equal to the ip' mass. The error in this measurement is 
then dominated by the error in the knowledge of the mass of the ip', and generates 
an uncertainty in L re f which is given by: 




(3.21) 
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5L ref = L ref x —= — x 8Mm = 0.6 mmt2me,s0.28 = 0.17mm. (3.22) 

Once the reference orbit length is known, the difference between the length of the 
orbit at any given energy and the length of the reference orbit can be determined 
using the beam position monitors (BPMs). The data from these BPMs must be 
used in a piecewise manner to determine the change in the orbit length, or by using 
a constrained fit using a detailed model of the Antiproton Accumulator lattice. The 
overall uncertainty in the orbit length is of the order ± 1 mm, which is the largest 
contributor to the uncertainty of the center of mass energy in E835. 



3.2 The Gas Jet Target 

The E835 target is a hydrogen gas jet which ejects clusters of hydrogen perpendicu- 
larly to the antiproton beam axis at a rate of ~ 1000 m/s |101j . With the hydrogen 
jet target on, a typical 50 mA beam lasts for 2 or 3 days before being depleted by 
a factor ~ 3; the beam lifetime is ~10 times longer with the gas jet turned off. 

In order to maintain a constant instantaneous luminosity during data taking, 
the density of the jet may be varied during operation to compensate for the loss 
of antiprotons in the beam. In a typical run the jet density is varied to keep the 
instantaneous luminosity in the range of 2 — 3 x 10 31 s _1 cm -2 . This corresponds to a 
minimum bias trigger rate of 3 MHz, which is close to the maximum sustainable by 
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the detector and data acquisition systems, and thus optimizes the run time alloted to 
E835. The constant instantaneous luminosity also allowed for better understanding 
of effects such as event contamination due to pileup, which is strongly luminosity 
dependent. A plot showing the variation of the gas jet density antiproton current, 
and instantaneous luminosity is shown in Fig. 3.5. 

The gas jet used in E835 is of a type known as a cluster jet because the core 
of the jet is made up of small droplets, or clusters, of condensed hydrogen. This 
cluster jet is created by allowing hydrogen, which is kept at high pressure and low 
temperature, to expand through a convergent-divergent nozzle. The expansion of 
the gas through the nozzle is isentropic, and the the sudden decrease in pressure 
and temperature caused by the expansion leaves the gas in a supersaturated state, 
favoring the formation and growth of a jet of clusters whose size varies from 10 7 to 
10 8 molecules. 

A schematic of the E835 gas jet nozzle is shown in Fig. 3.6. The nozzle is 
trumpet shaped, and has an opening angle of 3.5°, a divergent length of 8 mm, 
and a throat diameter of 37 /xm. As can be seen from the shape of the isentropes 
on the P-T diagram of hydrogen (shown in Fig. 3.7) the density of the jet may 
be maximized by allowing the expansion to begin at the highest possible pressure 
and the lowest possible temperature, which puts the hydrogen as close to the H 2 
saturation curve as possible. In order to keep the temperature low, a helium cryo- 
cooler was installed at the final stage of the hydrogen line, allowing operation with 
hydrogen gas temperatures as low as 20° K. Further manipulation of the point on the 
P-T curve where the expansion starts may also be done by changing the pressure; 
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Figure 3.5: Jet target density and antiproton beam current during ~ 30 hours of data 
taking at E cm = 3526 MeV. As the beam current steadily decreases, the jet density 
is increased to maintain an approximately constant instantaneous luminosity. The 
jet density is shown by the solid blue points, the antiproton current by the solid red 
curve, and the instantaneous luminosity by the open black circles |lU9j . 
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reduced pressure at the nozzle allows for correspondingly lower densities. The open 
circles in Fig. 3.7 show the operating points used by the E835 gas-jet target; these 
lie directly above the H 2 saturation curve, and give the range of jet densities shown 
in the upper curve. 




Figure 3.6: Detail of the gas jet nozzle and the formation of hydrogen clusters. 



Using the gas jet as a target inside the Antiproton Accumulator required setting 
up a system of turbo-molecular pumps in order to maintain the high vacuum in the 
beam pipe which is required to preserve a high quality p beam. The ability of this 
pumping system to maintain high vacuum is further enhanced by separating the 
cluster jet stream from the remaining gas exiting the nozzle. A differential pumping 
scheme is therefore used (as shown in Fig. 3.8), in which the jet crosses a series of 
chambers which are independently evacuated. These chambers have installed into 
them a series of ten turbo-molecular pumps (TMPs), eight of which have a capacity 
of 1000 liters per second, and two of which have a capacity of 3500 liters per second. 
As TMPs have a low compression ratio for hydrogen, two additional TMPs, three 
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Figure 3.7: The P — T phase diagram for H2. The operating points utilized by 
E835 are shown by open circles and the jet density corresponding to those operating 
points are shown by filled circles |§3*] . 
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Figure 3.8: Layout of the hydrogen jet target and pumping system. 
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positive displacement blowers and two roughing pumps were arranged in a cascade 
configuration in over to avoid limiting the pressure in the high vacuum zone of each 
pump due to the rough vacuum (see bottom of Fig. 3.8). The chambers which 
are upstream from the interaction zone are labelled Jl, J2, and J3, and the TMPs 
in these chambers remove the part of the gas which does not clusterize into the 
core of the jet. The chambers which are downstream from the interaction zone are 
labelled Rl, R2 and R3, and the TMPs in these chambers are used to remove the 
core jet once it has crossed the interaction zone. These must be removed as only 
a small percentage of the protons in the clusters interact with the antiprotons in 
the beam. By using this system, we can reduce the number of interactions of the p 
beam outside of the interaction zone to 5% of that inside the interaction zone. 

In order to maximize the cluster density during data taking, the control system 
which regulated the temperature of the nozzle and the pressure of the hydrogen line 
was automated. This system allowed for regulation of the pressure to within 0.5 psi 
and the regulation of the temperature to within 0.05° K, with a response time of 10 
s. Thus the densities of the hydrogen could be varied from ~ 1 x 10 14 atoms/cm -3 to 
~ 3 x 10 14 atoms/cm -3 over the lifetime of an antiproton stack. The jet diameter in 
the interaction region is 6 mm, which is only slightly larger than the beam diameter 
at this point (~ 5 mm). This diameter is regulated by the geometry of the skimmer 
between the second and third vacuum chambers. The system of collimators which 
cross the axis of the gas jet confine the direction of the jet to be perpendicular to 
the beam axis to within 2°. The gas jet beam pipe is stainless steel with a thickness 
of 0.18 mm in the region where the secondary particles pass into the main detector, 
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which will be described in the next section. 

3.3 The E835 Detector 

The E835 Detector is designed to detect electromagnetic final states, 7, e + and 
e~ , with two or more of these particles resulting from the decay of charmonium 
resonances. The detector is also designed to perform at a high interaction rate 
because pp total cross sections are very large, being ~ 70 mb in the y/s = 3-4 GeV 
region. A complete and detailed description of the E835 detector may be found in 
the recently published in Nucl. Inst. Meth. A j^j; much of the following discussion 
is based on that paper. 

The schematic of the E835 detector system is shown in Fig. 3.9. It has a cylin- 
drical geometry about the antiproton beam axis, with full coverage of the azimuthal 
angle for all its central components. The system consists of the central and for- 
ward electromagnetic shower calorimeters to measure the energy and momentum 
of electrons, positrons and photons, a system of inner detectors to detect charged 
particles, a Cerenkov counter to discriminate electrons and positrons from heavier 
charged particles, and a luminosity monitor to measure the interaction luminosity. 
The inner detectors are made up of three plastic scintillator hodoscopes, four layers 
of drift tubes (straws), scintillation counters for forward angle veto, and two scin- 
tillating fiber detectors. The inner detectors are contained in a cylinder of radius 
17 cm and length 60 cm; and their total thickness is less than 7% of a radiation 
length for particles crossing at normal incidence. Other than the luminosity moni- 
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tor, these detectors are all highly segmented to allow for a higher rate of signal and 
are equipped with time-to-digital converters (TDCs) to reject out of time signals. 

The detector's polar acceptance for charged particles is 6 = 15° — 65°, and 
10.6° — 70° for photons in the central calorimeter (CCAL). A forward calorimeter 
(FCAL) extends the photon acceptance down to 2°, but it is, as in the present 
measurements, used primarily to provide a veto. 



E835 EQUIPMENT LAYOUT (Y2K) 




Figure 3.9: Side view of the E835 detector system. 
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3.3.1 The Hodoscopes and the Veto Counters 

There are three scintillator hodoscopes in the inner detector; these are labeled HI, 
H2, and H2', as well as a set of forward veto counters. These are all segmented, and 
each unit is made of a plastic scintillator connected by a light guide to a phototube. 
The Bicron 408 plastic which is used for the scintillators has an index of refraction 
of 1.58 and a density of 1.03 g/cm~ 3 . 

The three hodoscopes HI, H2, and H2' are all designed to be axially symmet- 
ric around the antiproton beam. They are segmented azimuthally into planes of 
constant angle 0. The innermost hodoscope, HI, consists of 8 plastic scintillators 
which form a cone shaped structure around the segment of the beam pipe which 
is attached to the jet target body. The thickness of the HI scintillators is 2 mm; 
they provide full coverage in the azimuthal angle 0, and coverage from 9° to 65° 
in the polar angle 9. The HI light yeild is about 10-20 photoelectrons for a single 
minimum ionizing particle. HI is located at distance of 2.2 cm from the beam axis. 
It is the innermost element of the entire E835 detector. 

The outermost hodoscope, H2, is a cylindrical device of radius 17 cm, which 
is made of 32 segments. These segments are 60 cm long, 3 cm wide, and have a 
thickness of 4 mm. Like HI, H2 provides full geometric coverage in the azimuthal 
angle, and has an acceptance between 12° and 65° in the polar angle. The light 
yeild of H2 is on the order of 50-100 photoelectrons per minimum ionizing particle, 
much higher than the light yeild of HI. This allows H2 to give the best dE/dx 
measurement of all the E835 hodoscopes. 

In between HI and H2 is the third hodoscope, H2', which is generally similar to 
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H2 in shape, but is made up of 24 segments rather than 32. These segments are 
40.8 cm long and 4 mm thick, and are located at a distance of 7 cm from the beam 
axis. Like the other two hodoscopes, they provide complete azimuthal coverage, 
and coverage in the polar angle 9 between 9° and 65°. The cracks between the 
segments of H2' are deliberately not aligned with the cracks in HI and H2 in order 
to reduce the leakage of particles. It was added primarily to improve the charge 
veto for neutral triggers, and to improve the dE/dx measurement. 

A forward veto counter also exists which consists of a holed disk placed per- 
pendicularly to the beam pipe as the end cap of the inner detector cylinder. It is 
made up of 8 trapezoidal scintillators of 2 mm thickness, which form an annulus. It 
provides full azimuthal coverage, and coverage in polar angle 6 of between 2° and 
12°. It is used only as a charged particle veto in the forward acceptance region. 

All three hodoscopes, HI, H2, and H2', as well as the forward veto counter 
FV, are used to trigger on charged particles, to act as a veto on neutral triggers. 
For example, the coincidence between HI and H2 generates a first level trigger for 
charged events, and the coincidence between HI, H2' and FV generates a veto signal 
for neutral events. The hodoscopes are also used to measure dE/dx. Finally, they 
form the first piece of the electron identfication algorithm known as electron weight 
which will be described later. 

3.3.2 The Straw Chambers 

After the hodoscopes, the next components of the E835 inner detector system are 
the two straw chambers. These are cylindrical chambers placed with the antiproton 
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beam along their axis, and are comprised of aluminized mylar straws which act as 
proportional drift tubes. These straw chambers are used to determine the azimuthal 
angle of charged particles passing through the inner detector. They give full 
coverage in the azimuthal, and coverage between 15° and 58° in the polar angle 
for the first chamber and between 15° and 65° for the second. Each of the two 
straw chambers is made of two layers of straws, which are staggered azimuthally 
with respect to each other in order to resolve left-right ambiguity. Each layer of the 
chamber is made up of 64 straws. The inner straw chamber has a radius of 54 mm 
and is labelled STR1, and the outer straw chamber has a radius of 120 mm and 
is labelled STR2. The length along the beam axis of STR1 is 182 mm, while the 
length of STR2 is 414 mm. 

The straws are designed to have a low mass to minimize multiple scattering and 
photon conversions, and to have fine granularity to limit occupancy and increase the 
azimuthal angle resolution. Their thickness is 0.11% of a radiation length at a polar 
angle of 90°. These drift tubes are self supporting between two grooved flanges that 
allow gas to flow continously through the tubes. The gas in the tubes is a mixture 
of Argon, C 4 H 10 , and (OCH 3 ) 2 CH 2 in the ration 82:15:3. The tubes themselves are 
made of mylar with a thickness of 80/xm, and are coated on the inside with a layer 
of aluminum approximately 1000 atoms thick; these form the cathodes of the drift 
tubes. The anode of each drift tube is made up of a gold plated tungsten wire placed 
along its axis. These wires had diameters of 5.0-5.4 mm in STR1, and 11.1-12.1 mm 
in STR2. They are crimped at the end of each tube to gold plated copper pins. The 
voltage between the tubes and the wires was 1320 V and 1530 V for the inner and 
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outer chambers, respectively. The drift velocity at these operating voltages was ~ 
40 /xm/ns. 

The readout electronics of the straw chambers are designed to withstand high 
rates in the chambers. They use a custom analog bipolar integrated circuit, with 
an ASD-8B chip and an 8 channel amplifier-shaper-discriminator with fast peak- 
ing time (6-7 ns) and good double pulse resolution (25 ns) to avoid pile-up. The 
straw electronics have a signal amplitude of about 20 mV/fC, and the total power 
dissipation per channel is 23 mW. The front end electronics are mounted on the 
downstream flange of each chamber to minimize oscillations and pickup, and due 
to limited space only Surface Mounting Devices (SMDs) were used. Signals from 
the straw chambers are sent to the E835 counting room, which is located in the 
AP50 building directly above the Antiproton Accumulator. In the counting room 
the signals are processed by 32-channel LRS multihit Time-to-Digital (TDC) 3377 
converters used in common-stop mode. 

The particle detection efficiency of a single straw goes from ~ 100% in the 
vicinity of the wire, to ~ 80% close to the aluminum cathode surface. The measured 
efficiency of track reconstruction with at least two layers of straws is 97%, with an 
efficiency of about ~ 90% per layer. Using both chambers, the angular resolution 
of a track is approximately 9 mrad. The straw chambers are described fully in a 
dedicated paper published in Nucl. Inst. Meth A |102j . The straw chambers were 
not used for the analysis in this thesis, and the azimuthal angles of the e + e~ pairs 
were determined using the Central Calorimeter (CCAL), which is described in a 
later section. 
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3.3.3 The Scintillating Fiber Tracker 

The final elements of the E835 inner detector system are made up of two scintillating 
fiber trackers. The purpose of these trackers is to provide a measurement of polar 
angle 9 for charged particles. The detectors are made of two concentric layers of 
scintillating fibers wound around two coaxial cylindrical supports. The inner tracker, 
SciFl, has 240 fibers per layer, gives full azimuthal coverage, and has a coverage of 
between 15° and 55° in the polar angle 9. The outer tracker, SciF2, has 430 fibers 
per layer, gives full azimuthal coverage, and has a coverage of between 15° and 65° 
in the polar angle 9. The radii of the two cylidrical supports around which each 
of the trackers are wound are 85.0 mm and 92.0 mm for SciFl and 144.0 mm and 
150.6 mm for SciF2. 

The scintillation light from these fibers is detected by solid state photosensitive 
devices called Visible Light Photon Counters (VLPCs). These VLPCs were chosen 
because of their very high quantum efficiency (~ 70% for 550 nm photons). The 
fibers themselves are positioned onto the cylinders in a series of machined U-shaped 
grooves of pitches 1.10/1.19 mm, and 1.10/1.15 mm for the inner and outer SciFl 
and SciF2 support cylinders, respectively. These grooves are machined so that their 
depth varies linearly with the azimuthal angle in order that the fiber can overlap 
itself after one turn without any change in polar angle 9. The starting azimuthal 
angle of each fiber is offset from that of the neighboring fibers in order that the 
fibers do not overlap in any way as they are drawn axially away from the cylinder. 
The fibers are aluminized at one end to increase the light yield and reduce signal 
dependence on azimuthal position. On the other end they are thermally spliced to 
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clear fibers which are 4 m long, and which bring the light to the VLPCs, which are 
kept in a cyrostat at a temperature of 6.5° K. 

The signals from the VLPCs are amplified by QPA02 cards and then sent to 
discriminator-OR-splitter modules which provide an analog and a digital output for 
each input channel, together with the digital OR of all inputs. The analog signal 
is then sent to an Analog- To-Digital (ADC) converter, while the digital signal is 
sent to a Time-To-Digital (TDC) converter. The outer scintillating fiber tracker 
signals were grouped together into 19 bundles of adjacent fibers, and the digital OR 
of the signals from each of these bundles is sent both to a TDC and to the first-level 
trigger logic of the experiment. The scintillating fiber tracker has an efficiency which 
is better than 99% on average in the angular region between 15° and 50°, and better 
than 90% in the region between 50° and 65°. The intrinsic angular resolution of the 
fibers is (0.7 ± 0.1) mrad. 

The typical signal generated by a track through one fiber is 180 mV high and 80 
ns wide, corresponding to a collected charge of ~ 0.2 nC. During calibration, the one 
p.e. equivalent in ADC counts was measured using an LED test for each channel 
in the final readout configuration. The pulse charge in ADC counts generated by 
a minimum ionizing particle was obtained by studying a high statisics sample of 
hadronic tracks 10 3 events/fiber). 

The detection efficiency of the scintillating fiber tracker was measured by using 
the e + e~ tracks from J/ip and ip' decays and pp elastic scattering events. For each 
track an associated hit was sought in the scintillating fiber detector about a given 
software threshold (~ 0.2 m.i.p.) and within a polar angle window of ±50 mrad. The 
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results showed almost 100% efficiency, with the exception of the backward region 
between 50° and 65° in polar angle 9, where there is less redundancy as each track 
intercepts fewer fibers. 

One background which is particularly important for E835 is e + e~ pairs with a 
small opening angle, which are generated by photon conversion or by Dalitz decays 
of neutral pions, and which simulate single tracks. The scintillating fiber tracker 
can be used to help identify such e + e~ pairs in two ways; one using pulse height 
and the other using granularity. In the first case, when the opening angle of the 
e + e~ pair is so small that just one set of adjacent hit fibers is produced, the energy 
deposit is likely to be big. In the second case, when the pair separation is large, an 
extra set of adjacent hit fibers is produced. 

The intrinsic time resolution of the scintillating fiber tracker was evaluated by 
selecting tracks which hit two fibers belonging to adjacent bundles, and thus were 
read out by two separate TDC channels. The rms time resolution is the standard 
deviation divided by \/2 of the distribution in time of the signals of tracks which 
registered in two bundles; it was determined to be approximately 3.5 ns, due mostly 
to the decay time of the scintillator. Further information on the scintillating fiber 
tracker may be found in Refs. |l()3j |l()4j . The scintillating fiber tracker was not 
used in the analysis in this thesis for e + e~ polar angle determination. This was 
done using the CCAL, which will be described in a later section. 
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3.3.4 The Cerenkov Counter 

The E835 Cerenkov Counter is used to identify the lightest charged particles, elec- 
trons and positrons, in a background of much more numerous charged hadrons. This 
counter works on the principle that Cerenkov light is emitted when a charged parti- 
cle passes through a medium with a velocity greater than the velocity of light (c/n) 
in the medium. Thus a particle must have a velocity (5 greater than 1/n where n 
is the index of refraction of the medium. Since the energy of a particle moving at 
speed P is given by: 



E = 7 mc 2 = mC " (3.23) 



the condition for Cerenkov light being produced is 



mc 2 



E > j— = (3.24) 

Since this threshold energy for the production of Cerenkov light is proportional 
to the mass of the particle, and there are more than two orders of magnitude between 
the mass of the electron/positron and that of the next heaviest particle, the pion, 
it is possible to choose a medium with an appropriate index of refraction n which 
can discriminate between e + and e~ and all other charged states over a large range 
of energies. The E835 Cerenkov counter uses two gases as its Cerenkov media. It is 
divided into two cells in polar angle 9, with different gases in each cell to optimize 
the electron detection efficiency and the differentiation between electrons and the 
lightest charged hadrons, the 7r ± s. 
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The overall structure of the Cerenkov detector is that of a cylindrical shell with 



an inner radius of 17 cm and an outer radius of 59 cm. It is designed to completely 
enclose the inner detector within its inner radius, and extends past it along the 
beam axis in the forward direction. The Cerenkov detector has full coverage in the 
azimuthal angle <fi and coverage from 15° to 65° in the polar angle 9. The detector is 
divided internally into eight wedges, each with two separate gas-tight cells separated 
by a septum which is positioned inside the Cerenkov Counter at a polar angle of 
34° or 38°, depending on the distance from the beam axis. Each of the two cells is 
divided interally into eight wedges, each subtending an azimuthal angle of 45°. A 
cross section of one of these wedges is shown in Fig. 3.10. 



Figure 3.10: Schematics of the Cerenkov counter, showing one octant and its mirrors 
and angles. The septum at 34°/38° divides each chamber into its forward and 
backward-angle cells. The spherical, plane and ellipsoidal mirrors are shown in red. 
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Each of the wedges in the small 6 cell is filled with CO2 gas, which has an index 
of refraction n = 1.00041. This corresponds to a Cerenkov energy threshold for 
charged pions of 4.875 GeV, which is greater than the largest pion energy observable 
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by E835. The large 9 cell is filled with either Freon 12 (CF 2 C1 2 ), which has an index 
of refraction of 1.00108, corresponding to an Cerenkov energy threshold for charged 
pions of 3.003 GeV, or with Freon 13 (CF 3 C1), which has an index of refraction of 
1.0072, corresponding to a pion energy threshold of 3.677 GeV. Freon 13 was always 
used in the year 2000 run. The Cerenkov light which is emitted by charged particles 
travelling above the threshold velocity forms a cone around the direction of travel 
of the particle. The half angle 9 C of this cone is given by the relation: 

9 C = cos~ l (l/n(3) (3.25) 
For (3 « 1 this gives a 6 C of 1.64° for C0 2 , 2.17° for Freon 13, and 2.66° for Freon 

12. 

The 8 wedges in the large 9 cell are each equipped with one plane and one 
spherical mirror to reflect Cerenkov light into a photomultiplier. The wedges in 
the small 9 cell are each equipped with an ellipsoidal mirror which has as its two 
foci the interaction point and a photomuliplier for that wedge. The arrangement 
of the mirrors and photomultipliers for both cells are shown in Fig. 3.10. Both 
the small 9 mirrors and the overall mechanical structure of the Cerenkov Counter 
are constructed of carbon fiber epoxy composites covered with a layer of plexiglass 
to improve their reflectance. This was done to minimize the overall weight of the 
structure and the amount of material traversed. The large 9 mirrors, however, 
were made of glass. In these large 9 cells light coming from the direction of the 
interaction region was reflected twice, first by the spherical mirror and then onto 
the plane mirror and from there into the photomultiplier. In both the geometry 
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of the large 9 and small 6 cells, light coming from the direction of the interaction 
region, which is approximately a cube with side 5 mm, is focused by the mirrors 
onto the photomultiplier windows. 

The Cerenkov Counter mirrors are coated with a layer of aluminum which is 
approximately 100 nm thick, and are furthermore protected against oxidation by 
thin films which are MgF 2 on the ellipsoidal mirrors and the plane mirrors, and 
Si0 2 on the spherical mirrors. The thickness of the carbon- fiber /plexiglass mirrors 
was 4.3 mm; and each had a weight of 2 kg. The spherical glass mirrors each had 
a thickness of 3.0 mm and a weight of ~ 1 kg. The plane glass mirrors were 1 mm 
thick. The Cerenkov Counter had an overall light collection efficiency of 0.84-0.90 
for the C0 2 chamber, and an efficiency of 0.84-0.98 for the large 9 chamber filled 
with Freon 13. For Freon 12 this efficiency dropped to the range 0.75-0.88. The 
reflectance of these mirrors was found to be approximately constant for wavelengths 
of light ranging from ~200-300 nm. 

The signals from the Cerenkov photomultipliers were directly amplified by a 
factor of approximately 10. These outputs were split in two, with one half sent to 
the trigger logic, and the other to an ADC to record the integrated charge. A signal 
from at least one of the Cerenkov wedges to the trigger is required for an event 
to be tagged as an electron candidate. The efficiency of this tagging was done by 
examining samples of the reactions pp — > J /ip — > e + e~ and pp — > J/tp — > e + e~7 and 
comparing the numbers of events with at least one electron tagged (le) with the 
number of events with both electrons tagged (2e). From the ratio r = 2e/le between 
these, the detection probability for a single electron is given by e = 2r/ (1 +r) which 
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was found to be (98.1 ± 0.5)% for both the direct J /if) decay and the \c2 decay. 

The average number of photoelectrons per electron track was also measured 
using data from the reactions pp — ► J/ip — > e + e _ and pp — > Xc2 — > J/^l — * 
e + e~7. The Cerenkov signal charge distribution was used to estimate the light 
yield. The average charge may be converted to an average number of photoelectrons 
provided one knows the photomultiplier gain, or the number of ADC counts per 
photoelectron. This was determined by measuring the pedestal and dark current 
separation for each photomultiplier; the separation between the dark current peak 
and the pedestal corresponds to one photoelectron. For the small 9 cell the average 
number of photoelectrons per electron track was found to be dependent on the polar 
angle 9 of the electron, ranging from 19.8 ± 2.3 photoelectrons for electron tracks in 
the 9 = 15° - 20° region, to 11.8 ±1.1 photoelectrons for tracks in the 9 = 32° - 34° 
region. For the large 9 cell, the average number of photoelectrons per track showed 
a much lower variation in polar angle. It ranged from 8.2 ± 1.6 in the 9 = 34° — 38° 
region to 6.8 ± 1.1 in the 9 = 58° — 65° region. A fuller description of the Cerenkov 
detector may be found in Refs. |105j |106| . 

3.3.5 The Central Calorimeter 

The Central Calorimeter (CCAL) is one of the most important elements in the E835 
detector system. It is composed of a cylindrical array of 1280 lead-glass Cerenkov 
counters, and is used to measure the energy and position of photons, electrons and 
positrons. It has full coverage in the azimuthal angle and coverage in the polar 
angle 9 between 10.6° and 70.0°. The main features of the detector are described 
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below; more details may be found in a number of refereneces |lU7j |108j lO^j jllOj . 

The 1280 lead-glass blocks which make up the CCAL are shaped and oriented to 
form a projective geometry in which each block is pointed towards the interaction 
region. This cylindrical array of blocks may be grouped into rings and wedges. A 
CCAL ring is made up of 64 blocks positioned at the same polar angle 6, while 
each of the 64 CCAL wedges covers 5.625° in azimuth, and is made up of 20 blocks 
positioned at the same azimuthal angle 0. A diagram of the CCAL showing its ring- 
and-wedge structure may been seen in Figure 3.11. Each CCAL wedge is contained 
in a light-weight stainless steel container with exterior surfaces which are 1.735 mm 
thick, while the separators between counters between each block in the wedge are 
0.254 mm thick stainless steel. The thicker material separating the wedges leads to 
2% of the surface being inactive in the azimuthal direction; in the polar direction it 
is only 0.5% due to the thinner separations. A detailed description of the orientation 
of the CCAL blocks may be found in Table 3.2 PH] [TOT] . 

The CCAL blocks are made of a Schott F2 type lead glass. They have an index 
of refraction of 1.651 for light with wavelength 404.7 nm, and a density of 3.61 
g/cm 3 . Their composition by weight is 42.2% lead, 29.5% oxygen, 21.4% silicon, 
4.2% potassium, 2.3% sodium, and 0.15% arsenic. The glass has a radiation length 
of 3.141 cm, and its transmittance through 10 cm varies by wavelength from 95.5% 
for light in the 385-394 nm region to 99.4% for light in the 585-594 nm region. For 
light with smaller wavelengths the transmittance is significantly lower; it is 56.9% 
for the 335-344 nm wavelength region. 

A Hammamatsu photomultiplier was glued to the end of each block to collect 
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Central 




Distance from 


Ring 


Length 





A6 


Interaction Region 




f cml 


( deeree i 


( deeree i 


( CIYl) 


1 


37.80 


67.387 


5.226 


72.44 


2 


38.65 


62.259 


5.031 


75.87 


3 


39.88 


57.342 


4.803 


80.07 


4 


41.50 


52.664 


4.552 


85.08 


5 


43.54 


48.246 


4.284 


90.96 


6 


46.06 


44.101 


4.007 


97.79 


7 


48.98 


40.234 


3.728 


105.62 


8 


50.00 


36.644 


3.451 


114.54 


9 


50.00 


33.327 


3.183 


124.66 


10 


50.00 


30.273 


2.925 


136.07 


11 


50.00 


27.472 


2.679 


148.89 


12 


50.00 


24.908 


2.449 


163.26 


13 


50.00 


22.567 


2.233 


179.34 


14 


50.00 


20.434 


2.033 


197.28 


15 


50.00 


18.493 


1.848 


197.29 


16 


50.00 


16.730 


1.678 


197.29 


17 


50.00 


15.130 


1.552 


197.30 


18 


50.00 


13.679 


1.380 


197.30 


19 


50.00 


12.364 


1.250 


197.30 


20 


50.00 


11.174 


1.131 


197.30 



Table 3.2: Geometrical characteristics of the CCAL blocks. The 64 blocks contained 
in a ring share the same features. 
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Figure 3.11: View of a ring (left) and two opposite wedges (right) of the central 
calorimeter. 

the Cerenkov light produced by electromagnetic showers. These photomultipliers 
were chosen for efficient photoelectron collection rather than fast timing, and also 
for their relative insensitivity to magnetic fields. Photomultipliers with four differ- 
ent diameters were used in different rings of the CCAL in order to maximize the 
efficiency of light collection from the different sizes of blocks. 

The output of the photomultipliers is first sent to a summer box, where 5% of 
the signal from each photomultiplier was diverted to the CCAL summer circuit to 
be used for triggering purposes. The remaining signal was sent over a 310 ns coaxial 
RG-58 delay cable to an electronic board that reshapes the signal to fit into a 100 
ns gate for ADC integration. The long cable allows the trigger electronics time to 
make a decision about the signal. The reshaping is required in part because the 
long cable disperses and attenuates the pulse tail to over 600 ns. This was done 
by means of the Splitter-Shaper-Discriminator Circuit (SSD), or simply, the shaper, 
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whose schematic is presented in Figure 3.12. There are three versions of the shaper 
circuit for rings 1-16, rings 17-18, and rings 19-20. By narrowing the pulses, the 
charge may be collected in the 100 ns ADC gate. After ADC integration, a portion 
of the reshaped signal is sent to a discriminator, followed by a TDC. The short ADC 
gate and the TDC time information are needed to reduce the event contamination 
due to pile-up. 

As photons and electrons entering the CCAL are likely to produce a signal in 
more than one block, signals from the CCAL are analyzed by a software procedure 
called a clusterizer, the purpose of which is to determine the energy and position of 
a photon or electron that generates a signal in a group of neighboring blocks known 
as a cluster. The clusterizer first searches for cluster seeds, whcih are defined as 
groups of 3 x 3 CCAL blocks where the central block has an energy deposit which 
is greater than the deposit in any of its eight neighbors. The central energy deposit 
must also exceed a set seed threshold, generally 25 MeV. The sum of all energy 
deposits in the 3x3 grid is called the cluster energy, and this must also exceed a 
set threshold to be considered a cluster, generally 50 MeV. These thresholds may 
be lowered to 5 and 20 MeV for the seed and cluster thresholds respectively for 
analyses which are particularly dependent on the detection of low energy photons. 
Clusters may be grouped into three categories: isolated, shared, and split. 

Isolated Clusters 

Isolated clusters are defined as those which do not have a second seed block within 
a 5 x 5 grid of blocks around the original seed. For an isolated cluster, the position 
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Figure 3.12: Schematic diagram of an SSD channel for the shaping of the input 
signal from a CCAL block. 
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of the incoming particle is calculated as the energy-weighted average of the centers 
of the blocks in the 3x3 cluster grid. Defining (x, y) as the cluster coordinates in 
units of blocks, in the (6>(wedge),0(ring)) directions, where the origin is the center 
of the seed block, we have: 



where Ei is the energy of the ith block, fa is the fraction of the energy of the ith 
block taken by the jth cluster, and Xi = —1, 0, 1 and y; L = —1, 0, 1 are the coordinates 
of the centers of the nine blocks of the cluster in block units. 

In order to take into account the energy which is lost by the cracks between 
the blocks, a correction then is made which takes into account block and crack 
positions, and the electromagnetic shower profile, which is described by a function 
whose parameters were measured with a test beam at the Brookhaven National 
Laboratory. The correction to the shower position is computed as follows: 



x — ^i=ifiEiXi/Ti i=1 fiEi 



(3.26) 



(3.27) 



x 1 = N x [A w a w (l - e~ x/a ™) + B w b w (l - e"^)] 



(3.28) 



N x = 0.5/[A w a w (l - e 



'■ 5/a ™) + B w b w (l - e~°- 5/b ™)} 



(3.29) 
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y' = N y [A r a r {\ - e~ y/ar ) + B r b r (l - e"^)] (3.30) 



N y = 0.5/[A.o r (l - e- ' 5 ^) + - e- - 5 ^)] (3.31) 

where x'{y') is the distance in the wedge (ring) direction, in block units, from the 
center of the seed block. The constants A w , A r , a w , a r , b w , and b r were measured 
empirically with e + e~ decays of the J/ip resonance and are listed in Table 3.3. The 
effectiveness of this correction can be seen in Fig. 3.13. 



A r 


724.4 


a r 


0.03208 


A w 


706.5 


d w 


0.03969 


B r 


123.6 


b r 


0.1860 


B w 


102.6 


b w 


0.1715 


a 


0.0614 


Q 


7.367 


c h 


0.0857 


c h 


19.690 


D 1 


0.14736 


d\ 


48.908 


D 2 


0.15935 


d 2 


12.761 



Table 3.3: Constants used in calculating the position and energy of CCAL blocks. 



The corrected energy of the cluster is given by: 



E corr - E meas /[{1 - C h{l) e-^ C W) x (1 - D ie ~^ - D 2 e-^)} (3.32) 
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where x*,y* are the distances from the cluster center to the closest edges of the 
seed block and the constants are listed in Table 3.3. The ring faces are staggered, and 
Ch(ci) corresponds to the high(low) 6 block edge. The effectiveness of this correction 
can be seen in Fig. 3.13, where the dip in the ratio of measured to expected energies 
in the region around a CCAL crack is removed and there are fewer misidentified 
events after the correction is applied. 




Figure 3.13: Ratio between the measured energy of the cluster and the energy 
predicted by the 2 body kinematics in J/ip — > e + e~ events, before (top) and after 
(bottom) energy correction, as a function of the distance from a CCAL crack. 
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Shared Clusters 

Shared clusters occur when a second seed is found in a 5 x 5 grid centered on the seed 
of the original cluster. When this occurs, the clusterizer assumes that two particles 
have impacted near each other in the CCAL, and it generates two clusters, with one 
3x3 cluster around each seed. As some of the blocks in each 3x3 cluster will 
overlap, the energy of those blocks which belong to both clusters is shared between 
the two. The amount of energy which is assigned to each cluster from the shared 
blocks is determined by an iterative process which includes a method for correcting 
the position and the energy of the two clusters. This method is an extension of the 
position and energy determination algorithm described above for isolated clusters. 

There is also the possibility that a cluster may be due to two adjacent showers 
but may not have a second seed in the neighboring blocks. We then have to examine 
other features indicating that it has likely been formed by two different showers so 
close to each other so as not to create two distinguishable seeds. One such case 
occurs when a neutral pion decays into two photons of approximately equal energies 
in the laboratory frame, also known as symmetric decay, in which case the opening 
angle between the two photons will be a minimum. This minimum opening angle 
can be as low as ~ 1.5 block widths at the upper energy range for pions which can 
be produced in E835. As this separation is too small to produce shared clusters a 
different clusterizing procedure must be used. This involves the cluster mass, which 
is defined as 
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Figure 3.14: Distribution of the cluster mass for electrons and positrons using data 
from J/ip decays (dashed), and from photons from a 7r°7r° sample (solid). The large 
peak at M c i > 100 MeV is due to coalesced photons from the symmetric decay of 
the 7T°. 



where the sum runs over the 25 blocks of the 5x5 grid around the seed, is the 
energy deposited in the i-th block, and pi is equal to Erfi, r\ being the unit vector 
from the interaction region to to the center of the i-th block. This gives a measure 
of the invariant mass of the cluster, which should be minimal (ideally < 0.5 MeV) 
for clusters produced by single photons or electrons. By examining electron and 
positron showers from J/ip — > e + e~ decays, it can be shown that the cluster mass 
of a single electron shower in the CCAL does not exceed ~ 100 MeV, and the same 




(3.33) 
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Figure 3.15: Distribution of the ir° asymmetry (data), defined as (Ei — E 2 )/(E 1 + 
E 2 ), where E\ and E 2 are the energies of the two photons. The distribution is 
plotted without the split clusters (shaded) and including the split clusters (open). 
The distributions are expected to be uniform. 
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is generally true for single photon showers. This is illustrated in Fig. 3.14, where 
the cluster mass distributions of electron clusters from the reaction J/ip — > e + e~ 
are compared with cluster mass distributions from neutral pions from the reaction 
pp — > 7r°7r°, where the higher cluster mass (peak of ~ 140 MeV) from symmetric 
pion decay is clearly observable. 

These split clusters have seeds which are too close together to be resolved with 
the shared cluster algorithm, but use instead a slightly different method in which 
the second seed is chosen as the block with the largest deposit among the four blocks 
at the corner of the original seed. The clusterizer then proceeds to calculate a total 
energy and position for each block using a similar technique to that used for shared 
clusters. Figure 3.15 shows the effectiveness of the split cluster algorithm. The 
7T° asymmetry, defined as (E 1 — E 2 )/(E 1 + E 2 ), is expected to be uniform in the 
data over the range of the asymmetry from to 1. As can be seen from the figure, 
however, the unformity of this distribution is recovered only by including the shared 
clusters. 

By taking into account cluster shape and correcting for cracks and split/shared 
clusters, the CCAL provides an angular resolution of cr e m 6 mrad in the polar angle 
9 and cr^ ~ 11 mrad in the azimuthal angle <fi. The angular resolution is measured 
using a clean sample of electrons and positrons from J ' /ip e + e~ decays, which are 
selected using the hodoscopes and the Cerenkov counter. The energy resolution is 
measured using a similar J/i[) — > e + e~ sample, where one can predict the energy of 
each electron from its direction. The rms energy resolution is found to be: 
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<?e(E) 
E 



s/E(GeV) 



6% 



+ 1.4% . 



(3.34) 



The CCAL is calibrated with both e + e~ events from J /if) and if)' decays, and 
with 7r°7r° events. 7r°7r° calibration is generally superior because the yield of 7r°7r° 
is copious for all pp center of mass energies unlike e + e~ which is only copious at 
J /if> and if)' energies, and also because the gammas from these decays populate all 
counters of the E835 detector. 

Using 7r°7r° events, a clean sample is selected for each stack using the following 
criteria: given four photon 7r°7r° events the combinations of two gammas to form each 
7T° are taken by minimizing the value of a/ (2A6) 2 + (A0) 2 , where A(f> = tx— |0i — <f> 2 \ 
and AO = 9 lpred - 9 lmeas , where 9 lmeas < 9 2meas and 9 lpred is obtained from 6» 2meas 
by kinematics. Here 9 and <f> are the laboratory frame polar and azimuthal angles. 
Once these two two-photon combinations have been chosen to make up the two 7r°s 
in the sample, further criteria are applied. These include: 

• Requiring minimum energy and seed thresholds for all clusters. For energy 
and vertex calibration these must be 50 and 25 MeV respectively, while for timing 
calibration they must be 20 and 5 MeV. 

• Requiring all four photons to be completely detected within the CCAL. 

• Allowing a maximum A(f> (also called acoplanarity) accepted for both n°s. This 
maximum A<f> is 10 mrad for energy and vertex calibration and 15 mrad for timing 
calibration. 

• Allowing a maximum A9 (also called akinematics) accepted for both 7r°s. This 
maximum A9 is 30 mrad for energy and vertex calibration, and 32 mrad for timing 
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calibration. 

• Accepting only a certain range of invariant masses, which for n° — > 27 decays 
may be expressed simply as M T = {2E 1 iE 1 20- — cos6 open )) 1 / 2 where 9 open is the angle 
between the two photons in the lab frame, and E^i, E l2 are their energies in the lab 
frame. For energy calibration \M n — 135 MeV | < 40 is required, while for vertex 
and timing calibration \M n - 135 MeV | < 35 is required 

• Accepting only events with exactly four CCAL clusters. 

Since the event selection depends weakly on energy calibration, it is repeated 
after the energy calibration, and then recalculated. 

The CCAL energy calibration is done by finding the gain constant gi of each 
one of the 1280 CCAL channels. The gain constant is defined as the ratio between 
the energy deposited in the block, and the value determined by the ADC. The gain 
constants are determined by the following procedure. 

The predicted energy of each 7r°, Ej, is calculated from its polar angle using two- 
body kinematics. This is then comapred with the measured energy, Mj , of each one 
of the two 7r°'s, which is given by the sum of the energy deposited in all the blocks 
belonging to the clusters of the two decay photons: 



where j is the index of one of the two pions, % runs over the blocks of the j-th 
pion, and is the ADC value recorded for every involved block. Using Ej , the 
energy calculated by using the measured polar angle of the pion, a x 2 is defined as 



n 




(3.35) 



i=i 
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2 A (Mj - Ej) 2 

x =2^ — 2 ' ( 3 - 36 ) 

3=1 j 

where N is the number of 7r°'s and Oj the estimated uncertainty on Ej. The 
values of the gain constant that minimize the x 2 , calculated analytically, are inserted 
back in Equation 13.361 The procedure is iterated until convergence is reached. The 
gain constants that minimize the \ 2 are given by: 

For the thin, marginal ring 20 which cannot be calibrated with the 7r°7r° sample a 
punch-through method is used instead. It takes into account minimum ionizing par- 
ticles that do not shower and release Cerenkov light in a single block, proportionally 
to the block length. 

Once the above CCAL energy calibration is completed, the measured mass of 
the r\ is compared with the actual mass. A scale correction, typically 1-2%, is then 
applied to all gain constants to match the measured and actual values. This scale 
correction is found to be 2% for electron clusters by examining the mass of the J/ip 
and comparing it to its expected value. A similar 1% correction for gamma clusters 
is found using a GEANT Monte Carlo simulation examining the mass of the r] c . 

Timing information is added to every CCAL cluster by considering the TDC 
channels (used in common stop mode) of the two blocks of the cluster with the 
highest ADC counts. If either of these channels has recorded a hit which is within 
10 ns from the mean event time, then the cluster is classified as in-time. If the 
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recorded hits are not within 10 ns from the mean event time, then the cluster is 
classified as out-of-time. If neither TDC channels recorded a hit at all, then the 
cluster is classified as undetermined. 

The timing algorithm is tested with a 7r°7r° sample from which split clusters have 
been excluded, and for which the cluster energy threshold is set to 20 MeV, rather 
than 50 MeV, in order to examine the algorithm response at low as well as high 
energies. The number of clusters for which time information is determined rises 
sharply as the cluster energy increases from 20 to 50 MeV, since the higher cluster 
energies make it more likely that the TDC hit threshold is exceeded. The fraction 
/ of ADC signals with a corresponding TDC signal can be expressed as a function 
of ADC counts as: 



where slope and thresh are empirically determined parameters. The result of 
the timing algorithm is that over 99.5% of the 7r°7r° clusters with cluster energy 
above 50 MeV are found to be in-time, and the root mean square deviation in time 
for in-time clusters is found to be 1.3 ns. 

When all of the pulses associated with a cluster are below discriminator thresh- 
old, a pulse from a prior interaction will cause it to be mistakenly identified as 
out-of-time. Thus the frequency of out-of-time clusters increases with instantaneous 
luminosity, and this loss of CCAL efficiency must be corrected for later during data 



/ = 0.5 x (erf (slope x (ADC - thresh)) + 1.0), 
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analysis (see Sec. 3.4.3). 

For monitoring and testing purposes, a network of plastic-polymer optical fibers 
was installed to transmit light from a nitrogen laser to the back of each lead glass 
block. The laser produces 3 ns wide pulses with wavelength 337.1 nm which are 
directed to a scintillator that produces light at ~ 430 nm, a wavelength which is 
transmitted efficiently by the fibers. This scintillation light is passed to a rectan- 
gular lucite mixing bar and distributed to optical fibers coupled to the lucite which 
are directed to each wedge of the CCAL. Within the wedge container, the light en- 
counters a secondary mixing bar which further distributes it to 20 more fibers, one 
coupled to the back of each counter. A light-tight aluminum box holds the laser, 
the scintillator, the major mixing bar, and two PIN diodes used to measure the 
intensity of the laser pulse. This system was used throughout the setting up and 
running of E835 to test the CCAL counters and to monitor their gain. 
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3.3.6 The Forward Calorimeter 

A second electromagnetic shower calorimeter known as the forward calorimeter, or 
FCAL, covers the forward region around the beam pipe. It has full coverage in the 
azimuthal angle and coverage in the polar angle 6 between 3° and 11°. 

The FCAL is made up of an array of 144 rectangular SF2 lead-glass blocks of 
three different sizes arranged around the beam pipe as shown in Figure 3.16. The 
SF2 lead-glass used in the FCAL has an index of refraction of 1.673, and a radiation 
length of 2.76 cm. The blocks vary in thickness from 13-21 radiation lengths. Each 
FCAL block is glued to a photomultiplier, which is then connected to electronic 
shaper boards similar to those used by the CCAL. The FCAL signals are sent to 
both ADCs and TDCs. The spatial resolution of the FCAL is approximately 2 cm. 

To obtain an initial calibration for the FCAL, each block was placed in a muon 
beam downstream from the Fermilab MTEST area. For every block, the high voltage 
was adjusted so that the signals from through-going muons, measured using an LRS 
QVT and weighted by the length of each block, were equal. 

After this initial calibration had been done, a sample of six-photon events from 
one of the reactions pp — > 7r°7r°7r — > 67, pp — > n°'K ri — > 67, or pp — > ir°r)r) — > 67, 
taken at the center of mass energy of the Xci resonance, was used to calibrate the 
FCAL. For this calibration 5 of the 6 gammas were required to be detected in the 
CCAL, with the final gamma detected in the FCAL. Two pairs of CCAL clusters 
were required to form an invariant mass consistent with either a 7T° or an 77, and 
the remaining CCAL cluster is paired with the single FCAL cluster. Without using 
the FCAL cluster energy, the event is fit using SQUAW to one of the hypotheses 
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Figure 3.16: Front face (perpendicular to the beam axis) of the forward calorimeter. 



7T 7r 7r°, 7r 7r°?7, or TT°r]r]. The energy resolution of the FCAL, <Je{E) for a photon of 
energy E was found to be: 



(*e{E) 



6% 



+ 4% 



(3.39) 



E y/E{GeV) 

For the analysis in this thesis, the FCAL is used primarily as a veto. The FCAL 
is described more fully in Ref. [Ill] . 
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3.3.7 The Luminosity Monitor 

For an experiment such as E835, which measures excitation curves of resonances 
by varying the center of mass energy of the pp interaction and measuring the cross 
sections at each energy, it is crucial to get an accurate measure of the luminosity. 
To take full advantage of the excellent beam energy resolution afforded by the An- 
tiproton Accumulator, it is necessary that these luminosities be measured with an 
accuracy in the relative luminosity better than 1%, and an accuracy in absolute 
luminosity with an uncertainty of ~2-3%. This is done with the E835 luminos- 
ity monitor, which was built and designed by our Northwestern research group. A 
full description of the previous E760 version of the luminosity monitor is found in 
Ref. j!12j . The present version is described in Ref. [93] . 

The luminosity monitor is based on the detection of proton recoils from low mo- 
mentum transfer (t ~ 0.0077-0.0106 (GeV/c) 2 ) pp elastic scattering events. This 
typically corresponds to a polar angle 9 of the scattered p of ~ 0.20°, and a polar 
angle 9 of the scattered proton of ~ 89.2°. At these tiny momentum transfers, 
the elastic scattering is overwhelmingly Coulombic, allowing accurate absolute nor- 
malization. The recoil proton energies for these scattering events is also small (for 
a = 90° — 9 < 6°,T(p) < 15 MeV for the highest p energies used in E835), so 
that they may be detected in solid state detectors of excellent energy resolution and 
stability. 

The luminosity monitor is made up of three solid state detectors mounted in 
a steel vacuum vessel below the interaction point, at a polar angle of 9 r pa 86.5°. 
One of these detectors lies directly beneath the beam axis, while the other two lie 
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symmetrically on either side. The central detector is movable in the direction of the 
p beam. 

The solid angles of the detectors were precisely defined by a machined tungsten 
mask, with rectangular openings, with each dimension known to ±0.0001 inch. The 
masks for the three detectors had dimensions 

(beam left) : 0.3886" x 0.7889" 
(central) : 0.2752" x 1.7713" 
(beam right) : 0.3884" x 0.7886". 

The detector areas were thus known to better than 0.04%. The mask specifica- 
tions were confirmed with measurements using an a source. The three detectors 
provided not only for threefold redundancy in luminosity monitoring to safeguard 
against detector failure, but they also allow precise and continuous monitoring of 
horizontal displacements of the p beam by a comparison of the event rates in the 
three detectors. A schematic of the luminosity monitor is shown in Figure 3.17, and 
its orientation with respect to the interaction region and attachment to the gas-jet 
target is shown in Figure 3.18. The luminosity monitor is attached to the gas-jet by 
means of a conical vacuum enclosure called the horn. This is suspended just below 
the interaction region perpendicularly to the beam. Attached to the bottom of the 
horn is a pan which holds the solid state detectors, the central one of which may 
be moved along the beam axis from the outside, and the left and right detectors 
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which were fixed at angles qll = 3.496 ± 0.005° and o>r = 3.511 ± 0.005° during the 
running of E835, where a = 90° - 9. 
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Figure 3.17: Schematic of the E835 luminosity monitor. 



The luminosity monitor detectors are silicon surface barrier (500 fxm deep) and 
Si-Li drift (3000 /im deep), with areas ~ 1 cm x5 cm, so that the area exposed by 
each mask is fully contained within each detector. 

The number of protons, N, passing through each one of the luminosity monitor 
detectors is given by 

N = £ [ ^-dtl , (3.40) 

where C is the luminosity, AQ is the solid angle of acceptance, and da/dfl is 
the differential cross section for pp elastic scattering. The detectors' silicon surface 
barriers are sufficiently thick (> 500 /im) to completely stop the recoil protons, 
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Figure 3.18: The E835 luminosity monitor as connected to the gas-jet target. 



which have an energy of just a few MeV. 

The absolute luminosity can be calculated by using the relation, 



_ N 

(da/dt)(dt/dQ)(dQ) (3 ' 41) 



where N is the number of proton recoils counted, (da/dt) is the known differential 
cross section( [HE] ) and dfl is the solid angle subtended by the detector. 
The differential cross-section for the pp elastic scattering is given by: 



da/dt = \F c (t)e iS + F n {t)\ 2 (3.42) 
The Coulomb amplitude F c (t) and phase 5 are exactly calculable. The nuclear 
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amplitude is given by the relation: 

F n (t) = a T (p + t)e- b ^ 2 /(AV^h) (3.43) 

Thus, in order to determine absolute luminosity, three parameters must be de- 
termined. These are the total cross section o^, the 'slope' parameter b of forward 
nuclear scattering, and the parameter p = Re/(0)/Im/(0). 

In its E760 version, the luminosity monitor was designed to also determine the 
parameters cr^, b, and p for pp elastic scattering in the energy region of interest. 
Absolute luminosity for all E835 measurements could then be obtained by using 
these parameters. The elastic scattering differential cross sections da/d\t\, versus 
squared momentum transfer, \t\, were measured using the movable detector with 
better than ±0.5% precision in E760 |114j . the predecessor experiment to E835, 
for six values of the p beam momentum. These cross sections are shown in Fig. 
3.19. Analysis of these cross-sections led to the determination of the cry, b, and p 
parameters at various antiproton momenta. The cross sections shown were fitted 
with smooth curves, resulting in the following parametrization: 



a T = 34.48 + 89. 7/T - 70 , b = 13.64 - 0.2p, p = -0.12 + 0.03p (3.44) 
where a? is measured in mb, and b in (GeV/c)~ 2 . 

These parameters, which are valid for lab momentum p from 2-8 GeV, are used 
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to calculate da/dt for any beam momentum at momentum transfer squared, \t\, 
corresponding to the acceptance of the luminosity monitor detector. Once da/dt 
is known, the number of protons N incident on the luminosity monitor leads to an 
absolute luminosity according to equation 3.41. 
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Figure 3.19: Measured pp differential cross section, da/dt, versus squared momen- 
tum transfer, \t\, at six values of the p-beam momentum (shifted vertically by an 
amount Act to distinguish them on the plot). Fits to the data are shown. 

An example of the proton recoil spectra obtained by one of the luminosity mon- 
itor detectors plotted in a logarithmic scale is shown in Fig. 3.20. iV is determined 
from the spectra by taking the number of counts in the signal region and subtract- 
ing the background, which is well approximated by an exponential (which appears 




100 



as a straight line in the logarithmic plot). The signal region is defined separately 
for each center of mass energy. A luminosity measurement is carried out every 2 
minutes during the running of E835, so that a 120 second average of instantaneous 
luminosity is available during data taking. 




5 6 7 8 
Recoil proton energy (MeV) 



Figure 3.20: Typical spectrum recorded by one of the three solid state detectors of 
the luminosity monitor. 



The luminosity detector response was highly stable throughout the year 2000 
run. As can be seen from Fig. 3.20, the smooth exponential background occurs at 
the level of ~ 2 — 4%, and it can be reliably subtracted from the spectrum. The 
typical count rates in the detectors during data taking were of the order of 10-20 Hz. 
Statistical error in a typical luminosity measurement was very small, of the order 
of 0.3% for 100 nb _1 of data. The dominant source of error was the error in the 
differential cross sections paramatrization of equation 3.44, which was estimated to 
be < ±2.5% [TT3] . 
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The FORTRAN code used to read the counts from the luminosity monitor and 
calculate a luminosity off-line is found in Appendix B. The algorithm used in this 
software to determine the background level samples 21 channels above and 21 chan- 
nels below the signal window. A fit is then made to an exponential and this is 
subtracted from the total number of counts detected in the signal window. A sim- 
plified version of the luminosity monitor code was also used for on-line monitoring, 
in order to give readings of the instantaneous luminosity every two minutes. In this 
version, the readings from the 21 background channels were averaged for each side 
of the window, and a straight line was drawn between them which was used as the 
background level. The on-line code also gave readings of the left-right beam asym- 
metry: by comparing relative rates in the "beam-left" and "beam-right" detectors, 
it was possible to determine whether the p beam was displaced with respect to the 
nominal axis. This was important because during running it was found that the 
beam could occaisionally undergo horizontal shifts of several millimeters perpendic- 
ular to the beam axis. These shifts lead to errors in the luminosity measurement 
due to shadowing of the intersection volume by the aperture at the entrance to the 
luminosity monitor horn. Displacement of the beam from its central orbit causes 
asymmetry in the two fixed detectors. For orbit displacements of magnitude < 1.5 
mm, the asymmetry (L — R)/(L + R) is calculated to be less than 2%, and no cor- 
rections were applied to the average luminosity. For displacements of magnitudes 
between 1.5 mm and 3.0 mm, the corrections to the luminosity were calculated to 
be of the order 10%, but these corrections were generally not used; instead the beam 
orbit was adjusted back to its central value if the asymmetry was found to be greater 
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than 2%. 

The detectors were regularly calibrated over the course of the experiment with 
a 244 Cm alpha source, which is inserted pneumatically into the horn between the 
interaction region and the detector pan. The typical energy resolution for the de- 
tectors was found to be ~ 40 keV for 5.742 MeV alphas. Before the pan assembly 
was attached to the horn, the surface area of each detector exposed by the mask 
was determined by mounting a calibrated (±0.4%) 241 Am alpha source supplied by 
the National Bureau of Standards. The solid angles subtended by the areas of the 
detectors exposed by the masks were then determined by surveying the distance 
between the effective pp interaction point and the detector surfaces; this was mea- 
sured to be 149.60 ± 0.03 cm. The best determination of the recoil angle a for each 
detector was obtained from the centroid of the recoil peak. 



3.4 The Trigger System 

The E835 detector has independent triggers for both charged and neutral final states. 
The charged trigger is primarilly used to identify e + e~ pairs from the electromagnetic 
decay of the J ftp or the ip', and is used in this thesis in the study of the reactions 
pp — > (h c ) — > J/ip+X — > e + e~ +X. The neutral trigger is primarilly used to identify 
reactions in which the final decay products are entirely photons. It is used in this 
thesis for the study of the reaction pp — > (h c ) — > 7^7 — > 37. There is also a specific 
trigger for the completely hadronic reaction pp — > cc — > <fxf) — > K + K~K + K~ , and, 
of course, this trigger was not used for any analysis in this thesis. In addition to 
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these, several monitoring and efficiency study triggers were used in E835. Most 
important of these is the random trigger, also called Random Gate, which is run 
at a frequency of 10 Hz to study accidental signals and event pileup. The random 
gate was used in the analyses in this thesis primarilly in order to correct for dead- 
time due to event pileup. The triggers for both the charged and neutral final states 
are formed independently, and then fanned into a final trigger state along with the 
random gate. A schematic illustrating the general trigger layout is presented in 
Figure 3.21. 
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Figure 3.21: Schematic showing the process of E835 trigger construction. 
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3.4.1 The Neutral Trigger 

The neutral trigger is formed by using the CCAL, complemented by the hodoscopes 
to detect and veto on the presence of charged tracks. Two distinct event categories 
are selected by the neutral trigger; the first one is characterized by two large deposits 
of energy in the CCAL at opposite or almost opposite positions in the azimuthal 
angle 0, and the second one is characterized by a large fraction of the energy in the 
center of mass system being deposited in the CCAL. The triggers PBG1 and PBG3 
belong to the first category and the triggers ETOT-HI and ETOT-LO belong to 
the second category. 

The large number of CCAL blocks (1280) is reduced for trigger purposes to a 
more managable number of signals by associating neighboring blocks into groups. 
The number of CCAL elements is first reduced from 1280 to 160 signals consisting 
of the sum of 9 counters in the same ring (forming what is called a super-wedge) 
overlapping with adjacent sums within the same ring. These super- wedge signals 
are then split with one path being discriminated at ~ 100 MeV and the other 
path being summed in the ring direction to form 40 superblocks. Each superblock 
consists of nine wedges and five rings as shown in Fig. 3.22, with the exception of the 
8 superblocks at high 9, which span only four rings. These 40 superblocks overlap 
each other in order to prevent trigger inefficiencies for particles hitting the peripheral 
elements of a superblock. The superblocks divide the full CCAL into eight parts 
in the azimuthal angle 0, and five parts in the polar angle 9. A small fraction of 
the energy from every block is diverted to the neutral trigger; these signals from all 
of the blocks within a superblock are then summed to form 40 superblock signals, 
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which are in turn used as the inputs of discriminator modules. 




Figure 3.22: View of the front face of the CCAL from the interaction region. Suc- 
cessive superblocks, delimited by tick lines (solid and dashed, alternatively), overlap 
each other to prevent inefficiencies when a shower hits the borders of the superblocks, 
as indicated by E — > [1] and E — > [2]. 

The 40 signals from the superblocks are discriminated and sent to custom pro- 
grammable CAMAC memory-lookup units (MLUs). The Neutral MLU (N-MLU) 
then produces two level-1 triggers known as PBG1 and PBG3 for large invariant 
mass two-body events. Two additional triggers for multiphoton events known as 
ETOT-HI and ETOT-LO are also formed. All these modules are configured via 
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CAMAC in such a way that each output can be any logical combination of the 
incoming signals. The output levels are determined at the trailing edge of a strobe 
signal in order to synchronize the output signals for different triggers. The strobe 
for the N-MLU is based upon having at least two of the super-wedge signals above 
the discriminator threshold. 

The various level-1 CCAL triggers used in the neutral trigger are defined as 
follows: 

• The PBG1 trigger requires a signal in two superblocks whose position is con- 
sistant with two-body kinematics. It is designed for selecting out 77 and e + e~ final 
states. Events from decays such as pp — > 7r°7r° — > 47 are also selected through PBG1 
because of the small opening angle between the two photons from each 7r° decay. 

• The PBG3 trigger has similar requirements to the PBG1 trigger, but with 
a looser requirement on coplanarity: the two hit superblocks do not need to be 
exactly opposite each other in the azimuthal angle 0; rather, the superblock positions 
adjacent to the opposite position are also allowed. PBG3 can thus select events from 
such reactions as pp — > J/^7 — ► e + e~7, or in the neutral case pp — > r} c j — > 37. 

• The ETOT-HI trigger requires that 80% of the energy in the center of mass is 
deposited into the CCAL. It is designed to select multiple photon events. 

• The ETOT-LO trigger requires that 70% of the energy in the center of mass 
is deposited into the CCAL; as such it is a looser version of ETOT-HI. 

Once an event passes one of these CCAL triggers, the signals from the hodoscopes 
are used to determine whether the event is neutral or not. The hodoscope trigger 
signals are defined according to the following logics: 
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• The FV-OR trigger is defined as the logical "OR" between the eight elements 
of the forward veto counter. 

• The HlxH2' trigger is defined as the logical "AND" between any element of 
the HI hodoscope and one of the four corresponding channels of H2'. Four elements 
of H2' are required to competely cover one element of HI as the cracks in the two 
hodoscopes are not aligned. This trigger also allows cases where more than one HI 
element is matched with a signal from one of its corresponding H2' elements. 

• The H2>2 trigger indicates that two or more elements of H2 are hit. 

• Finally the charged veto CV is defined as the logical "OR" between FV-OR 
and HlxH2'; CV = FV-OR x HlxH2'. Its negation is used to declare an event as 
neutral for triggering purposes. 

Once any of the PBG1, PBG3, ETOT-HI, or ETOT-LO triggers are formed in 
the CCAL, they are then associated with the hodoscope trigger logics into the four 
final neutral triggers: 

Neutral-PBGl = PBG1 x HI x H2' x FV-OR = PBG1 x CV 
Neutral-ETOT = ETOT-HI x HI x H2' x FV-OR 
ETOT-NoVeto = ETOT-HI x H2 > 2 

Neutral-ETOT-LO = ETOT-LO x HI x H2' x FV - OR = ETOT-LO x CV 

where the x symbol denotes the logical "AND" operator, and the overline de- 
notes the logical "NOT" operator. 

The Neutral-PBGl and Neutral-ETOT triggers are the ones which are primar- 
illy used to record neutral events, while the ETOT-NoVeto and Neutral-ETOT-LO 
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triggers are mostly used for monitoring and efficiency purposes. The events selected 
by these triggers are then further processed by the software trigger PRUDE (Pro- 
gram Rejecting Unwanted Data Events) for online reconstruction, filtering, sorting, 
and logging. 

The formation time of the PBG triggers is 200 ns and that of the ETOT triggers 
is 155 ns. The formation time of the neutral-veto triggers alone is approximately 
40 ns. The efficiencies of these various triggers are summarized in Table 3.4 for 
various physics channels studied. The hardware for the E835 neutral trigger is 
shown schematically in Fig. 3.23. A detailed description of the neutral trigger may 
be found in Refs. \HM OT- 
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Table 3.4: Efficiencies of the neutral trigger for charmonium (PBG1) and multipho- 
ton (ETOT-H1) final states as a function of center of mass energy. 



Total Energy 
Summer 



1280 
CCAL 
signals 



(2.5%) 



Integrator 



Total Energy 
Discriminator 



20 ring 
sums 



Level I 
summer 



(95%) 



All 1280 
signals 



(2.5%) 



160 ring/wedge 
sums (8 per ring) 



Level II 
summer 



160 ring/wedge 
sums (8 per ring) 



40 super-block 



sums (95%) 



(5%) 



Min Bias 
discriminator 



Shaper 



Integrator 



Super-ring Energy 
Discriminator 



Neutral 
OR 



MLU Strobe 



ETOT-HI 



Neutral MLU 



ADC 


TDC 




Trigger 
Fanout 


t t 



PBG1 



PBG3 



Gatemaster 



Master MLU 



ETOT-LO 



110 



3.4.2 The Charged Trigger 



The charged trigger is designed mainly to accept e + e~ events from the decays of 
charmonium states such as J/ip and ip' |117j . As stated earlier, there is a separate 
trigger to detect 0</> states; it was not used in the analysis in this thesis and will not 
be discussed here. 

The charged trigger selects events with one or two electron tracks and generates 
a neutral veto. The input signals to the charged trigger are the outputs from the 
hodoscopes HI, H2, and H2', the Cerenkov counter, and the scintillating fiber de- 
tector, although the latter is not used in the present analysis. The crucial signal for 
the charged trigger for e + e~ events is that from the Cerenkov counter. A schematic 
of the charged trigger is shown in Figure 3.24. 
le 



2e 
2ch 



NH2=2 



NH2<6 



COPL 



FCH 



u 



C MLUl 



C_MLU2 
C MLU3 



fromN MLU 



w 

H 

< 



e+e- Final 
Trigger 




To DAQ 



Figure 3.24: Schematic of the E835 charged trigger. 



The charged trigger hardware consists of both CAMAC and NIM modules, di- 
vided into three successive stages; the discrimination stage, the single logic stage, 
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and the final trigger stage. In the discrimination stage, the signals from the ho- 
doscopes and the Cerenkov are discriminated by six LRS 4413 discriminators, and 
the 80 output discriminated signals are sent to the trigger logic and to TDCs. The 
outputs from the scintillating fiber detector are discriminated seperately with 32 
discriminator modules. The second stage of the charged trigger consists of the 
single-logic modules; track reconstruction, multiplicity, coplanarity and also the 00 
logic for the 00 trigger which outputs to its own separate MLU. The neutral-veto 
and forward-veto modules which also use information from the hodoscopes have 
been described in the previous section. 

The track reconstruction module gives the number of charged tracks and the 
number of electrons. A charged track is defined as the coincidence between one 
element of HI and a corresponding element of H2. An electron for trigger purposes 
is defined as a charged track where there is also a coincidence with a signal from the 
corresponding Cerenkov octant. The outputs of the track reconstruction module are 
the signals le, 2e, and 2ch, which represent the identifcation of one electron, two 
electrons and two charged tracks. The track reconstruction logic is done with LRS 
NIM modules in approxmately 70 ns. 

The multiplicity module provides the numbers of H2 scintillators and fiber bun- 
dles which recorded signals; these are labelled NH2 and HSF respectively. The 
multiplicity module is comprised of LRS 4532 logic modules that give output sig- 
nals proportional to the number of active input signals. The multiplicities that are 
selected from this module for the final charged trigger are NH2=2, NH2<4, NH2=3 
or 4, NSF=3 or 4, and the combined multiplicities NH2 + NSF = 7 or 8. The 
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mulitplicity module requires approximately 65 ns to output a signal. 

The coplanarity module requires at least two hits in H2 which are back to back 
within 22.5°, i.e. 180° ± 22.5°. This module uses only the signals coming from the 
H2 hodoscope and requires that if an element of H2 is hit then at least one of the 3 
opposite elements of H2 are also hit. The coplanarity module requires approximately 
50 ns to output a signal. 

Once the track reconstruction, multiplicity and coplanarity signals have been 
formed by the single-logic modules, they are sent to a custom CAMAC programmable 
MLU known as the C-MLU (charged memory-lookup unit). There is also a separate 
0-MLU for the <£><£> trigger signals. The outputs from these MLUs are combined with 
the neutral trigger output from the N-MLU and sent to the master memory-lookup 
unit (M-MLU) . The output of the M-MLU is then strobed into the data acquisition 
system (DAQ), which writes the event. 

The strobe signal for the C-MLU is the logical OR of the H2 hodoscope signals. 
Since the timing of this signal is critical, it is reshaped using an Ortec 934 constant 
fraction discriminator, which leads to a jitter of 2 ns or less. Input signals that 
arrive at least 8 ns before the strobe are successfully latched. To allow for a 5 ns 
jitter, the inputs are set to reach the MLUs 15 ns before the strobe and are 30 ns 
wide. The strobe itself is 10 ns wide. 

The transit time through the single-logic stage of the charged trigger is 96 ns 
for the e + e~ branch. The transit time through the final trigger stage to the output 
of the M-MLU is about 90 ns (including some delay needed for synchronization 
of the neutral logic). The total transit time through the whole trigger is thus 
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about 200 ns. The total efficiency of the charged trigger is e = (90 ±1)%; this is 
the number which is used as the trigger efficiency for all charged analyses in this 
thesis jnn]. This efficiency was calculated by examining a set of ip' — ► e + e~ events 
which were collected using a simplified trigger requiring only one electron track, 
and then checking to see how many times the second electron track passed the full 
charged trigger |118j . 

3.4.3 The Random Gate Trigger 

Because of the high instantaneous luminosity at which E835 is run, one of the main 
sources of inefficiency is the dead-time due to the pile up of different events. The 
event pileup is difficult to deal with because it varies with the instantaneous luminos- 
ity, which, despite the efforts made to limit its fluctuations during the data taking, 
may be different at various center-of-mass energy points. Thus when performing 
an energy scan, it is critically important to estimate the amount of the luminosity 
dependent pileup contamination and correct for it according to the instananeous 
luminosity of each data taking run. This is done primarilly by using information 
recorded through the random gate trigger. 

During the data taking, this random gate trigger is run at a frequency of 10 
Hz with priority over the physical triggers. While it is active the gate of the data 
acquisition is opened and the event signals from all channels of the CCAL, FCAL 
and Charged Veto (hodoscopes) are recorded. Off-line a Monte Carlo simulation is 
run and the efficiency of the detector system is determined by overlapping random 
gate events with generated Monte Carlo events; during this overlap the Monte Carlo 
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content of each ADC and TDC channel of the detectors is summed to the random 
gate content of the same channel, and afterwards the clusterization and event selec- 
tion is performed on the composite Monte Carlo/random gate signal. In this way 
it can be determined how likely a random gate event is to cause a genuine event to 
be lost. That information, combined with the overall rate of random gate events, is 
then used to determine the efficiency 

For charged channels the detection efficiency as a function of instantaneous lu- 
minosity was found to be: 

e « (1 - 0.04 x L inst ) (3.45) 

where the instantaneous luminosity Ci nst is measured in units of 10~ 31 cm _2 s _1 . 

This efficiency was found using the random gate data collected at xo energies, 
overlapped with GEANT Monto Carlo generated pp — > J/ipl e + e~7 events j!19j . 
This correction includes the effect of random gate events on both the e + e~ signals, 
and the 7 signals. 

For all-neutral channels, such as pp — > 7r°7r° — > 47, the luminosity dependent 
correction to the efficiency is greater: 

e « (1 - 0.075 x C inst ) (3.46) 

Fig 3.25 shows the result of the random gate analysis for the charged channel 
described above. The luminosity dependent corrections to the efficiency in the above 
equations are calculated from the slope of the efficiency vs. luminosity plot for that 
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reaction. 

The overall efficiency at each energy point is determined by averaging the Ran- 
dom Gate efficiency of its runs, (e) run , and taking a weighted average of the runs 
according to their integrated luminosity C run : 

a x e = ^ — l — — x ^ £mn x (e) run . (3.47) 

l^run L run rm 

The uncertainty on e is 



* •run 



V 7— x JE t^x^w] 2 - ( 3 - 48 ) 

2-^irun ^run \ mm 



3.5 The Data Acquisition System 

Data acquisition and recording in E835 is done from four independent data ac- 
quisition streams running in parallel: the event data stream, the data stream for 
parameters of the p beam, the luminosity monitor data stream, and the scaler data 
stream to monitor the rates of selected counters and triggers. Fig. 3.26 illustrates 
the schematic of the system. 

The event data stream is based on the DART j!20j data acquisition system, which 
was developed by a collaboration between the Fermilab Online System Department 
and several experiments that operated during the 1996-97 fixed target run. DART 
provides a common system of hardware and software which has been configured 
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Figure 3.25: The probability that an e + e~ event survives the pileup contamination 
versus the instantaneous luminosity, on a run by run basis |119j . The probability was 
calculated using random gate events collected at the Xo energy region. The black 
circles represent the efficiencies calculated using GEANT, and the red triangles 
represent the same efficiencies scaled to account for the difference between GEANT 
event and data events. 
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Figure 3.26: Layout of the hardware of the event data acquisition. The ADC and 
TDC readouts, where the signals first enter the data acquisition system, are located 
in the CAMAC crates on the upper right of the diagram. Here, DYC stands for 
Damn Yankee Controller, DPM's are the Dual Ported Memory units, SGI Indy 
and SGI Indigo are two of the Silicon Graphics computers used for run control and 
monitor display respectively, and PTI are Performance Technologies Interfaces 109J . 
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and extended to meet the specific requirements of the E835 experiment. The beam 
parameter information data stream is received from the FNAL Beams Division via 
ethernet. Both the luminosity monitor stream and the scaler data streams are input 
and received via standard CAMAC. 

Three Silicon Graphics (SGI) computers were used to coordinate, process, filter, 
log, and monitor the data streams. The run-control computer, called SGI-Indy, co- 
ordinates the independent data streams, communicates with the CAMAC branches, 
and performs the slow data acquisition functions. The event data acquisition, which 
consists of event building, filtering, and logging, is performed by an SGI Challenge- 
L computer with 12 processors. Finally, the monitor-display computer is an SGI 
Indigo, which monitors the detectors and displays events. 

3.5.1 The Event Data Acquisition 

The layout of the event data acquisition is illustrated in Fig. 3.26. The readout 
consists of 163 ADC (LRS4300) modules and 114 TDC (LRS3377) modules orga- 
nized in 14 CAMAC crates. They are further arranged into 3 CAMAC branches, 
two parallel and one serial, and addressed by the run-control computer through two 
SCSI 411 Jorway Interfaces. 

Data from the detectors are read out by the ADC and TDC modules, and trans- 
ferred through their front-end ECL ports to the DART-developed Damn Yankee 
Controller (DYC) modules. The DYC is an intermediate data buffer which stores 
the 16-bit input data in a 32-bit wide FIFO (first-in first-out). The DYC adds a 
header word including a word count, an error bit in the case of event overflow, and 
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4 bits of event synchronization. Data is sent out to two pairs of Access Dynamics 
DC2/DM115 modules over two RS-485 cables. 

Each ADC or TDC CAMAC crate is served by a DYC. The DYC responds to 
the request from the readout electronic modules it serves, which are connected in 
a daisy-chain, by issuing a readout enable to the first module, starting the data 
transfer. When the first readout electronic module has finished sending data, it 
enables the following module in the chain by sending a pass signal. The end of the 
data transfer is signified by the receipt of the pass signal from the last module. 

While receiving data, the DYC modules issue a busy signal. The logical "OR" 
of the busy signals from all the DYC modules is sent to the trigger logic to inhibit 
triggers during the data reception. 

The DYCs are connected in a daisy-chain via two RS-485 cables (labelled as 
DART cables in Fig. 3.21) in two groups, each of them connected to one of the 
DC2/DM115 pairs. Data transmission to the DC2/DM115 pairs is controlled by a 
wait signal, established by the destination buffer, and a permit token. On receipt of 
a permit token, and in the absence of a wait signal, the DYC transmits the header 
word and the data from one event over the RS-485 cable. When the DYC has 
completed its transmission, it passes the permit token to the following DYC in the 
daisy-chain. The DYC module is able to simultaneously perform the data reception 
from the readout electronic modules and the data transmission to the DC2/DM115. 

Each DC2 is connected, via the VSBus of the VME crate, to two 32 MB and 
one 8 MB Dual Ported Memory (DPM) MM-6390 modules. Each DC2 fills the two 
32 MB DPMs with the data it receives from the DYCs using a ping-pong algorithm. 
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According to this algorithm the DC2 writes to the ping DPM with exactly N (the 
number which is set in the configuration database) events before it switches to the 
pong DPM. 

A process running on the Challenge, called the gateway process, reads data from 
the ping DPM via PTI(PT-VME64) while the DC2 writes to the pong DPM. When 
all the data from the ping DPM are read, gateway allows the DC2 to write to the 
ping DPM and turns to read out the pong DPM. Both the DC2 and the gateway poll 
each other to see if each has finished using a DPM. The communication between 
the DC2 and the gateway is done using the 8 MB DPM as a mailbox. Neither the 
DC2 nor the gateway interrupt a data transfer while polling. 

The gateway writes the two DC2 sets of data, which comprise N events, to 
two buffers in a reserved part of the memory of the Challenge-L, which has been 
designated as memory to be shared among processes. When finished writing a 
buffer, the gateway releases it to a filter queue before writing into the other buffer 
in the shared memory. An online filter process runs continuously on each CPU of 
the Challenge-L. Each filter queues for a gateway buffer; when one is available, an 
online filter processes the N events. First, the events are built taking data from both 
of the two DC2 data sets. Then, synchronization numbers of the DPMs and the 
DYCs are checked and the data integrity is verified for 0.1% of the events processed. 
Finally, an event header is added which includes a unique event number, date, and 
time of processing and pointers to each DYC header word for later use by the offline 
event reconstruction. Some basic event trigger classification and online analysis bit 
flags are also included. 
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The trigger information is decoded to classify the event. This classification 
determines which online filter analyses to perform, if any, and which logging stream 
the event is sent to if it is not rejected. 

The following online analyses are performed: a simplified CCAL determination 
of energies and angles for electromagnetic showers and formation of invariant mass 
pairs to identify n°, rj, and cc candidates, simple charged track reconstruction with 
electron identification and association with CCAL energy deposits for identifying 
events that fit exclusive e + e~ final states or J/ip + X — > e + e _ + X decay, and 
four-track identification for cfxp events. 

The online filters set analysis bit flags in the event header word and write sum- 
mary event information which is added to the end of the event. The analysis bit 
flags are compared to a list of trigger masks to determine acceptance or rejection of 
the event. 

Accepted events are copied to logging buffers reserved in the shared memory for 
eventual recording to tape and transfer to the monitor-display computer. There 
are three logging streams; which stream an event is written to is based on the 
event trigger classification and the analysis bit flags. Events for which the filter 
determines a large invariant mass and those to be used for calibration purposes are 
written to other logging buffers. When the logging buffers are full (they have a 
maximum size of 64 KB; the average event size is ~ 1 KB), they are released to 
the appropriate logging queues and to the monitor-display computer. There are 
ten logging processes, each one with its own queue and tape drive. Every logging 
process receives buffers from all of the filter processes. Generally, five tape drivers 
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record data in parallel to 8-mm EXABYTE tapes with a 5.0 GB capacity, while 
the other five frivers are left ready to be able to start as soon as one of the tapes 
in the first five drivers becomes full. The switching of tape drivers signifies the end 
of a single run of data taking and the beginning of a new one. Two of the five 
simultaneously running tape drivers were dedicated to neutral events and another 
two were dedicated to charged events; thus there were two charged tapes and two 
neutral tapes for each run. The final tape driver is used for efficiency study events. 

For a characteristic instananeous luminosity of 2.5 x 10 31 cm~ 2 s _1 , the trigger 
rate is ~ 4 KHz. The typical event size from the readout electronics is ~ 1.2 KB. 
Thus the two daisy chains and the gateway typically transfer ~ 4.8 MB per second 
to the online filters. The filters reduce the output event rate to ~ 1.2 KHz, which 
corresponds to writing ~ 1.4 MB per second to the logging streams. Each tape 
drive is capable of 1000 KB per second; the tape drives typically record data at 
a rate of 350, 300, and 250 KB per second for the neutral, charged, and efficiency 
study events, respectively. The resulting lifetime of the data acquisition as a whole is 
typically 93%, at the instantaneous luminosity of 2.5 x 10 31 cm _2 s _1 . The inefficiency 
comes from inhibiting the triggers during the ADC conversion time and from the 
data transfer to the DYC modules. 

3.5.2 Non-Event Data Streams 

The control and monitoring of the FNAL accelerator complex is performed through 
the Accelerator Control Network System, known as ACNET. Accelerator data for 
E835 is acquired by a Fermilab Beams Division ACNET computer, which is inter- 
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faced to the accelerator instrumentation and sends a data stream to the run-control 
computer. This data stream consists of the positions and intensities recorded by 
the BPMs, the longitudinal frequency spectrum of the beam, the gas-jet density, 
pressures within the accelerator, magnet currents, magnetic field strenghts, and the 
online determination of the beam energy, emittances and current. The run-control 
computer writes the accelerator data to disk depending on the status of the run 
control. The online beam momentum and beam current are unpacked and made 
available for monitoring purposes every two minutes. 

The CAMAC readout electronics of the luminosity monitor detector consists of 
a multiplexed ADC converter (EGG ORTEC-AD413A)and a FERAbus Histogram- 
ming Memory (EGG ORTEC-HM413). A self gating ADC channel is used to read 
each of the three silicon detectors. Each channel is read and histogrammed via the 
front-end ECL port by the HM413 module. Every two minutes a process running 
on the run-control computer reads the three data sets from the HM413 module 
via the CAMAC backplane, resets the memory, computes the luminosity by the 
algorithm described in Section 3.3.7, writes the result on disk, and sends it to the 
monitor-display computer. 

The readout of the luminosity monitor data stream shares a CAMAC crate with 
the readout from the scaler data stream. Three LRS4434 scaler modules receive data 
via their front-end ECL ports from the trigger and detector electronics. A Fermilab 
CAMAC pulser (RFD01) module sends a signal to the scalers to update every 10 
seconds. The data stream is sent every ten seconds via the CAMAC backplane to 
the run-control computer which writes it to disk and displays it on a monitor. 
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3.6 The Offline Event Reconstruction 

Tracking of charged particles by the inner detectors is done separately for the r/0 
and r/z projections. The straw chambers are used for the r/0 projection and and 
scintillating fibers together with the nominal position of the interaction vertex for 
the r/z projection. These processes are then associated using CCAL information. 

The azimuthal angle for charged tracks is determined using the signals from 
the straw chambers. First the drift distance is determined within the hit straw and 
a fit to a straight line through the layers is made using the geometrical parameters 
of the straws. For each straw the discriminator threshold is set at one primary 
ionization. 

The drift time, which is measured with 500 ps (rms) resolution by an LRS 
multihit TDC 3377 used in Common Stop mode, is defined as the delay of the straw 
signal with respect to the event time, the latter of which is defined by a strobe 
constructed from CCAL signals. For each straw, a reference time T is measured 
which depends on the propagation delays through the electronics and cables. The 
drift time distributions (T — To) should be uniform for constant drift velocity and 
full efficiency, since the particle density is uniform in the azimuthal angle <fi. Non- 
constant drift velocity is corrected for to obtain the distance X(T — T ) for each 
layer using the relation: 



where N is the total number of tracks, dN/dx (= constant) is the track density, 




(3.49) 
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and X is the perpendicular distance between the particle line of flight and the wire, 
which determines the azimuthal angle <fi and is independent of 9. For each track a 
straight line is fit through the hits to determine with a single-track rms resolution 
of ~ 9 mrad. This is comparable to the angular resolution of the CCAL, which 
was measured to be ~ 6 mrad in the polar angle and ~ 11 mrad in the azimuthal 
angle fTUTj . 

The polar angle 9 is determined for charged tracks using the signals from the 
scintillating fiber detector. A charged particle traversing the scintillating fibers hits 
up to three consecutive fibers in each layer, depending on the polar angle 9. Clusters 
are then defined as sets of adjacent hit fibers; those with more than three fibers are 
split into two or more smaller clusters. The z coordinate of the intersection between 
the track and the layer is taken as the mean of the z coordinates of the fibers in 
a cluster, weighted by the energy deposit in each fiber. This value of z is then 
converted to a 9 coordinate using the nominal coordinates of the interaction point. 
Once this is done, the clusters found in each layer are associated to form tracks, by 
requiring that their polar angles 9 be within a specified range. The clusters for each 
track and the interaction point determined using the CCAL, as previously described, 
are then fit to a straight line. 

A nominal interaction point for the CCAL analysis is required in order to com- 
pute the laboratory angles of the clusters, and thus the momentum of each particle. 
The interaction region is located at the intersection of the beam and the gas jet. 
While the jet location does not change, the beam position and size varies from stack 
to stack. For every stack, the (x, y) coordinate of the center of the interaction region 
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is determined from the 7r°7r° sample. The desired precision of 0.1 mm is obtained 
by using 10,000 events per stack. 

The laboratory coordinate system is defined as follows: The z axis is the beam 
direction, the z coordinate of the gas-jet center is zero, the intersection point of the 
(pointing) counter axes of the CCAL is (x,y,z) = (0, 0, z c ), and the x — z plane 
passes through the axes of all of the counters of CCAL wedge number 1. 

For a two-body reaction, trigonometry gives the following formula for the acopla- 
narity, A<p = 7r — \<f>i — fal, for a vertex at (x , y , 0) where <f> is the azimuthal angle 
of one of the 7r°s and R is the counter radius in the x — y plane: 

XnSind) — VnCOSd) 

tan(A(f) 2) = — — (3.50) 

R — x cos(p — y sin(p 

The distribution of A0 for each stack is fit to the above equation to obtain xo 
and ?/o- As these parameters depend only on the 7r° directions, (xo, yo) is relatively 
insensitive to the energy calibration, and we determine it using the gi (described 
previously in the CCAL section) from the previous stack, prior to the energy cali- 
bration, which depends strongly on (x ,y ). 

3.7 Electron Identification 

Proper electron identification is crucially important for most of the analyses per- 
formed in this thesis; in particular the identifcation of electrons from J/tp decay. 
For all analyses of interest to the E835 experiment it is necessary to identify elec- 
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tromagnetic decays of charmonium in the presence of a hadronic background which 
is between 10 8 and 10 11 times larger, requiring not only excellent electron identi- 
fication, but also the ability to distinguish between 7 and n° signals. Individual 
electrons must be distinguished from other "tracks" producing Cerenkov light, in- 
cluding e + e~ pairs due to photon conversion in the beam pipe and from Dalitz 
decays of n°s. 

The detectors that provide useful electron identification information are the HI, 
H2, and H2' scintillator hodoscopes, which measure dE/dx; the scintillating fiber 
detector, which also measures dE/dx and can identify small angle electron pairs; 
the threshold Cerenkov counter, which is designed to count only electrons and can 
identify electron pairs using pulse height; and the CCAL, in which cluster shapes 
are different for single electrons, hadrons, and merged electron pairs. 

Using all the available information, a likelihood ratio has been developed for 
E835 analyses. It is known as the electron weight (EW), which combines signals 
from all these detectors. A second method, known as the k-nearest neighbor rule, is 
used to check the assumptions of the EW method. It is the electron weight method 
which is used primarilly in this thesis to identify electrons. 

The electron weight of any given charged track is calculated as follows. For each 
charged track, characterized by measured quantities x = (xi, ...,x n ), we seek to 
distinguish between the hypothesis e that it is an electron vs. the hypothesis b that 
it is not an electron. The optimal (Neyman-Pearson) method is the construction of 
the likelihood ratio p(x), derived from the probability density functions F e (x) and 
F b (x), under the two hypotheses respectively. For a single track we thus have: 
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P(x) 




(3.51) 



We define P e and P& as the usually-unknown overall efficiencies for electrons e 
and non-electrons b respectively, resulting from the trigger and selection efficiencies. 
It is thus convenient to rewrite the above equation as: 



If P e and Pb were known, the cut p > 1 would give the optimal discrimination 
criterion. Since they are unknown, the reduced likelihood ratio is used, which is 
defined as: 



Since it is desirable to maximize the signal to background ratio rather than to 
minimize the absolute number of misassignments, a suitable minimum value of p' is 
empirically determined, depending on the analysis being performed. This reduced 
likelihood ratio p' is refered to as the electron weight. The measured quantities x 
which are included in the electron weight are: 

• HI pulse height corrected by polar angle 9. 

• H2 pulse height corrected by polar angle 9. 

• H2' pulse height corrected by polar angle 9. 

• Cerenkov counter pulse height corrected by polar angle 9 and mirror photo- 




(3.52) 




(3.53) 
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electron yield, with the correction done independently for all 16 mirrors. 

• sg and s<f,, the CCAL cluster second moments in wedge and ring directions, 
calculated according to the relations: 



Z 3 r , w=1 E(r,w) x (w- w ) 2 

■ s » = (3 ' 54) 



E(r,w)x(r-r ) 2 

where the index r runs over the CCAL rings, the index w runs over the CCAL 
wedges, and E(r,w) is the energy in block (r,w). 

• F 35 , the ratio of the energy deposited in a 3 x 3 block matrix divided by that 
in a 5 x 5 block matrix, both surrounding the highest energy block in the cluster. 

• F 2 4, the ratio defined similarly to the F 35 ratio, but using matrices of sizes 2 
x 2 and 4x4. 

• M d) the cluster mass as defined previously in the discussion of the CCAL (see 
Section 3.3.5). 

For the CCAL variables s^, s$, F 35 , F 24 , and M d , all of the blocks surrounding 
that with the highest energy are used, whether or not the cluster is split. 

To compute the electron weight, / e (x) and /&(x) are taken as the normalized 
products of the measured distributions of each variable, for clean samples of elec- 
tron and non-electron tracks. Although components of the electron weight, such as 
hodoscope signals, are obviously correlated in the presence of a real electron track, 
it is assumed that the electron weight variables x are uncorrelated, in order to fac- 
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tor / e (x) and /&(x) into individual probability distribution functions. The electron 
weight EW is then calculated as: 

^=ftK=nf§=n* (3, 6) 

where the Wi are the individual weights of which the final electron weight is the 
product. 

The individual probability distribution functions are determined from events 
selected by the trigger. The tracks must consist of a CCAL cluster associated with 
at least two out of three hodoscope hits and a Cerenkov signal. For the electron 
probability distribution functions a clean sample of J / ip and Xc2 events is taken with 
consistent event topology and acceptable kinematic fit probability. For the non- 
electron background probability distribution functions events are used from runs 
taken well outside the regions of the known charmonium resonances, specifically 
at energies between 3590 and 3660 MeV (which were taken during the rj' c search 
in the 1996/7 E835 run). The probability distribution functions themselves are 
described in detail in Ref. An illustration of the ability of the electron weight 
to distinguish between e + e~ events and non-electron background events is shown 
in Figure 3.27. The unshaded histogram in this figure shows the electron weight 
distribution of events from the reaction pp — > xi — > J/'i'l ~ > e + e~ r y. The shaded 
histogram shows the electron weight distribution of background events from the 
same runs. Background events were defined as those in which M(e + e~) < 2.8 GeV. 

Since all events of interest have either zero or two electrons in the final state 
the product of electron weights is most often used in E835 analyses which examine 
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log(EWl xEW2) 



Figure 3.27: Comparison of electron weights of e + e~ events and non-electron back- 
ground events for all data taken in the \i region. The dashed line represents the 
electron weight cut at EW1 x EW2 > 1, which was used in the analysis of most of 
the decay channels in this thesis. 



132 



reactions containing the decays J/ip — > e + e~ and ip' — > e + e~. Some of the analyses 
in this thesis also use requirements on the electron weights of individual electron 
tracks in order to give a particularly clean sample. The level of the electron weight 
cut for each particular analysis is done empirically to yield the maximum signal to 
background ratio; examples of such analyses are done in the next chapter for the 
various e + e~ channels. Once a cut level has been established an independent J/i/j 
sample is used to measure the efhency of the electron weight cut, and a sample of rj' c 
search events are used to determine the background rejection power. The efficiencies 
of the various electron weight cuts used in the analyses in this theses are discussed 
in the next chapter. 

3.8 Simulations 

Several simulations are used in the analysis of data from E835. The most important 
of these is the GEANT 3.2 simulation used for simulated data sets and to compute 
efficiencies and response properties for individual detectors. When very large sim- 
ulated data sets are required, mainly for all-neutral channels, events may also be 
simulated using a parametric description of the CCAL response which runs 10-20 
times faster than GEANT and allows for enough event generation to study large 
backgrounds such as the 7r°7r° feeddown background to 77 events. Further discus- 
sion of this technique may be found in Ref. [UH]; for the purposes of this thesis the 
GEANT simulation was used exclusively to compute the efficiencies of various cuts. 
A separate set of kinematic simulation programs was also written in FORTRAN 
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specifically for the analyses performed in this thesis, in order to calculate geometric 
acceptances with large statistics. 



Chapter 4 

Data Selection and Analysis 



4.1 Decay Channels 

The present search for the h c i}Pi) resonance was done by scanning in the range y/s 
= 3523 - 3529 MeV in the following decay channels: 

(1) pp -> -> J/^ + X, J/^ -> e+e- 

(2) -> (^0 -> J/rj) + 7r°, J/^ -> e+e", vr° -> 77 

(3) pp — > J/^ + 7, J/V — > e + e~ 

(4) pp -> (iPx) -> J/^ + 7r°7T , J/ip -> e+e-, 7r°7r° -> 2(77) 

(5) pp -> (^1) -> 7^7 -> 37 

The decay channels (1-4) containing J/ip are the most promising because our 
detector is optimized for e + e~ and photon detection. The first channel J /if) + 

134 
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X includes all pp events which produce a J/ip and any number of gammas with 
an appropriate total energy. It is sensitive to any resonance which could produce 
J/ip as a decay product, and thus has the largest statistics of all the individual 
J/ip containing decay channels. The second channel, the isospin violating reaction 
J/ip + n , is important because it was the channel in which E760 claimed observation 
of the l Pi in its scan of the 3523 - 3529 MeV energy range. The third channel is 
used only as a control channel, since the decay 1 P 1 — > J/ip^f is strictly forbidden 
by C-parity conservation. The fourth channel, 1 P 1 — > J/ip + 2ir° does not violate 
isospin, but we note that the phase space space for this reaction is small. The final 
channel (pp — s- 1 P± — > r) c ~f) was investigated because it is the most promising of all 
decays which do not involve a J ftp. This channel will be discussed separately at the 
end of this chapter. 

The data used in this analysis was taken by experiment E835p (year 2000 run) 
at eighteen energy points in the region 3524.79 < y/s MeV < 3527.29, which is 
approximately ± 1 MeV around the value of \/s = 3526.2 MeV at which the res- 
onance enhancement was reported by the E760 experiment. In order to constrain 
the background, data were also taken at three energies (y/s = 3523.3, 3528.7 and 
3529.1 MeV) which were about 1 MeV outside these limits. The total luminosity 
invested in the present 1 P± search in E835p was 50.5 pb" 1 . This is to be compared 
with 38.9 pb -1 invested in the E835 1 Pi search, and 15.9 pb -1 invested in the E760 
measurement. The typical mass resolution was ^(v^) ~ 0.3 MeV, or FWHM(y / s) 
0.7 MeV, which was due entirely to the momentum spread of the antiproton beam 
which was measured by its Schottky noise spectrum. Table 4.1 lists the parame- 
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ters of the present scan, and those of E835 (year 1997 run) and E760 are listed in 
Table 4.2. Energies shown in these tables are calculated by taking an average over 
an entire stack of the measured energies in each run, weighted by their respective 
luminosities. These weighted averages generally differ from the nominal energies of 
the stack, being higher by ~ 300-500 keV. 

The hardware trigger for the first four decay channels was designed to accept all 
events with a J/ip — > e + e~ decay. It required exactly two charged particle hits in 
the corresponding elements of the inner and outer scintillator hodoscopes, with at 
least one particle identified as an electron by the associated Cerenkov signal. Inde- 
pendently, it was required that two large energy clusters in the CCAL be separated 
by more than 90° in azimuth, and have an invariant mass greater than 2.2 GeV. 
Additional cuts were imposed in off-line event selection. 

The total e + e" data set from all 50.5 pb~l taken in the E835p l P\ search is 
shown in the top panel of Figure 4.1. As can be seen from that figure, additional 
cuts are required to identify J/ip, as well as n°, cleanly. These cuts are discussed 
below. 

4.2 Analysis Cuts 

The cuts used for this analysis fall into three main categories; cuts for electron iden- 
tification, cuts for gamma identification, and probability cuts for overall hypothesis 
testing. 

All channels involving J/ijj require cuts for clean identification of electrons. This 
is done by the use of the electron weight cut. As described in Sec. 3.7, the electron 
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Stack No. 


Nominal \fs 


Measured y/s 


c 


£(corr) 


&beam 




^ivie v j 


^ivie v ) 


l nD ) 


l nD ) 


^ivie v ) 


11-23 


3523.3 


3523.33 


3058.6 


2811.0 


0.390 


11-39 


3524.7 


3524.79 


2033.0 


1820.7 


0.297 


11-21 


3525.1 


3525.17 


3709.6 


3303.2 


0.300 


11-20 


3525.4 


3525.46 


4305.1 


3840.4 


0.309 


11-16 


3525.7 


3525.74 


3316.0 


2807.9 


0.376 


11-37 


3525.8 


3525.88 


1668.1 


1506.5 


0.256 


11-38 


3525.8 


3525.89 


1638.3 


1475.6 


0.288 


11-18 


3525.95 


3526.02 


3742.4 


3353.0 


0.294 


11-15 


3526.2 


3526.21 


3674.5 


3233.2 


0.442 


11-24 


3526.25 


3526.25 


510.5 


464.5 


0.346 


11-22 


3526.2 


3526.28 


2627.9 


2384.6 


0.280 


11-36 


3526.2 


3526.29 


1647.5 


1458.5 


0.268 


11-35 


3526.2 


3526.29 


1488.5 


1333.1 


0.291 


11-41 


3526.2 


3526.30 


1947.3 


1775.3 


0.253 


11-25 


3526.25 


3526.32 


952.7 


875.3 


0.304 


11-26 


3526.25 


3526.42 


2469.5 


2238.0 


0.328 


11-17 


3526.5 


3526.57 


3707.5 


3249.5 


0.291 


11-19 


3526.8 


3526.89 


3036.7 


2720.2 


0.311 


11-40 


3527.2 


3527.29 


1223.3 


1098.7 


0.265 


11-28 


3528.5 


3528.61 


1279.3 


1135.7 


0.364 


11-27 


3528.9 


3529.11 


2467.0 


2233.4 


0.386 


Total 






50503.3 


45118.3 





Table 4.1: Stack numbers, nominal energies, measured center of mass energies, 
luminosities, dead-time corrected luminosities and beam widths for each E835p data 
stack. The dead-time corrected luminosity is calculated for each stack using the 
random gate method described in Sec. 3.4.3. 
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E760 






E835 








Stack 


v ° 


c 


Stack 


^/s(meas) 


c 


£(corr) 


IN UlllUt:! 


^IVie V ) 


l nD ) 


IN UIIlUcI 


^IVie V J 


l nD ) 


l nD ) 


3 


3522.6 


980 


1-61 


3524.6 


3717 


3441 


5 


3523.5 


490 


1-58 


3525.2 


2903 


2768 


2 


3524.0 


783 


1-62 


3525.5 


3532 


3259 


1 


3524.3 


823 


1-55 


3525.7 


3477 


3237 


9 


3525.0 


1041 


1-60 


3525.8 


2976 


2750 


13 


3525.6 


1310 


1-63 


3525.8 


3821 


3509 


18 


3525.9 


885 


1-59 


3525.9 


1196 


1103 


14 


3526.1 


1364 


1-19 


3526.1 


1992 


1841 


20 


3526.1 


980 


1-56 


3526.1 


2291 


2117 


11 


3526.2 


1337 


1-57 


3526.2 


3309 


3041 


22 


3526.2 


876 


1-64 


3526.5 


3234 


2976 


19 


3526.2 


940 


1-51 


3526.5 


976 


904 


16 


3526.3 


1017 


1-54 


3526.6 


649 


601 


21 


3526.5 


911 


1-52 


3526.9 


3094 


2871 


15 


3526.6 


1137 


1-53 


3527.5 


1396 


1310 


17 


3527.2 


1016 


1-50 


3529.1 


2328 


2177 


Total 




15890 






38899 


36064 



Table 4.2: Center of mass energies and luminosities for each E760 data stack. Center 
of mass energies, luminosities, and dead-time corrected luminosities for each E835 
data stack. The dead-time corrected luminosity is calculated for each stack using 
the random gate method described in Sec. 3.4.3. 
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Gold ee Files — All P 1 Events 




all e + e 

88 1 5 events 



2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 

MASS e + e" (GeV) 




MASS 77 (GeV) 

Figure 4.1: All e + e~ events from the 21 stacks taken during the E835p X P\ search. 
The top figure shows the invariant mass of the two electrons for all e + e~ events. 
The J/ip peak is barely visible. The middle figure shows the e + e~ invariant mass 
for all e + e~ + 77 events. The J/ijj peak is observable. The bottom figure shows the 
77 invariant mass for all e + e~ + 77 events. The 7r° peak is clearly visible. In all 
cases the need for further cuts to define the e + e~ and 77 peaks clearly is indicated. 
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weight takes into account pulse heights in the hodoscopes and Cerenkov counters, 
and the transverse energy distributions in the CCAL clusters. 

Several levels of electron weight cuts were examined in this analysis. The results 
are shown in Figure 4.2 (with any number of gammas), Figure 4.3 (with two gam- 
mas), Figure 4.4 (with one gamma), and Figure 4.5 (with four gammas). As can be 
seen from these figures, a clear J/ip peak at 3.1 GeV is visible with the minimum 
electron weight cut generally used in E835 and E835p, i.e. EW1 x EW2 > 1. As can 
be seen from Figures 4.2-4.5, tightening the selection to require EW1 x EW2 > 10 
has the effect of making the J/ ip peak sharper, but at some cost to the total number 
of events. Requiring both electron weights to each be separately greater than one 
is an even stronger cut. This is shown as the third choice of electron weight cuts 
in these figures. As shown later, for most channels, choosing the combined EW(1) 
x EW(2) > 1 proved to be most efficient at giving a clean J/ip sample when used 
together with probability cuts as described later. 

The gamma cuts were based mostly on the number of on-time clusters observed 
in the CCAL which were not associated with one of the two electrons. For the 
J/ip + X channel, no number-of-clusters cuts were imposed, and any number of 
gammas were allowed. The data for this channel is shown in Figure 4.2. For the 
J ftp + 7T°, J/ip + 7, and J/ip + 7r°7r° channels, 2, 1, and 4 gamma clusters were 
required. These data are shown in Figures 4.3, 4.4, and 4.5, respectively. Figure 4.3 
also shows the invariant mass spectrum of the two gammas; the 7r° peak is already 
clearly observable even before additional cuts are made. One additional on-time 
cluster was also allowed in the J/ip + 7r°, J/tp + 7, and J/tp + tt°-k° to account for 
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a possible bremsstrahlung gamma as long as it made an angle of < 10° with one of 
the charged particles, had an energy < 100 MeV, and did not make a n° (i.e. have 
invariant mass within the limits 135 ± 35 MeV) with any of the other gammas. 
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All Events With Electron Weight Cuts 




2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 



M(e+e-) (GeV) 

Figure 4.2: Effect of electron weight cuts on the M(e + e~) distribution for e + e~ 
events with any number of gammas. The tightest cuts (shaded green) are those 
requiring both EW1 and EW2 to be separately greater than one. The effect of the 
tighter cuts is to reduce the number of non- <//?/> events at the lower end of the e + e~ 
invariant mass spectrum. The small number of e + e~ events with an invariant mass 
in the 3.4-3.6 GeV region are ascribable to direct production or the proton form 
factor. 
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Figure 4.3: Effect of electron weight cuts on the M(e + e~) distribution for e + e~ + 77 
events. Top figure shows the e + e~ invariant mass for three choices of electron weight 
cuts. Bottom figure shows the corresponding 77 invariant mass distributions. 
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e + e + 7 Events With Electron Weight Cuts 
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Figure 4.4: Effect of electron weight cuts on the M (e + e ) distribution for e + e + 7 
events. 
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Figure 4.5: Effect of electron weight cuts on the M(e + e ) distribution for e + e +4 / y 
events. 
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4.2.1 Probability Cuts 

Probability distributions were constructed for each final state hypothesis and these 
are shown in Figures 4.6, 4.7, 4.8, and 4.9 for the J/ip + X, J/ip + 7r°, J/ip + 7, and 
J/ip + 7r°7r° channels respectively. As can be seen from the first three figures, the 
probability distribution begins to flatten out at ~ P = 0.05, and thus in the final 
selection events with P > 5% were accepted for each of these three decay channels. 
The probability plot is also shown for the J/ipir ir (4 cluster) channel, but in this 
case only two events survive the P > 5% cut. 

The effect of various combinations of electron weight and probability cuts are 
shown in Figures 4.10-14 for the J/ip + X, J/ip + 7r°, J/ip + 7, and J/ip + 7r°7r° 
channels, respectively. Figures 4.11 - 4.14 illustrate the fact that when combined 
with probability cuts, EW(1) x EW(2) > 10 does no better in isolating J/ip or 7r° 
than EW(1) x EW(2) > 1, but costs some in net counts. In the J/ip + X case, 
however, the tighter electron weight cut (EW1 x EW2 > 10) is found to give a 
significantly cleaner J/ip mass spectrum, and was therefore used. 

The probability distributions for the J/ip+X and J/ip+j decays were calculated 
using a five constraint fit; the constraints being the total energy and momentum of 
the event and the invariant mass of the e + e~ pair which makes up the J/ip. For the 
J/ip + 7T° channel, the probability distributions had six constraints, the additional 
constraint being the invariant mass of the two gammas which make up the n°. 
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Figure 4.6: Probability distribution for J/ip + X for all e + e events with EW1 x 
EW2 > 10. Events with P(J/ip + X) > 0.05 (dotted line) were used in the final 
analysis. 
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Figure 4.7: Probability distribution for J/ip + n for all e + e~ + 77 events with EW1 
x EW2 > I. Events with P(J/ip + n°) > 0.05 (dotted line) were used in the final 
analysis. 
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Figure 4.8: Probability distribution for J/ijj + 7 for all e + e~ + 7 events with EW1 
x EW2 > I. Events with P(J/ip + 7) > 0.05 (dotted line) were used in the final 
analysis. 
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Figure 4.9: Probability distribution for the channel J/ip + 7r°7r° for all e + e + 47 
events with EWlx EW2 > I. Events with P{J/i) + vr°7r°) > 0.05 (dotted line) were 
used in the final analysis. Only two events survive this cut. 
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Figure 4.10: Effect of various electron weight and probability cuts on the M(e + e~) 
distribution for J/tp + X events. The top plot shows the effect on the e + e~ invariant 
mass distribution for two choices of electron weight cuts with p(J/ip + X) > 0.01. 
The bottom plot shows the same for p(J/ip + X) > 0.05. The tightest cuts (EW1 
x EW2 > 10, p(J/ip + X) > 0.05) show the cleanest J/ip peak and are used in the 
final analysis. In addition, a J/ip mass cut is made to remove events outside the 
region 2.9 < M(e + e~) < 3.3 GeV (dotted lines). 
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Figure 4.11: Effect of various electron weight and probability cuts on the M(e + e~) 
distribution for J/ip + ir events. The top plot shows the effect on the e + e~ invariant 
mass distribution for two choices of electron weight cuts with probability P(J/tp + 
7T°) > 0.01. The bottom plot shows the same for probability P(J/tp + 7r°) > 0.05. 
The cuts (EW1 x EW2 > 1, P(J/ip + vr ) > 0.05) are used in the final analysis. 
J/ip mass cuts are made to remove events outside the region 2.9 < M(e + e~) <3.3 
GeV (dotted lines). 
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Figure 4.12: Effect of various electron weight and probability cuts on the M(7i72) 
distribution for J/ip + tt° events. The top plot shows the effect on the 77 invariant 
mass distribution for two choices of electron weight cuts with P(J/tp + 7r°) > 0.01. 
The bottom plot shows the same for P(J/ip + tt°) > 0.05. 7r° mass cuts are made 
to remove events outside the region 90 < M (gamma'j) < 180 MeV (dotted lines). 
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Figure 4.13: Effect of various electron weight and probability cuts on the M(e + e~) 
distributions for J/ip+j events. The top plot shows M(e + e~) for P(J/i/;+ r y > 0.01), 
and the bottom plot for P(J/ip + 7 > 0.05). A final mass cut is made to remove 
events outside the region 2.9 < M(e + e~) < 3.3 GeV (dotted lines). 
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Figure 4.14: Effect of various electron weight and probability cuts on the M(e + e ) 
distribution for J/ip + 7r°7r° events. The top plot shows the M(e + e~) spectrum for 
P(J/V + vr°7T ) > 0.01 and the bottom plot for P(J/ip + 7r°7r°) > 0.05. As in the 
previous channels, a final mass cut (dotted line) is made to remove events outside 
the region 2.9 < M(e+e~) < 3.3 GeV. 
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Finally, motivated mainly by the tails present in the M(e + e~) distribution for 
the J ftp + X decay channel, a mass cut of 2.9 GeV < m e + e - < 3.3 GeV was made 
to define J/ip. Although this cut is almost unnecessary in the J/ip + n and J/tp + l 
channels, where it removes very few events, we also apply it for these decay channels 
for the sake of consistency with the J/ip+X decay channel. For the J/tp+ir° channel, 
an additional invariant mass cut was used for n°. As indicated in Figure 4.12, the 
cut used was 90 MeV < m 77 < 180 MeV. The same mass cuts were imposed on 
the J/ip + 7r°7r° channel (Figure 4.14), although they had no effect, since only two 
events remained in this channel once the electron weight and 5% probability cuts 
were used. 

4.3 Excitation Curves 

Having decided on the best cuts for event selection in each decay channel, the data 
for each stack was analyzed separately, in order to plot the excitation curves. In 
Table 4.3 the counts obtained in each final state are listed for each E835p energy 
point, along with the corresponding luminosities and beam widths. Table 4.4 shows 
the energy by energy event listing of the E760 experiment in the J/ip + tt° and 
J/ipl channels, as published in Phys. Rev. Lett.[?]. Table 4.5 shows the same 
for the J/ip + 7T° channel for the E835 (1997) run, using the same final cuts as 
were used for the current E835p analysis. Cross sections at each energy point were 
obtained by dividing the event counts iVj at each energy by the luminosity Ci and 
the overall efficiency e assumed to be constant over the < 5 MeV total range of 
y/s scanned. Account was taken of dead time corrections due to instantaneous 
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luminosity variations using the random gate corrections described in Sec. 3.4.3. 
The efficiencies of the various cuts on the four e + e~ channels are shown in Table 
4.6. In all trigger efficiency of 90% is taken, as descibed in Sec. 3.4.2. The 

efficiency of the probability cuts for each final state was obtained by GEANT MC 
simulation, as were the acceptances for each decay channel. For the J/ip+X channel, 
which may have several final states, the efficiency was calculated as the average of 
those for the various individual final states. For example, the acceptance of the two 
electrons was 57.4% for the J/ip + n°, 56.3% for the J/ip + n channel, and 59.4% for 
the J/ip + 7r°7r° channel. The J/ip + X acceptance of 57.7% is the average of these 
three. The electron weight and probability cut efficiencies for the J/ip + X channel 
were calculated similarly to obtain the overall efficiency for the J/ip + X channel. 
The overall efficiency (including acceptance) was as follows e(J/ip + X) = 39.1%, 
e(J/ip + 7t°) = 17.2%, e( J/ip + 7) = 25.6%, and e(J/ip + tt° + tt°) = 5.3%. 

Using the total efficiencies, the cross sections <7j = Ni/de) are obtained. These 
correspond to a (final state) x B(J/ip — > e + e~), with final state = J/ip + X, J/tp + j, 
or J/ip + 7r°7r°. These cross sections are listed for the J/ip + X, J ftp + ir° and 
J/ip + 7 channels in Table 4.7. Luminosities listed in Table 4.7 are after they 
have been corrected for dead time using the random gate. The table does not 
list (i(J/ip + 7r°7r°), since only two counts were observed in the entire scan. For 
J/ip + X the results from the tight EW selection are listed. The total efficiency for 
the J/ip + 7T° channel in the E835 (year 1997 run) data was 16.2%. For the E760 
data, the total efficiencies (as reported in Ref. [?]) are 32.4% and 44.6% for the 
J/ip + 7T° and J/ip + 7 channels respectively. 



158 



The excitation curves corresponding to the cross sections in Table 4.7 for the 
decay channels J/tp + X, J/tp + n° and + 7 are shown in Figures 4.15-4.18. For 
the sake of clarity the data for points with Ay's < 0.15 MeV have been combined 
in all plots, although the subsequent determination of upper limits for branching 
ratios was done with the uncombined data. Figures 4.15 shows the excitation curve 
for the J/ip + X channel, using the two different choices of electron weight cuts 
described before. In neither case any enhancement is observed. A constant cross 
section fit gives a = 26.6 ± 1.4 pb (x 2 /d.o.f. = 1.30) for the loose EW selection, and 
a = 25.2 ± 1.4 pb (x 2 / d -°- f - = i- 75 ) for the tight EW selection. 

Figure 4.16 shows the excitation curve for the J/ip + ir° channel, in which E760 
found evidence for the observation of 1 P\ with a mass of 3526 ± 0.15 ± 0.2 MeV. 
As can be seen from the figure, we do not find any evidence for an enhancement in 
this channel anywhere in the scan region. The data cross sections are fitted with 
a constant of 13.6 ± 1.4 pb with an excellent xVd-o.f. of 0.31. Figure 4.17 shows 
a comparison of the excitation curves for the J/ip + 7r° channel for E835p, E835, 
and E760 data. The excitation curve for the 1 P\ for the forbidden decay channel 
J/ip + 7 is shown in Figure 4.17. Again, no enhancement is observed, and J/ip + 7 
cross sections are fitted with a constant 8.8 ± 0.9 pb (x 2 /d.o.f. = 0.66). Figure 4.18 
shows a comparison of the excitation curves for the J/ip + tt° channel for E835p, 
E835, and E760 data. 
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Table 4.3: Number of events observed for each E835p data stack in the channels 
J/ip + X, J/ip + 7, J/ip + 7r°7r°, and J/ip + 7r°. Both tight and loose electron weight 
cuts for the J/ip + X channel are shown. Data stacks with At/s < 0.15 MeV are 
combined for later analysis as follows: stacks marked * are combined to a single 
point at 3525.95 MeV, stacks marked * are combined to a single point at 3526.22 
MeV, and stacks marked * are combined to a single point at 3526.50 MeV. 
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Table 4.4: Number of events observed for each E760 data stack for the channels 
J /ib + 7 and J/ip + 7r°. The cross sections have been obtained by using overall 
efficiencies e = 32.4% for the J/ip + n° channel and e = 44.6% for the J/ip + 7 
channel, calculated as the products of the acceptances and analysis efficiencies given 
in the published paper. 
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J/lp + 7T° 


(Pb) 


3524.6 


3441 


6 


10.8 


3525.2 


2768 


9 


20.1 


3525.5 


3259 


6 


11.4 


3525.7 


3237 


17 


32.4 


3525.8 


2750 


10 


22.4 


3525.8 


3509 


1 


1.8 


3525.9 


1103 


1 


5.6 


oozo.l 


Zll 1 


4 


11. I 


3526.2 


3041 


3 


6.1 


3526.5 


2976 


4 


8.3 


3526.5 


904 


2 


13.7 


3526.6 


601 





0.0 


3526.9 


2871 


4 


8.6 


3527.5 


1310 


4 


18.8 


3529.1 


2177 


3 


8.5 


Total 


36064 


74 




Efficiency 






16.2% 



Table 4.5: Number of events observed for each E835 data stack for the J/ip + ir { 
channel. The overall efficiency of this channel for the E835 data is 16.2%. 
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Channel 


Cut 


Efficiency 


Common 


Trigger 


90% 




J/^ mass cut 


98% 


(All but J/tjj + X) 


EW1 x EW2 > 1 


r\ r~ fry 

95% 


J/^J + X 


p(J/4> + X)> 5% 


90% 


loose 


EW1 x EW2 > 10 


92% 


tight 


EW1 > 1, EW2 > 1 


85% 




Acceptance 


57.7% 




Total (loose) 


42.1% 




Total (tight) 


39.1% 


J /if) + 7T° 


P (J/lp + 7T°) > 5% 


71% 




iV( 7 ) = 2 


93% 




7T° mass cut 


95% 




Acceptance 


32% 




lotai 


1 / .z/o 


J/ip + 7 


p( J/^ + 7) > 5% 


77.5% 




7V( 7 ) = 1 


93.7% 




Acceptance 


42.2% 




Total 


25.6% 




p(J/ip + 7T 7r°) > 5% 


56.8% 




7V( 7 ) = 4 


92.3% 




Acceptance 


12.0% 




Total 


5.3% 



Table 4.6: Efficiencies of cuts used in decay channels containin; 



■g W- 
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-v/sfmeas) 

V V / 

(MeV) 


££(corr) 
(nb- 1 ) 


EiV 
X 


for 

7T° 


J/ip+ 

I T 1 

7 


a (pb) 

+x 


for 

+ 7T° 


JM(-> e+e-) 
+7 


3523.33 
3524.79 


2811.0 
1820.7 


22 
18 


6 
3 


6 
8 


on n+4.3 

20. 01^3 

or 0+6.8 
zo - z -5.9 


12.4_^5 

n «+7.3 

9-6-6.1 


o q+4.6 

o-^.o 

17 - 2 -5.8 


3525.17 


3303.2 


21 


6 


6 


16.3™ 


10.61H 


7 i +3.9 
' • i -2.6 


3525.46 


3840.4 


43 


11 


11 


28.6™ 


16.71H 


11.2™ 


3525.74 
3525.95 


2807.9 
6335.1 


26 
68 


8 

15 


8 
8 


on o+3.6 

— 3.0 

22.4™ 


17.1™, 

—4.2 

10.9™ 


9.5™ 

—2.5 

4.9™ 


3526.22 


11524.5 


130 


31 


30 


28.6™ 


15-2™ 


9.711;? 


3526.50 


5487.5 


61 


12 


11 


29.1111 


12.5™ 


8-4™ 


3526.89 


2720.2 


24 


6 


4 


22.6™ 


12.811;? 


5-7™ 


3527.29 
3528.61 
3529.11 


1098.7 
1135.7 
2233.4 


6 

14 
17 


2 
2 

5 


2 
3 
4 


14.0™ 

31-5™ 
19.5™ 


10.6™ 
10.2™ 
13.0™ 


7 i +8.0 

10.3™ 
7-0™ 


Efficiency 


(e) 




39.1% 


17.2% 


25.6% 



Table 4.7: Cross sections a (final state) x B(J/ip — > e + e~) observed for the 
finalstates J/ip + X, J/ip + n , and J/ip + 'j. Data taken at Ay's < 0.15 MeV have 
been combined. Errors in the individual cross sections are statistical (Poisson). 
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o 

CJ 



50 



£ 45 



40 



55 - 



50 - 



25 



20 



1 5 
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J/V + X Excitation Curve 



Tight cuts — 450 Events 
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Loose cuts — 512 Events 
Average a = 26.6 ± 1 .4 pb 
X 2 /d.o.f. = 1 .30 
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Figure 4.15: Efficiency corrected excitation curve for the J/ip+X channel with loose 
electron weight cuts (EWlxEW2>10) (closed circles) and tight electron weight 
cuts (EW1>1, EW2>1) (open circles - slightly displaced). The cross sections are 
consistent despite different cuts. Stacks with energies less than 0.15 MeV apart have 
been combined. 
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Figure 4.16: Efficiency corrected excitation curve for the J/ip + ir channel from the 
E835p data. Stacks with energies less than 0.15 MeV apart have been combined. 
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Figure 4.17: Efficiency corrected excitation curve for the J/ip + 7 channel from 
the E835p data. Stacks with energies within 0.15 MeV of each other have been 
combined. 
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30 



20 



10 
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Figure 4.18: Comparison of excitation curves for the J/ip + 7r° channel for E835p 
(top plot), E835 (second plot), and E760 (third plot). 



4.4 Backgrounds 
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As Table 4.7 and Figs. 4.15 - 18 show, there is no apparent enhancement in the 
excitation curves for the three reactions, and good fits are obtained assuming con- 
stant cross sections: a(J/i^ + X) = 25.2 ± 1.4 pb, a(J/ip + tt°) = 13.6 ± 1.4 pb, and 
a(J/t/j+^f) = 8.8±0.9pb. The observed sum of cross sections a(J/ip + n°) +a(J/^p + 
7) = 22.4 ± 1.9 pb essentially exhausts the observed a(J/ip + X) = 25.2 ± 1.4 pb. 

These observations suggest that the cross sections observed in all three decay 
channels arise primarilly from continuum production and tails of distant resonances. 
As illustrated by Table 4.8 and Fig. 4.19, an estimate of continuum <r(J/tp + vr°) is 
provided by its value of 11.9 ± 1.9 pb which we obtain from our scan in the region 
y/s = 3415 ±15 MeV containing the Xo resonance which is forbidden to decay in this 
channel by C-parity conservation. To obtain these results, cuts identical to those 
used in the l P\ scan were used. From the known branching ratios for xj —> + 7 
we calculate the total contribution at 3526.2 MeV from the tails of the xj resonances 
as a x (J/ip + r y) = 5.5 ±0.9 pb, of which ~ 2.9 pb comes from the tail of X2, ~ 2.1 pb 
comes from the tail of xi, an d ~ 0.5 pb comes from the tail of xo- It is estimated by 
Monte Carlo simulation that an additional contribution of 2.2 ± 0.3 pb arises from 
J/ip + n° events in which a photon is lost, bringing the total estimated background 
to a{J/ , i])+ r y) = 7.7±1.0 pb. Thus, both <r( J/ip + ir°) and a(J/ip+^f) observed in the 
present measurements are indeed well accounted for by the 'background' processes 
described above, as is cr(J/ip + X) which is made up primarilly by the sum of the 
two channels. This leaves little room for an observable resonance enhancement in 
either J/ip ± 7r° or J/ip ± X channels. 
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Stack 


^(meas) 


C 


£(corr) 


Events 


cr(J/ip + vr°) 


Number 


(MeV) 


(nb- 1 ) 


(nb- 1 ) 


observed 


(Pb) 


II-47-b 


3339.5 


659 


585 


1 


10.9 


II-47-a 


3365.0 


1424 


1287 





0.0 


11-10 


3384.4 


1630 


1464 


2 


8.6 


11-31 


3384.8 


3371 


3038 


4 


8.4 


11-42 


3392.0 


1431 


1276 





0.0 


II-34-b 


3400.1 


1483 


1352 


2 


9.4 


II-12-b 


3404.7 


76 


68 





0.0 


II-14-b 


3406.1 


2481 


2243 


5 


14.2 


II-32-b 


3409.1 


1134 


1057 





0.0 


TT 1 O 

11-12-a 


3410.2 


1613 


1461 


2 


8.7 


11-33 


3413.8 


2959 


2616 


4 


9.7 


II-9 


3415.0 


2352 


2131 


8 


23.9 


11-13 


3415.9 


2801 


2498 


5 


12.7 


II-7 


3417.9 


1456 


1345 


2 


9.5 


II-32-a 


3422.1 


2512 


2250 


3 


8.5 


11-11 


3425.9 


1820 


1631 


6 


23.4 


II-34-a 


3430.1 


1438 


1260 





0.0 


II-8 


3469.9 


2513 


2302 


1 


2.8 


Total 




33153 


29864 


45 


< a >= 9.6 ± 1.4 


Efficiency 


(e) 








15.7±0.4% 



Table 4.8: Number of events observed and corresponding cross sections for each 
E835 data stack in the \o region for the J/ip + 7r° channel. The 'a' and 'b' labels 
denote stacks which were split into two energies by decelerating the p beam between 
runs. The overall efficiency of this channel for the E835 data in this energy region 
is 15.7% ± 0.4%. 
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X stacks — p p — > J / 1^ + 7T 
Vs = 3400 - 3430 MeV 
Average a = (1 1 .9 ± 1 .9) pb 
X 2 /d.o.f. = 0.85 
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X 2 /d.o.f. = 0.56 
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Figure 4.19: Efficiency corrected excitation curve for the J/i/j + tt channel from the 
E835p data for the xo energy range. The average cross section is 9.6 ± 1.4 pb over 
the whole Xo region (solid line), and 11.9 ± 1.9 pb in the central Xo region (dashed 
line), as compared to 13.6 ± 1.4 pb in the 1 P i energy range. The total efficiency of 
this channel at y/s = 3415 MeV is calculated using GEANT to be 15.7% ± 0.4%, as 
opposed to 17.2 ± 0.4% for ^ = 3526 MeV. 
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4.5 Upper Limits 

As shown in Figures 4.15-4.18, there is no apparent enhancement in the excitation 
curves for any of the three reactions pp — > J/x/j + X, pp — > J/tpn° or pp — > J/tpj, 
and fits assuming constant cross sections yield good x 2 per degree of freedom. 

In order to quantify the above observations in terms of upper limits for the 
identification of /i c ( 1 P 1 ), the data for pp — > ( 1 Pi) — > J ftp + X, pp — > ( 1 Pi) — > 
J/ip + n°, and — > ( 1 Pi) — > J/ip + r y (only for control purposes) have been analyzed 
by the maximum likelihood method for fit to a constant background and a Breit- 
Wigner resonance. The likelihood function was defined as 

i=i 

where 

^(Vs) = A • e • + y fi(y/s)<7Bw(Vs)dy/s\ (4.2) 

The expected number of counts in the ith scan point is Vi(y/s) with center of 
mass energy y/s, while the observed number of events is rij. £j is the integrated 
luminosity and e is the overall efficiency for the observation of the final state. The 
integral provides for the convolution of the Breit Wigner cross section with the 
(Gaussian) center of mass beam energy distribution function fi(\/s) as determined 
by the Schottky noise spectrum. For the Breit Wigner resonance, three different 
choices of width T = 0.5 MeV, T = 1.0 MeV, and T = 2.0 MeV were assumed. The 
data at = 3523.33 MeV, 3528.61 MeV, and 3529.11 MeV were used to determine 
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the level of the constant backgrounds which was 11.1 ±3. 1 pb for J/ip + ir , 22.1±3.0 
pb for J/ip + X, and 8.2 ± 2.3 pb for J/tp'y. These fixed backgrounds were therefore 
assumed and the mass of the resonance, Mr, was varied in 0.1 MeV steps. For each 
trial the 90% confidence upper limit for the product branching ratios B in x B out 
were determined by the maximum likelihood fit, where B in = B(pp — > R) and 
B out = B(R — > finalstate) (where finalstate can be one of J/ip + X(J/ip — > e + e~), 
J/ip + 7T°(J/ip — > e + e~) or J/ip + ^(J/ip — > e + e~), all of which include the factor 
B(J/ip — > e + e~) = 0.06. The resulting 90% confidence upper limits for the three 
final states are plotted in various combinations in Figs. 4.20 - 4.27 for three choices 
of resonance width: T = 0.5 MeV, T = 1.0 MeV and T = 2.0 MeV. 

Figure 4.20 shows the upper limits associated with T = 0.5, 1.0 and 2.0 MeV 1 P\ 
resonances for the J/vp + X channel from the E835p data. As can be seen from the 
plot, meaningful upper limits can only be set in the region ~ 3525-3527 MeV, where 
the bulk of the data was taken. Outside these regions, the upper limits rise rapidly, 
as the fits are unconstrained by actual data (there being only one data point ~ 2 
MeV lower, and two points ~ 2 MeV higher). As may be expected, the smallest 
width (r = 0.5 MeV) tends to give the largest upper limits. 

Figures 4.21-4.23 show the upper limits for the J/ip + n channel separately 
for the three experiments E835p, E835, and E760, each for the three choices of 
resonance width. The background for the E835 scan was obtained from the data 
at 3524.6 MeV and 3529.1 MeV and that for E760 from data at 3522.6 MeV and 
3523.5 MeV. The upper limit curves in all of these plots have the same basic shape 
as in the E835p J/ip + X case, but the overall level of the upper limits is lower in 
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the E835p plot (Figure 4.21) than in either the E835 plot (Figure 4.22) or the E760 
plot (Figure 4.23). This is essentially due to the larger luminosity invested in the 
E835p scan (50.5 pb" 1 in E835p as compared to 38.9 pb _1 in E835 and 15.9 pb _1 
in E760). 

In Figures 4.24-4.26, we compare the upper limits for the J/vp + ir° channel from 
the three experiments for each of the three choices of the resonance width, 0.5 MeV 
in Figure 4.24, 1.0 MeV in Figure 4.25, and 2.0 MeV in Figure 4.26. These figures 
clearly demonstrate the superiority of the E835p scan over that of the E835 and 
E760 scans in establishing more stringent upper limits. Finally, Figure 4.27 shows 
the upper limits for the C-parity violating reaction pp ^ Pi — > J/ip^ for the E835p 
data. As in the previous plots, upper limits are shown for 0.5, 1.0, and 2.0 MeV 
resonance widths. 

In summary, the present measurements do not show evidence for the presence of 
a resonance in the mass region 3524.8 MeV < M R < 3527.3 in any of the reactions 
studied: pp — > (h c ) — > J/ip + X , pp — > (h c ) — > J/ip + ir°, and the C parity forbidden 
reaction pp — > (h c ) — > J/ip + 7. 

In terms of the 90% confidence upper limits, we find that for assumed 1 Pi total 
widths larger than 0.5 MeV, the product branching ratios are: B(pp — > h c ) xB(h c — > 
J/il> + X(J/il> -> e+e - )) < 2.3 x 10~ 7 for M( x Pi) in the region 3524.6 MeV - 3527.0 
MeV, and < 1.7 x 10" 7 for M( x Pi) in the smaller region 3526.2 ± 0.25 MeV. The 
corresponding limits for the l P\ — > J/ip + it decay channel are 1.8 x 10~ 7 and 
0.9 x 10~ 7 , and those for the l P\ — > J/ip + 7 decay channel are < 1.6 x 10 -7 and 
0.5 x 10~ 7 . These upper limits, along with upper limit values for B in x S out at the 
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E760 claimed l P\ energy and at < M(x.j) >, the centroid of the x states, are listed 
in Tables 4.9 and 4.10. 

In the E760 experiment it was stated that for the claimed observation of 1 P\ 
at 3526.2 ± 0.2 MeV the product branching ratio, B(pp -> h c ) x B(h c -> J/tp + 
7T°( J/ij; -> e+e")) = (2.3 ± 0.6) x 10~ 7 for T(h c ) = 0.5 MeV and (1.7 ± 0.4) x lO" 7 
for T(h c ) = 1.0 MeV. These correspond to 90% confidence upper limits of 3.4 x 10~ 7 
and 2.3 x 10 -7 . Our present limits are therefore nearly a factor two smaller in both 
cases. 
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Figure 4.20: 90% CL upper limits for B in B out (pp P x -> J/^ + X, J/ip -> e + e"), 
calculated using the events from the tight J/ip + X data selection. The three curves 
show the results for r tot ( 1 P 1 ) = 0.5, 1 and 2 MeV. The stack energies are marked by 
the short lines at the bottom of the plot. The first vertical dashed line represents 
the center of gravity of the x states, and the second line shows the 1 P± mass value 
claimed by E760. 
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Figure 4.21: 90% CL upper limits for B in B out {pp P 1 -> J/ip + vr , J/ip -> e+e - ) 
for data from E835p. The three curves show the results calculated for T( 1 P 1 ) = 0.5, 1 
and 2 MeV. The first vertical dashed line represents the center of gravity of the x 
states, and the second line shows the 1 P± mass value claimed by E760. 



177 



x 10 



CD 

CL 



o 



f 0.4 



0.3 



0.2 




• T = 0.5 MeV 

O T = 1 MeV 

* T = 2 MeV 



3523 



3524 



5525 



3526 



3527 



3528 3529 
Energy (MeV) 



Figure 4.22: 90% CL upper limits for B in B out {pp P 1 -> J/tj) + tt°, J/ip -> e+e") 
for data from E835. The three curves show the results calculated for T( 1 P 1 ) = 0.5, 1 
and 2 MeV. The first vertical dashed line represents the center of gravity of the x 
states, and the second line shows the 1 P± mass value claimed by E760. 
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Figure 4.23: 90% CL upper limits for B in B out {pp P x -> J /if) + tt°, J /if) -> e + e~) 
for data from E760. The three curves show the results calculated for T( 1 P 1 ) = 0.5, 1 
and 2 MeV. The first vertical dashed line represents the center of gravity of the 
X states, and the second line shows the 1 P 1 mass value claimed by E760. The 
horizontal line in the middle of the plot represents the E760 90% CL upper limit for 
Bi n B ou i. 
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Figure 4.24: 90% CL upper limits for B in B out (pp P 1 -> J/ijj + tt , J/ip -> e + e~) 
for T( 1 P 1 ) = 0.5 MeV. The three curves shown show the results from the data 
from E760, E835, and E835p. The first vertical dashed line represents the center 
of gravity of the x states, and the second line shows the 1 Pi mass value claimed by 
E760. 
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Figure 4.25: 90% CL upper limits for B in B out (pp P x -> J/ip + n°, J/ip -> e + e~) 
for r( 1 P 1 ) = 1 MeV. The three curves shown show the results from the data from 
E760, E835, and E835p. The first vertical dashed line represents the center of gravity 
of the x states, and the second line shows the 1 Pi mass value claimed by E760. The 
horizontal line in the middle of the plot represents the E760 90% CL upper limit for 
Bi n B ut- 
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Figure 4.26: 90% CL upper limits for B in B out (pp P 1 -> J/ip + tt°, J/ip -> e+e") 
for r( 1 P 1 ) = 2 MeV. The three curves shown show the results from the data from 
E760, E835, and E835p. The first vertical dashed line represents the center of gravity 
of the x states, and the second line shows the 1 P\ mass value claimed by E760. 
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Figure 4.27: 90% CL upper limits for B in B out (pp P x -> J/ip + 7, J/ip -> e+e~). 
The three curves show the results for T( 1 P 1 ) = 0.5, 1, and 2 MeV. The stack energies 
are shown at the bottom of the plot. The first vertical dashed line represents the 
center of gravity of the x states, and the second line shows the 1 P 1 mass value 
claimed by E760. 
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Channel 


bnergy 


1 7~~) TX 7"* 1 J-l 

P\ Width 






(MeVJ 


(MeVJ 


LL Upper Limits 


J/ip + X 


3525.27 


0.5 


1.6 x 10" 7 






1 


1.3 x 10~ 7 






2 


1.1 x lO" 7 


111 i "V 

J/ip + X 


3526.2 


0.5 


1.7 x 10 






1 


l.l x 10 ' 






2 


l.O x 10~ 7 


JH) + 7T° 


3525.27 


0.5 


1.4 x 10~ 7 






1 


l.l x 10~ 7 






2 


0.8 x 10~ 7 


J/l/> + 7T° 


3526.2 


0.5 


l.l x 10~ 7 






1 


0.8 x 10~ 7 






2 


0.7 x 10~ 7 


J/ip + 7 


3525.27 


0.5 


1.4 x 10~ 7 






1 


0.8 x 10~ 7 






2 


0.5 x 10~ 7 


J/^ + 7 


3526.2 


0.5 


0.7 x 10~ 7 






1 


0.4 x 10~ 7 






2 


0.3 x 10~ 7 



Table 4.9: B in B out 90% CL Upper Limits for the J/ip + X, J/ip + vr°, and J/ip + 7 
channels, fit for 0.5, 1, and 2 MeV width 1 Pi resonances at the centroid of the xj 
states and the E760 observed 1 Pi energy. 



184 





90% Confidence Upper Limits xlO 7 


Mass Region (MeV) 
(Width) 


B{pp — > h c ) x B(pp — > h c ) x B(pp — »■ h c ) x 
B(h c -> J/V> + 7T°) B(h c ^ J/i/j + X) B(h c ^ J/ip + 7) 


3526.2 ± 1.1 
(r = 2.0 MeV) 
(r = 1.0 MeV) 
(r = 0.5 MeV) 


< 0.9 < 1.3 < 0.5 

< 1.2 < 1.8 < 1.0 

< 1.8 < 2.3 < 1.6 


3526.2 ± 0.25 
(r = 2.0 MeV) 
(r = 1.0 MeV) 
(r = 0.5 MeV) 


< 0.7 < 1.0 < 0.3 

< 0.8 < 1.1 < 0.5 

< 1.1 < 1.7 < 0.8 



Table 4.10: 90% confidence upper limits from the present measurements for the 
mass regions 3526.2 ± 1.1 MeV and 3526.2 ± 0.25 MeV. 



Chapter 5 



The pp Pi 

—> Vcl —> 37 Reaction 

5.1 The 77 c 7 Decay Channel of /i c ( 1 Pi) 

An additional 1 P\ decay channel which was examined in the E835p data was the 
reaction pp — ^ Pi — > rj c j — > (77)7. This channel was examined even though it 
lacks a J/^ -* e + e~ signal for clear charmonium identification. This is because 
the radiative decay to r\ c is expected to be the dominant decay mode of the 1 P\ 
being the spin singlet version of the radiative decay Xci — > J/ipl for which the PDG 
gives a branching ratio of ~ 32%. Unfortunately, i] c (0~ + ) cannot decay into e + e~, 
but decays almost entirely hadronically. This makes its identification difficult in a 
detector optimized for the detection of electrons and photons. There is, however, 
a two photon decay of the r] c , with a branching ratio of ~ 3 x 10~ 4 which can be 
identified in the E835 detector. In fact, if the products of the estimated branching 
ratios for both the J/ipir° — > e + e~77 and 7^7 — > (77)7 decays of the 1 P 1 are 
compared, they should be of roughly the same order. This is because the J/ipn 
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decay is suppressed by approximately two orders of magnitude relative to 77 c 7, but 
B(r] c — > 77) = 3 x 10~ 4 is suppressed by roughly two orders of magnitude relative to 
B(J/ip — > e + e~) = 6 x 1CT 2 . Thus the expected decay rate of l P\ — > ^7 — > (77)7 
is comparable to 1 P i — > J/ipn° — > (e + e~)77. 

5.2 GEANT Monte Carlo Study of the rj c <y Decay 
Channel 

A study of the 7^7 channel was done using the full E835p GEANT Monte Carlo 
to examine the feasibility of detecting a signal in this channel. Figure 5.1 shows 
the two photon invariant mass spectrum from the GEANT Monte Carlo simulation 
of the reaction pp — s- 1 Pi — > 7^7 — > 37. The resulting 77 c signal shows a spread in 
-^(7172) over 2.98 ± 0.20 GeV due to detector effects. This range was eventually 
chosen as the r/ c mass cut for this channel. 

Figure 5.2 shows the 7173 invariant mass spectrum from the GEANT Monte 
Carlo generated events for the r/ c 7 — > (7172)73 decay. The Monte Carlo, of course, 
does not contain tt , 77, and 7/ which are produced in great numbers in pp annihila- 
tions, as shown later. To remove them from the data, cuts were made in the data 
for M(7i7 3 ) < 1.1 GeV and M(7 2 7 3 ) < 1.1 GeV. This is indicated by the dashed 
line in Fig. 5.2. The expected angular distribution of photons from the r] c decay 
as calculated by GEANT Monte Carlo simulation is shown in Figure 5.3, and as 
expected it is flat over the acceptance region. Although the angular distribution 
of the two gammas coming from the decay of the r\ c is expected to be essentially 
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isotropic, the gamma background in the detector is known to be strongly peaked in 
the forward and backward regions. This forward-backward peaking is known to arise 
from the very large photon yield from n° decays from pp — > 7r°7r° events. Cuts are 
made in the data to remove events with | cos(0*), cos(^| > 0.4 in order to eliminate 
the large background from feeddown from 7r°7r° events. The cuts are indicated by 
dashed lines in Fig. 5.3. 

Finally, the total energy of the GEANT Monte Carlo event in the laboratory 
frame, E^lab), for the events which survive the invariant mass and angular cuts 
is plotted in Figure 5.4 with respect to the center of mass energy of the radiative 
gamma (73) which is expected to be monoenergetic with E cm (^) = 504 MeV. These 
two quantities are shown again in projections in Figure 5.5, where the detector effects 
can be clearly seen: they spread the events around the nominal values of total energy 
(6.626 GeV), and E cm (^ 3 ) = 504 MeV for the assumed M( 1 P 1 ) = 3.5262 GeV. Cuts 
are therefore made in the data sample to remove events which do not have 6.376 
GeV < E lab (tot) < 7.000 GeV, and 454 MeV < E cm (^) < 580 MeV. These cuts are 
non-symmetric about the nominal values due to the 1% energy correction applied 
to each 7, which is described in chapter 3. 
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M(7i7 2 ) ~ ^7 Monte Carlo Events 
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Figure 5.1: The 7^2 invariant mass distribution for all 3 7 events in the GEANT 
monte carlo simulation of the 1 P\ — > ?y c 7 — > 37 decay channel. % mass cuts are 
made for events with invariant mass M(7i7 2 ) < 2.78 and M(7i7 2 ) > 3.18 GeV. 



189 



M(7-i7 3 ) ~ ^7 Monte Carlo Events 
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Figure 5.2: The 7x73 invariant mass distribution for all 3 7 events in the GEANT 
Monte Carlo simulation with 2.8 GeV < M (7172) < 3.2 GeV. 
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costf*(70 - 3y E835p GEANT MC Events 




Figure 5.3: Top: The cos(0*)(7i )2 ) distribution for all 3 7 events in the GEANT 
Monte Carlo simulation of the reaction pp — s- 1 Pi — > rj c j — > (7172)7- The vertical 
dashed lines show the cuts made to remove events with | cos(0*)(7i)| > 0.4. Bottom: 
The same plot made with a strictly kinematic Monte Carlo with 100 times the 
number of generated events. 
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E tot ( I a b) vs. E 3 (cm) - 3y E835p Monte Carlo Events 
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Figure 5.4: The i? tot (lab) vs. -E cm (73) distribution for all 3 7 events in the GEANT 
Monte Carlo with 2.78 GeV < M( lll2 ) < 3.18 GeV, and M( 7l 7 3 ), M( l2 lz) > 11 
GeV and | cos(0*)(7i, 72) | < 0.4. The box outlines the region in which 6.376 GeV 
< E(tot) < 7.000 GeV and 0.454 GeV < £(73) < 0.580 GeV. 
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E lol (lab) - 3y GEANT Monte Carlo Events 
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Figure 5.5: Top: The E tot (lab) distribution for all 3 7 events in the GEANT Monte 
Carlo with 2.78 GeV < M( 7l72 ) < 3.18 GeV, and M( 7l 7 3 ), M( 72 7 3 ) > 1.1 GeV 
and I cos(0*)(7i,7 2 )| < 0.4. Bottom: The E cm (^) distribution for the same GEANT 
Monte Carlo events. 
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5.3 Data and Analysis Cuts for the 7^7 Decay 
Channel 

The ?7 C 7 decay channel was studied for the 21 energy points of the 1 Pi search using 
the neutral trigger data collected in E835p. Most of the cuts in the data are based 
on GEANT Monte Carlo simulations as described in Section 5.2. Here we refer only 
to the figures illustrating the Monte Carlo results. All three gamma events from 
this data sample were tested for their highest possible two gamma invariant mass 
combinations. These gammas were designated as 71 and 72 . The remaining gamma 
was, of course, 73. Events with M(7i7 2 ) > 2.75 GeV were used in the data analysis. 
The 7i7 2 invariant mass spectrum for the entire data sample is shown in Figure 5.6. 
(see Figure 5.1 for the corresponding Monte Carlo results) Note that at this stage 
no enhancement is visible at the r\ c region around 2.98 GeV. 

Figure 5.7 shows the 7173 invariant mass spectrum. As can be clearly seen in the 
figure, a large fraction of these photons make a %° (M = 135 MeV) or an 77 (M = 547 
MeV). Figure 5.8 shows M(7 2 73) after a cut has been made for M(7i7 3 ) < 1.1 GeV 
to remove 7r°, 77, and rj . Now one can see rj (M = 958 MeV) and some other residual 
structure. The cut at ^^(71^73) < 1.1 GeV rejects all of these. Figure 5.9 shows 
M(7i7 2 ) after the M(7i7 3 ) and M(7 2 73) cuts have been made. No r\ c enhancement 
is visible, even now. 

Although the angular distribution of the two gammas coming from the decay of 
the rj c is expected to be isotropic (see Figure 5.3 for the Monte Carlo results), the 
gamma background in the detector is known to be strongly peaked in the forward and 
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backward regions because of mainly feed-down from the prolific yield of the reaction 
pp — > 7r°7r°. This can be seen in Figures 5.10 and 5.11, which show the angular 
distributions of 71 and 72 respectively. Cuts are made to remove background events 
with —0.4 < I cos(^*)(7i, 7 2 )| < 0.4. The invariant mass spectrum M(7i7 2 ) after 
the cuts are made to cos(6 l *)(7i, 72) is shown in Figure 5.12. 

Lastly, the total energy of the event in the lab frame for the remaining events in 
the E835p data is plotted in Figure 5.13 with respect to the center of mass energy of 
the radiative gamma (73). Cuts are made in the data sample to remove events which 
do not have 6.376 GeV < Ei a b(tot) < 7.000 GeV, and further cuts on the radiative 
gamma are made to remove events which do not have 454 MeV < E cm (^) < 580 
MeV. The levels of these cuts and their efficiencies are determined from the GEANT 
Monte Carlo (see Figure 5.4 for Monte Carlo results). 
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Figure 5.6: The 7172 invariant mass distribution for 3 7 events in the E835p 1 P\ 
stacks with 2.8 GeV < M 7l72 < 3.2 GeV. No enhancement is seen in the r] c mass 
region. 7 1; 72, and 73 are numbered according to their energies in the center of mass 
frame, with ji being the highest energy and 73 the lowest. Using this numbering 
system for the reaction pp — ^ Pi — > r] c j — > 37, the two 7s from the r\ c will be 71 
and 72 while the radiative 7 will be 73. 
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1^1(7,73) — 3y E835p Data Events 




V(y,y,} (CcV) 

Figure 5.7: The 7173 invariant mass distribution for 3 7 events in the E835p 1 P\ 
data. The 7r° and 7/ peaks are clearly visible. Mass cuts are made to remove events 
with invariant mass M(7i7 3 ) < 1.1 GeV. 
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Figure 5.8: The 7273 invariant mass distribution for 3 7 events in the E835p 1 P\ 
data. The vertical axis is plotted logarithmically so that the contamination from 
the rj' is visible as well as the tc° and 77 peaks. Mass cuts are made to remove events 
with invariant mass M(7 2 7 3 ) < 1.1 GeV. 
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I(7i7 2 ) ~~ 3y E835p Data Events 



3073 Events 




2.75 2.8 2.85 2.9 2.95 



3.05 3.1 3.15 



3.2 3.25 
U(y,y 2 ) (GeV) 



Figure 5.9: The 7172 invariant mass distribution for 3 7 events in the E835p 1 P\ 
data with 2.78 GeV < M( 7 i 72 ) < 3.18 GeV, and M( 7l 7 3 ), M( 72 7 3 ) > 1.1 GeV. As 
before, no r\ c enhancement is visible at 2.98 GeV. 
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Figure 5.10: The cos(0*)(7i) distribution for 3 7 events in the E835p 1 P\ data with 
2.78 GeV < M( 7 i7 2 ) < 3.18 GeV, and M( 7 i7 3 ), M(7 2 7 3 ) > 1.1 GeV. As shown in 
Fig. 5.3, the distribution is expected to be flat within the acceptance range of the 
detector. A cut is made to remove background events with | cos(0*)(7i)| > 0.4. 
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Figure 5.11: The cos(#*)(7 2 ) distribution for 3 7 events in the E835p l P\ data 
with 2.78 GeV < M( 7l72 ) < 3.18 GeV, and M( 7l73 ), M( 72 7 3 ) > 1.1 GeV and 
I cos(#*)(7i)| < 0.4. As shown in Fig. 5.3, the distribution is expected to be flat 
within the acceptance range of the detector. A cut is made to remove background 
events with | cos(0*)( 7 i)| > 0.4. 
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Figure 5.12: The M(7x7 2 ) distribution for 3 7 events in the E835p X P\ data 
with 2.78 GeV < M( 7l72 ) < 3.18 GeV, and M( 7l 7 3 ), M( 72 7 3 ) > 1.1 GeV and 
|cos(0*)(7i, 72 )| < 0.4. 
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E tot (lab) vs. E 3 (cm) - 3y E835p Data Events 
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Figure 5.13: The £' tot (lab) vs. E cm (^j^) distribution for 3 7 events in the E835p 
1 P 1 data with 2.78 GeV < M{ lll2 ) < 3.18 GeV, and M( 7l73 ), M( 72 7 3 ) > 1.1 GeV 
and I cos(0*)(7i, 72)1 < 0.4. E(tot) and £"(73) cuts have not been made yet. The 
box outlines the region in which 6.376 GeV < E(tot) < 6.876 GeV and 0.454 GeV 
< £(73) < 0.554 GeV. 
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5.4 rjcj Final Event Selection 

As a final check on the consistancy of the selected events the scatter plot of the 
7i7 2 invariant mass vs. the center of mass energy of the radiative gamma (73) is 
plotted in Figure 5.14. 38 events are left in the final selection, and their properties 
are given in Tables 5.1 - 5.4. The efficiencies of the various cuts are given in Table 
5.5, and their excitation curve is shown in Figure 5.15. The average cross section 
for the reaction pp — > r} c j is measured to be 16.8 ± 2.7 pb in the 3523-3529 MeV 
region. This is compared to the excitation curve from the 3 7 data in the xo energy 
region generated using the same set of cuts, which is shown in Figure 5.16. with an 
average cross section of 11.3 ± 3.4 pb. 
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M(y 1 7 2 ) vs. E 3 (cm) — 3y E835p Data Events 
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Figure 5.14: The M^i^) vs. -E cm (73) distribution for all 3 7 events in the E835p 
1 Pi data in the final event selection, including cuts on E tot (lab) and E cm {j^) as 
shown in Figs 4.38 and 4.39. 
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Table 5.1: Properties of the events in the final selection from the E835p l P\ — > 
Vcl — > 37 data channel (Part I). The center of mass energy of the third gamma is 
given as are the three possible invariant masses. 
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Table 5.2: Properties of the events in the final selection from the E835p l P\ — > 
Vcl — ► 37 data channel (Part II). The center of mass energy of the third gamma is 
given as are the three possible invariant masses. 
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Table 5.3: Properties of the events in the final selection from the E835p l P\ — > 
Vcl — ► 37 data channel (Part I). The angles of 71 and 72 are given, as well as the 
total event energy in the lab frame, and the center of mass energy of 73. 
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28 


5811 


134846 


-0.248 


-0.237 


6.649 


0.580 


28 


5812 


2006974 


-0.145 


-0.031 


6.552 


0.522 


35 


7005 


673075 


0.362 


-0.342 


6.484 


0.563 


35 


7005 


2551591 


0.258 


-0.330 


6.950 


0.502 


37 


7021 


2417964 


-0.008 


0.237 


6.953 


0.567 


38 


7025 


3092852 


-0.045 


0.389 


6.432 


0.482 


38 


7028 


318008 


0.267 


0.302 


6.741 


0.559 


39 


7031 


2763321 


0.277 


-0.363 


6.665 


0.532 


39 


7036 


2967823 


0.282 


0.067 


6.616 


0.506 


Expected 










6.626 


0.504 



Table 5.4: Properties of the events in the final selection from the E835p l P\ — > 
Vcl — * 37 data channel (Part II). The angles of 71 and 72 are given, as well as the 
total event energy in the lab frame, and the center of mass energy of 73. 



Cut 


Efficiency 


Trigger 
3 Cluster Cut 


90% ± 1% 
94% ± 1% 


Analysis Cuts 
2780 MeV < M( lll2 ) < 3180 MeV 
M(7 1 7 3 ),M(7 2 7 3 ) < 1100 MeV 
|cos(0J j2 )| < 0.4 
6376 MeV < E tot (\ah)< 7000 GeV 
454 MeV < E y3 (cm) < 580 MeV 
Total Analysis Cuts 


76.9% ± 0.8% 
37.8% ± 0.3% 
51.5% ±0.4% 
91.8% ±4.0% 
88.8% ± 4.0% 
(12.2 ±0.8)% 


Acceptance 


(48.7 ±0.1)% 


Total 


(5.0 ±0.3)% 



Table 5.5: Efficiencies of all cuts used in rj^ channel. 
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Figure 5.15: Efficiency corrected excitation curve for the r/ c + 7 channel. Stacks with 
energies less than 0.15 MeV apart have been combined. Numbers beside each point 
denote number of observed events at that energy. The straight line is the fit to a 
constant cross section. 
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Figure 5.16: Efficiency corrected excitation curve for the 7] c + 7 channel in the xo 
energy region. The straight line is the fit to a constant cross section. 
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Appendix A 

pp — ^ P\ — > + 7T° — > e + e _ 77 
Selected Events 

Tables A. 1-5 list the properties of all 107 events included in the final event selection 
for the channel pp ^ Pi — > J/ip + 7r° — > e + e _ 77, as described in Chapter 4. 
Properties included in the tables are the stack, run, and event numbers of each 
event, the two electron weights (EW(1), EW(2)), the invariant mass of the e + e~ 
pair M(e + e~), the invariant mass of the 77 pair M(77) which make a 7r°, the 
J/ip + 7T° probability as calculated by a 6 constraint kinematic fit, and the width of 
the beam for each run. 
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Table A.l: Properties of the events in the final selection from the E835p l P\ 
J/ipn° — > e + e _ 77 data channel (Part 1). 
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Table A. 2: Properties of the events in the final selection from the E835p l P\ 
J/ipn — > e + e~77 data channel (Part 2). 
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Table A. 3: Properties of the events in the final selection from the E835p l P\ 
J/ipn° — > e + e~77 data channel (Part 3). 
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Table A. 4: Properties of the events in the final selection from the E835p l P\ 
J/ipn — > e + e~77 data channel (Part 4). 
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Table A. 5: Properties of the events in the final selection from the E835p l P\ 
J/ipn — > e + e _ 77 data channel (Part 5). 



Appendix B 

Luminosity Monitor Software 

The following is the FORTRAN code used to read data from the E835 luminosity 
monitor and calculate the luminosities for each data run: 

PROGRAM LUMREAD 
C= READLUM.F 

C= ADC Channel 1: fixed detector 1 (FBL) 
C= ADC Channel 2: fixed detector 2 (FBR) 
C= ADC Channel 3: unused 
C= ADC Channel 4: movable detector (4 degree) 
C= Description: 

C= 1) Open Luminosity Data File, Read the Spectrum From File by DAFT. 

C= 2) Fill HBOOK Histograms (for Spectrums) 

C= 3) Print out the Integrated Luminosity for input run 

C= 4) Print out the Instantaneous Luminosity for input run 

C= 5) Print out the Elapsed Time and Live Time for input run 

C= Modified For 
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C= 



C= 



c= 



c= 



c= 



1) DAQ Spikes 

2) DAQ hangs 

3) Background calculation 

4) Cross Section calculation 

5) Add movable detector for the beam displacement correction 



IMPLICIT NONE 
include 'daft.inc' 

INTEGER NEVT,I,J,K,IEVENT, ii 
CHARACTER*80 FILEN 
CHARACTER*4 RUNN 
INTEGER BUFFERSIZE, STATUS 
PARAMETER (BUFFERSIZE = 32768) 
INTEGER BUFFER(BUFFERSIZE) 
INTEGER FDS(DAFTFDSIZE) 
COMMON /E835DAFT / fds 
INTEGER*2 LENRC 

INTEGER EVENTSIZE, ICYCLE, ISTAT, TIMESIZE 
PARAMETER (EVENTSIZE = 512) 

REAL DATARAW(4, EVENTSIZE), DATACOR(4, EVENTSIZE) 
INTEGER DATALUM (EVENTSIZE) 

REAL RDATALUM(EVENTSIZE),CDATALUM(EVENTSIZE) 

REAL RCOUNTS(4), RCOUNTPS(4), CCOUNTS(4), CCOUNTPS(4) 

REAL RBGCNT(4,42),CBGCNT(4,42), BGCNL(4,42) 

REAL RCOUNTBG(4), CCOUNTBG(4) 

REAL RCOUNTT(4), RCOUNTOT(4) 

DATA RCOUNTS/0.0, 0.0, 0.0, 0.0/ 

DATA RCOUNTOT/0.0, 0.0, 0.0, 0.0/ 
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DATA RCOUNTBG/0.0, 0.0, 0.0, 0.0/ 

DATA CCOUNTS/0.0, 0.0, 0.0, 0.0/ 

DATA CCOUNTBG/0.0, 0.0, 0.0, 0.0/ 

INTEGER STARTTIME(6), STOPTIME(6) 

INTEGER BL(4), BR(4), INDEX 

INTEGER ELTIME, LITIME, COUNTS 

REAL XI,YI,XIYI,XI2,BGSUM,AA,BB 

REAL LDDTEL, LDDTLI, RRIO, RELTIME, RLITIME 

DATA LDDTEL/0.0/ 

DATA LDDTLI/0.0/ 

LOGICAL DOFIT 

REAL INSLUMFIXL,INSLUMFIXR,INSLUMMOV 
REAL INTLUMFIXL,INTLUMFIXR,INTLUMMOV 
REAL RINTLUMFIXL,RINTLUMFIXR,RINTLUMMOV 
REAL INSLUM, INTLUM 
C= Beam Displacement 

REAL BEAMXM(451), ASYMM(451), BEAMXP(451), ASYMP(451) 
REAL LUMDIFF, LUMTRUET, LUMTRUES, X, OFFSET 
REAL CFR, CFL, CFM 

REAL INTLUMTRUEL,INTLUMTRUER,INTLUMTRUEM 
REAL INSLUMTRUEL,INSLUMTRUER,INSLUMTRUEM 
REAL BMMIN ,BMMAX, BPMIN, BPMAX 
DATA BMMIN/2.0774424E-02/ ! beam displacement -1.5 mm 
DATA BMMAX/1.016171/ ! beam displacement -5.0 mm 
DATA BPMIN/-1.7685434E-02/ ! beam displacement 1.5 mm 
DATA BPMAX/-1. 017238/ ! beam displacement 5.0 mm 
C= Cross section calculation 
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REAL PI, MP, ALPHA, HBC2, DETRAD, B, RHO, ST, PLABFILE 

INTEGER RUNNUM, RUNFILE 

DATA PI/3.141592654/ 

DATA MP/0.93827231/ 

DATA ALPHA/0.00729927/ 

DATA HBC2/0.38937966/ 

REAL AST, BST, NST, BBA, BBB 

DATA AST/34.48/ ! const for SigmaT calculation 

DATA BST/89.7/ 

DATA NST/-0.702/ 

DATA BBA/13.643/ ! const for b calculation 
DATA BBB/-0.2053/ 

REAL DETANGL, DETANGR, DETANGM, DETANG 
DATA DETANGL/3.496/ ! in degree, left detectors 
DATA DETANGR/3.511/ ! in degree, left detectors 
DATA DETANGM/ 1.704/ ! in degree, movable detectors 
REAL OMEGAL, OMEGAR, OMEGAM, OMEGA 
DATA OMEGAL/0.91469e-4/ ! solid angle, left detectors 
DATA OMEGAR/0.91377e-4/ ! solid angle, right detectors 
DATA OMEGAM/1.4572e-4/ ! solid angle, movable detectors 
REAL LUMCONSTL, LUMCONSTR, LUMCONSTM, LUMCONST, RATIO 
REAL SIGMAC, SIGMAN, SIGMAI, SIGMAT, DTDO 
REAL PLAB, ELAB, BETA, RS, ECM, KA, T, CGT, GT, GT4, CC 
REAL AA1, AA2, A A3, BB1, BB2, BB3 
real vaa(4), vbb(4) 
C** histogram stuff 

CHARACTER*40 HINAME 
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INTEGER HINO 

REAL ELTIME2, INSLUM1, INSLUM2, INSLUM4, INSLUMI, LUMEST 
REAL HMEMOR 

COMMON/PAWC/HMEMOR(1500000) 
C= Program starts here 

C= Open lookup data file for beam displacement correction 

OPEN(unit=2, name='ASYMMETRY.TABLEM', type='old') 
DO I = 1, 451 

READ(2,*) BEAMXM(I), ASYMM(I) 
END DO 
CLOSE(2) 

OPEN(unit=2, name= ' AS YMMETRY.TABLEP ' , type='old') 

DO I = 1, 451 

END DO 

CLOSE(2) 
C= Cross section calculation 
C== Input the RUN NUMBER first 

TYPE *, 'Enter The RUN Number: ' 

Accept *, RUNNUM 

OPEN(UNIT=2, NAME='../dat/runplab.txt\ TYPE='old') 
1 READ(2, *, END=10) RUNFILE, PLABFILE 
IF(RUNNUM.EQ.RUNFILE) THEN 

PLAB = PLABFILE 

CLOSE(2) 

GO TO 20 
END IF 
GO TO 1 
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10 CLOSE(2) 

TYPE *, 'This RUN can not be found in the runplab.txt file' 
TYPE *, 'Please input the Pbar momentum in GeV: ' 
Accept *, PLAB 
20 DO 60 I = 1, 3 
IF (I.EQ.l) THEN 

DETANG = DETANGL 

OMEGA = OMEGAL 
ELSE IF(I.EQ.2) THEN 

DETANG = DETANGR 

OMEGA = OMEGAR 
ELSE IF(I.EQ.3) THEN 

DETANG = DETANGM 

OMEGA = OMEGAM 
END IF 

DETRAD = DETANG*PI/180.0 ! Angle in radian 
B = BBA + BBB*PLAB 
RHO = 0.0 

ST = AST + BST*PLAB**NST 
ELAB = SQRT(PLAB**2 + MP**2) 
BETA = PLAB/ELAB 

RS = (SQRT(PLAB**2+MP**2)+MP)**2 - PLAB**2 
ECM = SQRT(RS) ! Ecm calculation 
KA = (ELAB+MP) / (ELAB-MP) 
C- t value evaluation 

T = 2.0*MP*DETRAD*DETRAD/KA ! rough estimation 
t = 2.0*MP*T 
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C- Calculate dS/dt, first, Sigma-c 
cgt = (1.0+(abs(t))/0.71) 
g t = 1.0/cgt/cgt 
gt4 = gt*gt*gt*gt 

cc = 4.0*PI*alpha*alpha/beta/beta/t/t 

sigmac = cc*hbc2*gt4 
C- Calculate dS/dt, second, Sigma-n 

aal = (1.0+rho*rho)*st*st*exp(-b*abs(t)) 

aa2 = 16.0*PI*hbc2 

sigman = aal/aa2 
C- Calculate dS/dt, third, Sigma-i 

aal = alpha*st*gt*gt*exp(-b*abs(t)/2.0) 

aa2 = beta*abs(t) 

bbl = alog(b*abs(t)/2.0) 

bb2 = alog(1.0 + 8.0/0.71/b) 

cc = 4.0*abs(t)/0.71 

bb3 = cc*alog(cc) + cc/2.0 

aa3 = alpha*(0.577+bbl+bb2+bb3) 

sigmai = aal*(rho*cos(aa3)+sin(aa3))/aa2 

sigmat = sigmac + sigman + sigmai ! total sigma 
C- calculate dt/dOmega 

aal = (t + 4.0*mp*mp) 

aa2 = t*aal*aal*aal/ka 

dtdo = sqrt(aa2)/4.0/PI/mp/mp 
C- get the lumconst 

lumconst = sigmat * dtdo * omega 

IF(I.EQ.l) THEN 



240 

LUMCONSTL = LUMCONST 
ELSE IF(I.EQ.2) THEN 

LUMCONSTR = LUMCONST 
ELSE IF(I.EQ.3) THEN 

LUMCONSTM = LUMCONST 
END IF 
60 CONTINUE 

C= Boundaries for three detectors 
BL(1) = 75 ! left fixed detector 
BR(I) = 165 

BL(2) = 75 ! right fixed detector 
BR(2) = 165 

BL(4) = 200 ! movable detector 
BR(4) = 320 
BL(3) = 150 ! unused 
BR(3) = 335 
C== Do some initializations 
DO I = 1, 4 
DO J = 1, 21 

BGCNL(I,J) = FLOAT(BL(I) - 11+J) 
END DO 
DO J = 1, 21 

BGCNL(I,J+21) = FLOAT(BR(I) - 11+J) 
END DO 
DO J = 1, 42 

RBGCNT(I,J) = 0.0 
CBGCNT(I,J) = 0.0 
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ENDDO 
ENDDO 
DO I = 1, 512 

RDATALUM(I) = 0.0 

CDATALUM(I) = 0.0 

DATALUM(I) = 
END DO 
DO I = 1, 4 

DO J = 1, 512 

DATARAW(I,J) = 0.0 
DATACOR(I,J) = 0.0 

END DO 
END DO 
C== HBOOK stuff here 
CALL HLIMIT(1500000) 

call hbookl(101,'left detector $',512,0.0,512.0,0.) 
call hbookl(102,'right detector $',512,0.0,512.0,0.) 
call hbookl( 103, 'unused $',512,0.0,512.0,0.) 
call hbookl( 104, 'movable detector $',512, 
+ 0.0,512.0,0.) 

call hbookl(105, 'Ratio = luminosity L / R $',100,0.0,100.0,0.) 
call hbookl( 106, 'Ratio = lumonosity L / M $',100,0.0,100.0,0.) 
C= Input the name of luminosity data file 
IF(RUNNUM.LE.352) THEN 
WRITE(FILEN, 70) RUNNUM 
70 FORMAT('/lumin/runs00010352/lum000',I3,'.dat') 
RUNN = FILEN(27:30) 
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ELSE IF(RUNNUM.GT.352.AND.RUNNUM.LE.605) THEN 
WRITE(FILEN, 71) RUNNUM 

71 FORMAT('/lumin/runs03530605/lum000',I3,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GT.605.AND.RUNNUM.LE.786) THEN 
WRITE(FILEN, 72) RUNNUM 

72 FORMAT('/lumin/runs06060786/lum000',I3,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GT.786.AND.RUNNUM.LE.835) THEN 
WRITE (FILEN, 73) RUNNUM 

73 FORMAT('/lumin/runs07910835/lum000',I3,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GT.835.AND.RUNNUM.LE.999) THEN 
WRITE (FILEN, 74) RUNNUM 

74 FORMAT('/lumin/runs08360999/lum000',I3,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GT.999.AND.RUNNUM.LE.1999) THEN 
WRITE(FILEN, 75) RUNNUM 

75 FORMAT('/lumin/runsl0001999/lum00',I4,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GE.2000.AND.RUNNUM.LE.2999) THEN 
WRITE (FILEN, 76) RUNNUM 

76 FORMAT('/lumin/runs20002999/lum00',I4,'.dat') 

RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GE.3000.AND.RUNNUM.LE.3999) THEN 
WRITE (FILEN, 77) RUNNUM 

77 FORMAT('/lumin/runs30003999/lum00',I4,'.dat') 
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RUNN = FILEN(27:30) 
ELSE IF(RUNNUM.GE.4000) THEN 
WRITE(FILEN, 78) RUNNUM 
78 FORMAT('../dat/lum00',I4,'.dat') 
RUNN = FILEN(13:16) 
END IF 
C== Open file 

STATUS = daftopenf(fds, filen, DAFTRDONLY, buffer, BUFFERSIZE) 
IF(STATUS.NE.O) THEN 

TYPE *, 'Error Opening File, Status = ', STATUS 

STOP 
ENDIF 

NEVT = 1000000000 
C= Read the start time 
TIMESIZE = 6 

LENRC = daftgeteventf(fds, STARTTIME, TIMESIZE, 

& DAFTBORSKIP.or.DAFTEORSKIP.OR.DAFTLABSKIP.or. 

& DAFTBOFSKIP.or.DAFTEOFSKIP) 

IF (STARTTIME(l) .GT. 99) STARTTIME(l) = 

IF(LENRC.LT.0)THEN 

TYPE *, 'Error During First Read, Status = ', LENRC 

STOP 
ENDIF 
C= Read the spectrum 

DO 200 IEVENT = 1, NEVT 
DO 100 I = 1,4 

LENRC = daftgeteventf(fds, datalum, EVENTSIZE, 
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& DAFTBORSKIP.or.DAFTEORSKIP.OR.DAFTLABSKIP.or. 
k DAFTBOFSKIP.or.DAFTEOFSKIP) 
IF(LENRC.EQ.6) THEN 
C= Read the stop time 
DO J = 1, 6 

STOPTIME(J) = DATALUM(J) 
ENDDO 

IF (STOPTIME(l) .GT. 99) STOPTIME(l) = 
TYPE *, 'You have reached the end of this file' 
status = daftclosef(fds) 
GOTO 300 
ENDIF 

IF(LENRC.LT.O)THEN 

TYPE *, 'Error During Reading Spectrum, Status = ', LENRC 

goto 300 
ENDIF 

DO J = 1, EVENTSIZE 

DATARAW(I,J) = float (DATALUM(J) ) 

DATACOR(I,J) = float(DATALUM(J)) 
END DO 
100 CONTINUE 

C= Read the ELTIME and LITIME 
TIMESIZE = 1 

LENRC = daftgeteventf(fds, ELTIME, TIMESIZE, 

& DAFTBORSKIP.or.DAFTEORSKIP.OR.DAFTLABSKIP.or. 

& DAFTBOFSKIP.or.DAFTEOFSKIP) 

LENRC = daftgeteventf(fds, LITIME, TIMESIZE, 
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& DAFTBORSKIP.or.DAFTEORSKIP.OR.DAFTLABSKIP.or. 
& DAFTBOFSKIP.or.DAFTEOFSKIP) 
IF(LENRC.LT.O)THEN 

TYPE *, 'Error During Time Reading , Status = ', LENRC 

STOP 
ENDIF 

C== Correction by live time for each spectrum 
DO I = 1, 4 

DO J = 1, EVENTSIZE 

DATACOR(I,J) = DATACOR(I,J)*FLOAT(LITIME)/FLOAT(ELTIME) 

END DO 
END DO 
DO 150 I = 1, 4 
DO J = 1, EVENTSIZE 

RDATALUM(J) = DATARAW(I,J) 

CDATALUM(J) = DATACOR(I,J) 
END DO 

C — book individual 120 second interval histogram 

WRITE (HIN AME , 777) IEVENT 
777 FORMAT ('intervall no. ',15) 

HINO = (IEVENT* 10 + I) + 1000 

CALL HBOOKl(HINO, HINAME, 512, 0., 512., 0.) 
C= Fill the spectrums 

DO 105 J = 1, EVENTSIZE 

CALL HF1(100+I,FLOAT(J),RDATALUM(J)) 
CALL HFILL(HINO, FLOAT (J), 0., RDATALUM(J)) 
105 CONTINUE 



246 

C= Sum the boundaries counts for background calculation. 
C- Raw counts 
DO J = 1, 21 

RBGCNT(I,J) = RBGCNT(I,J) + RDATALUM(BL(I)-11+J) 
END DO 
DO J = 1, 21 

RBGCNT(I,J+21) = RBGCNT(I,J+21) + RDATALUM(BR(I)-11+J) 
END DO 
C- Live time corrected counts 
DO J = 1, 21 

CBGCNT(I,J) = CBGCNT(I,J) + CDATALUM(BL(I)-11+J) 
END DO 
DO J = 1, 21 

CBGCNT(I,21+J) = CBGCNT(I,21+J) + CDATALUM(BR(I)-11+J) 
END DO 
C= Sum the peak counts 
C- Raw counts 

RCOUNTPS(I) = 0.0 
DO J = BL(I), BR(I) 

RCOUNTPS(I) = RCOUNTPS(I) + RDATALUM(J) 
END DO 

RCOUNTS(I) = RCOUNTS(I) + RCOUNTPS(I) 
RCOUNTT(I) = 0.0 
DO J = 1, 512 

RCOUNTT(I) = RCOUNTT(I) + RDATALUM(J) 
END DO 

RCOUNTOT(I) = RCOUNTOT(I) + RCOUNTT(I) 
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C- Live time corrected counts 
CCOUNTPS(I) = 0.0 
DO J = BL(I), BR(I) 

CCOUNTPS(I) = CCOUNTPS(I) + CDATALUM(J) 
END DO 

CCOUNTS(I) = CCOUNTS(I) + CCOUNTPS(I) 
150 CONTINUE 

if ((rcountps(2).gt.O) .and. (rcountps(4).gt.O)) then 
ratio = rcountps(l)/rcountps(2) 
call hfl(105,float(ievent), ratio) 
if (rcountps(4) .gt. 0.) then 

ratio = rcountps(l)*lumconstm/rcountps(4)/lumconstl 
call hfl(106,float(ievent),ratio) 
endif 
endif 

LDDTEL = LDDTEL + 1.0E-6*FLOAT(ELTIME) 
LDDTLI = LDDTLI + 1.0E-6*FLOAT(LITIME) 
ELTIME2 = 1.0E-6*FLOAT(ELTIME) 

INSLUM1 = RCOUNTPS(l)*1.0E+27/LUMCONSTL/ELTIME2 
INSLUM2 = RCOUNTPS(2)*1.0E+27/LUMCONSTR/ELTIME2 
INSLUM4 = RCOUNTPS(4)*1.0E+27/LUMCONSTM/ELTIME2 
C - in case movable detector is 

INSLUMI = (INSLUMl+INSLUM2)/2.0 

write(6, *) IEVENT, 1.0E-6*FLOAT(ELTIME), 1.0e-6*float(litime), 
$ INSLUMI 
200 CONTINUE 

C== Determine the background calculation method 
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C== Raw counts 
300 DOFIT = .TRUE. 
DO 305 I = 1,4 

IF(I.EQ.3) GOTO 305 
DO J = 1, 42 

IF(RBGCNT(I,J).LT.1.0) DOFIT = .FALSE. 
ENDDO 
305 CONTINUE 
C== Background calculation 
IF(DOFIT) THEN 
DO 306 I = 1,4 
IF(I.EQ.3) GOTO 306 
XI = 0.0 
YI = 0.0 
XIYI = 0.0 
XI2 = 0.0 
DO J = 1, 42 

RBGCNT(I,J) = ALOG(RBGCNT(I,J)) 
YI = YI + RBGCNT(I,J) 
XI = XI + BGCNL(I,J) 
XIYI = XIYI + RBGCNT(I,J)*BGCNL(I,J) 
XI2 = XI2 + BGCNL(I,J)*BGCNL(I,J) 
END DO 

AA = (XPXIYI - YI*XI2)/(XI**2 - 42.0*XI2) 
BB = (42.0*XIYI - YI*XI)/(XI**2 - 42.0*XI2) 
A A = EXP(AA) 
type *, 'aa =', aa, 'bb =', bb 
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vaa(i) = aa 

vbb(i) = bb 

DO K = BL(I), BR(I) 

RCOUNTBG(I) = RCOUNTBG(I) + AA*EXP(-BB*FLOAT(K)) 
ENDDO 
306 CONTINUE 
ELSE 

DO 308 I = 1,4 

IF(I.EQ.3) GOTO 308 
BGSUM = 0.0 
DO J = 1, 42 

BGSUM = BGSUM + RBGCNT(I,J) 
ENDDO 

RCOUNTBG(I) = BGSUM*FLOAT(BR(I)-BL(I)+1)/42.0 
308 CONTINUE 
END IF 

C== Live time corrected counts 
DOFIT = .TRUE. 
DO 405 I = 1,4 
IF(I.EQ.3) GOTO 405 
DO J = 1, 42 

IF(CBGCNT(I,J).LT.1.0) DOFIT = .FALSE. 
ENDDO 
405 CONTINUE 

IF (DOFIT) THEN 
DO 406 I = 1,4 

IF(I.EQ.3) GOTO 406 
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XI = 0.0 
YI = 0.0 
XIYI = 0.0 
XI2 = 0.0 
DO J = 1, 42 

CBGCNT(I,J) = ALOG(CBGCNT(I,J)) 

YI = YI + CBGCNT(I,J) 

XI = XI + BGCNL(I,J) 

XIYI = XIYI + CBGCNT(I,J)*BGCNL(I,J) 

XI2 = XI2 + BGCNL(I,J)*BGCNL(I,J) 
END DO 

AA = (XPXIYI - YI*XI2)/(XI**2 - 42.0*XI2) 
BB = (42.0*XIYI - YI*XI)/(XI**2 - 42.0*XI2) 
AA = EXP(AA) 
DO K = BL(I), BR(I) 

CCOUNTBG(I) = CCOUNTBG(I) + AA*EXP(-BB*FLOAT(K)) 
ENDDO 
406 CONTINUE 
ELSE 

DO 408 I = 1,4 

IF(I.EQ.3) GOTO 408 

BGSUM = 0.0 

DO J = 1, 42 

BGSUM = BGSUM + CBGCNT(I,J) 
ENDDO 

CCOUNTBG(I) = BGSUM*FLOAT(BR(I)-BL(I)+1)/42.0 
408 CONTINUE 
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END IF 

C== ======= Background calculation done ======== 

C== Output the statistics 

TYPE * ' ' 

TYPE *, ' ' 

TYPE *, ' LUMINOSITY MONITOR OUTPUT' 
TYPE *, ' =========================' 

TYPE *, ' RUN NUMBER = ', RUNN 
WRITE(6, 414) ECM 

414 FORMAT(26X, 'Ecm = ', F5.3) 
TYPE *, ' ' 

WRITE(6,415) (STARTTIME(I), 1=1,3) 

415 F0RMAT(1X,'RUN BEGIN DATE: ',I2.2,'/',I2.2,'/',I2.2) 
WRITE(6,416) (STARTTIME(I), 1=4,6) 

416 F0RMAT(1X,'RUN BEGIN TIME: ',I2.2,':',I2.2,':',I2.2) 
WRITE(6,417) (STOPTIME(I), 1=1,3) 

417 F0RMAT(1X,'RUN STOP DATE: ',I2.2,'/',I2.2,'/',I2.2) 
WRITE(6,418) (STOPTIME(I), 1=4,6) 

418 F0RMAT(1X,'RUN STOP TIME: ',I2.2,':',I2.2,':',I2.2) 
TYPE *, ' ' 

TYPE *, 'ELAPSED TIME FOR THIS RUN, ET =', LDDTEL, V 

IF(LDDTLI.LE.l.O) LDDTLI = LDDTEL 

TYPE *, 'LIVE TIME FOR THIS RUN, LT =', LDDTLI, V 

type *, 

+' ' 

TYPE *, ' ' 
type *, 
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+'Detector Counts in Luminosity Effective lum.' 
type *, 

+' recoil peak (nb-1) corr for beam ' 
type *, 

+' offset, LT (nb-1)' 
type *, 



write(6, *) 
DO 500 I = 1,4 

IF (I.EQ.3) GOTO 500 

lumest = rcountot(i)*0.60*(lddtli/lddtel)*1.0E-6/lumconstl 
write (6, *) i, rcountot(i), rcounts(i), rcountbg(i), 
$ rcountbg(i)/rcounts(i) 

write (6, *) i, (rcounts(i)-rcountbg(i))/rcountot(i), 
$ lddtli/lddtel, lumest 

RCOUNTS(I) = RCOUNTS(I) - INT(RCOUNTBG(I)) 
CCOUNTS(I) = CCOUNTS(I) - INT(CCOUNTBG(I)) 
500 CONTINUE 

RINTLUMFIXL = RCOUNTS(l)*1.0E-6/LUMCONSTL 
RINTLUMFIXR = RCOUNTS(2)*1.0E-6/LUMCONSTR 
RINTLUMMOV = RCOUNTS(4)*1.0E-6/LUMCONSTM 
INTLUMFIXL = CCOUNTS(l)*1.0E-6/LUMCONSTL 
INTLUMFIXR = CCOUNTS(2)*1.0E-6/LUMCONSTR 
INTLUMMOV = CCOUNTS(4)*1.0E-6/LUMCONSTM 
INTLUM = (INTLUMFIXL+INTLUMFIXR+INTLUMMOV)/3.0 
INSLUMFIXL = CCOUNTS(l)*1.0E+27/LUMCONSTL/LDDTEL 
INSLUMFIXR = CCOUNTS(2)*1.0E+27/LUMCONSTR/LDDTEL 
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INSLUMMOV = CCOUNTS(4)*1.0E+27/LUMCONSTM/LDDTEL 
INSLUM = (INSLUMFIXL+INSLUMFIXR+INSLUMMOV)/3.0 
C== Beam Offset Correction 

LUMDIFF = 2.0*(INTLUMFIXL - INTLUMFIXR) 
+ / (INTLUMFIXL + INTLUMFIXR) 
X = LUMDIFF 

IF(LUMDIFF.GT.BPMIN. AND. LUMDIFF. LT.BMMIN) THEN 
INTLUMTRUEL = INTLUMFIXL 
INTLUMTRUER = INTLUMFIXR 
INTLUMTRUEM = INTLUMMOV 
LUMTRUET = INTLUM 
LUMTRUES = INSLUM 
ELSE IF(LUMDIFF.LT.BPMAX.OR.LUMDIFF.GT.BMMAX) THEN 
IF(LUMDIFF.LT.O.O) THEN 

CFL = 0.9972-1.1062*X-0.00507*X**2-1.0338*X**3 

CFR = 0.9997+0.00749*X+0.0500*X**2 

CFM = 1.0007-0.08216*X+0.2775*X**2 
ELSE 

CFL = 0.9997-0.0103*X+0.0539*X**2 
CFR = 0.9939+1.1045*X-0.01369*X**2+1.01459*X**3 
CFM = 1.0008+0.0771*X+0.2833*X**2 
END IF 

INTLUMTRUEL = CFL*INTLUMFIXL 
INTLUMTRUER = CFR*INTLUMFIXR 
INTLUMTRUEM = CFM*INTLUMMOV 

LUMTRUET = (INTLUMTRUEL+INTLUMTRUER+INTLUMTRUEM)/3.0 
INSLUMTRUEL = CFL*INSLUMFIXL 
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INSLUMTRUER = CFR*INSLUMFIXR 
INSLUMTRUEM = CFM*INSLUMMOV 

LUMTRUES = (INSLUMTRUEL+INSLUMTRUER+INSLUMTRUEM)/3.0 
ELSE 

IF(LUMDIFF.LT.O.O) THEN 
DO I = 50, 451 

IF(LUMDIFF.GE.ASYMP(I)) THEN 

OFFSET = BEAMXP(I) + (BEAMXP(I-l)-BEAMXP(I)) 
+ *(LUMDIFF-ASYMP(I))/(ASYMP(I-1)-ASYMP(I)) 
GOTO 585 
END IF 
END DO 
ELSE 

DO I = 50, 451 

IF(LUMDIFF.LE.ASYMM(I)) THEN 

OFFSET = BEAMXM(I) + (BEAMXM (I- 1 )-BE AMXM (I) ) 
+ *(LUMDIFF-ASYMM(I))/(ASYMM(I-1)-ASYMM(I)) 
GOTO 585 
END IF 
END DO 
END IF 

585 IF(LUMDIFF.LT.O.O) THEN 

CFL = 0.9972-1.1062*X-0.00507*X**2-1.0338*X**3 
CFR = 0.9997+0.00749*X+0.0500*X**2 
CFM = 1.0007-0.08216*X+0.2775*X**2 
ELSE 

CFL = 0.9997-0.0103*X+0.0539*X**2 
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CFR = 0.9939+1.1045*X-0.01369*X**2+1.01459*X**3 
CFM = 1.0008+0.0771*X+0.2833*X**2 
END IF 

INTLUMTRUEL = CFL*INTLUMFIXL 
INTLUMTRUER = CFR*INTLUMFIXR 
INTLUMTRUEM = CFM*INTLUMMOV 

LUMTRUET = (INTLUMTRUEL+INTLUMTRUER+INTLUMTRUEM)/3.0 
INSLUMTRUEL = CFL*INSLUMFIXL 
INSLUMTRUER = CFR*INSLUMFIXR 
INSLUMTRUEM = CFM*INSLUMMOV 

LUMTRUES = (INSLUMTRUEL+INSLUMTRUER+INSLUMTRUEM)/3.0 
END IF 
TYPE *, ' ' 

COUNTS = INT(RCOUNTS(l)) 

WRITE(6, 592) COUNTS, RINTLUMFIXL , INTLUMTRUEL 

592 FORMAT(lX,'Left ', 6X, 18, 9X, F9.3, 5X, F9.3) 
COUNTS = INT(RCOUNTS(2)) 

WRITE(6, 593) COUNTS, RINTLUMFIXR, INTLUMTRUER 

593 FO RM AT ( IX , ' Right ', 5X, 18, 9X, F9.3, 5X, F9.3) 
COUNTS = INT(RCOUNTS(4)) 

WRITE(6, 594) COUNTS, RINTLUMMOV, INTLUMTRUEM 

594 FORMAT (IX, 'Middle ', 4X, 18 ,9X, F9.3, 5X, F9.3) 
TYPE *, ' ' 

WRITE(6, 600) LUMTRUET 

600 FORMAT(8X, 'Average Integrated Luminosity = ', F9.3, ' nb-1') 
WRITE(6, 601) LUMTRUES 

601 FORMAT(8X, 'Average Instantaneous Luminosity = ',2X,E9.3) 
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type *, 

+' ' 

TYPE *, ' ' 

TYPE *, 'Asymmetry =', LUMDIFF 

IF(LUMDIFF.GT.BPMIN.AND.LUMDIFF.LT.BMMIN) THEN 

TYPE *, 'Beam offset = less than 1.5 mm' 
ELSE IF(LUMDIFF.LT.BPMAX.OR.LUMDIFF.GT.BMMAX) THEN 

TYPE *, 'Beam offset = greater than 5.0 mm' 

TYPE *, 'WARNING: CALL MCR !!!' 
ELSE 

IF(OFFSET.GT.O.O) THEN 
WRITE(6, 590) OFFSET 

590 FORMAT(lX,'Beam offset = +', F5.3, ' mm') 

ELSE 

WRITE(6, 591) OFFSET 

591 FORMAT(lX,'Beam offset = ', F6.3, ' mm') 

END IF 
END IF 
TYPE *, ' ' 

WRITE(6, 780) LUMCONSTL*1.0E06 

780 FORMAT(lX,'Luminosity constant, left detector = ', 
$ F12.6, ' nb') 

WRITE(6, 781) LUMCONSTR*1.0E06 

781 FORMAT(lX,'Luminosity constant, right detector = ', 
$ F12.6, ' nb') 

WRITE(6, 782) LUMCONSTM*1.0E06 

782 FORMAT(lX,'Luminosity constant, movable detector = ', $ F12.6, ' nb') 
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CALL HRPUT(0, 'luminst.hbook', 'NT') 

OPEN(10, file='lum.vec', form='formatted', type='unknown') 
vaa(3) = 1 
vbb(3) = 1 
do i = 1, 4 

write(10, *) vaa(i), vbb(i), bl(i), br(i) 
enddo 
close(lO) 
STOP 
END 



